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Chapter 



Preface 



Who cares about calculus? Well, calculus is used in virtually every branch 
of science where a problem can be mathematically modeled and an optimal 
solution is desired. In applications, quite often you "know" or strongly suspect 
something based on your understanding of the general principles in that field. 
Calculus allows you to quantify your expectations. 

For example, you know if you drop a car from a helicopter it will fall to 
the ground in a short length of time, possibly causing a lot of damage to the 
vehicle. Now suppose you are directing a truck commercial for television. You 
really want to know exactly how much damage is caused by dropping the car 
from a specific height so you can prepare precise instructions for the chopper 
pilot and cameramen. This quantitative precision is where mathematics, namely 
calculus, can help. 

Image a football quarterback throwing a football to a receiver during game. 
Now suppose you want to design a program for a computer game that allows 
players to play football on their computer or game console (or cell-phone) . You 
need to know exact equations for the football's trajectory in its path to the 
receiver. Again, this is where calculus plays a role. 

Hopefully this book will help provide a bit of foundation for you to go on to 
design the best video game or commercial or whatever your future brings to 
you! 

Acknowledgements: We thank Minh Van Nguyen for his careful proofreading 
and many excellent suggestions for revisions. 
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Chapter 1 

The Integral 



Differential calculus starts with a "simple" geometric problem: given a curve 
and a line tangent to that curve, find the slope of the line. The subject then 
develops techniques used in solving this problem into a combination of theory 
about derivatives and their properties, techniques for calculating derivatives, 
and applications of derivatives. This book begins the development of integral 
calculus and starts with the "simple" geometric idea of area. This idea will 
be developed into another combination of theory, techniques, and applications. 
The integral will be introduced in two (completely different) ways: as a limit 
of "Riemann sums" and as an "inverse" of differentiation (also known as "an- 
tiderivative" ) . Conceptually, the limit approach to the integral is geometric or 
numerical, while the antiderivative approach is somewhat more algebraic. 

One of the most important results in mathematics, The Fundamental Theo- 
rem of Calculus, appears in this chapter. It connects these two notions of the 
integral and also provides a relationship between differential and integral calcu- 
lus. Historically, this theorem marked the beginning of modern mathematics, 
and it provided important tools for the growth and development of the sciences. 
The chapter begins with a look at area, some geometric properties of areas, and 
some applications. First, we will see ways of approximating the areas of regions 
such as tree leaves that are bounded by curved edges, and the areas of regions 
bounded by graphs of functions. Then we will find ways to calculate exactly the 
areas of some of these regions. Finally, we will explore more of the rich variety 
of uses of "areas" . The primary purpose of this introductory section is to help 
develop your intuition about areas and your ability to reason using geometric 
arguments about area. This type of reasoning will appear often in the rest of 
this book and is very helpful for applying the ideas of calculus. 

1.1 Area 

We know from previous experience how to compute the areas of simple geo- 
metrical shapes like triangles, circles and rectangles. Formulas for these have 
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been known since the days of the ancient Greeks. But how do you find the area 
under a "more complicated" curve such as y = x 2 where — 1 < x < 1? First, 
let's graph the function using Sage as follows^: 

sage: a = -1; b = 1 

sage: f = lambda x: x~2 

sage: Lb = [[b,f(b)], [b,0] , [a,0] , [a,f(a)]] 

sage: Lf = [[i/20, f (i/20)] for i in xrange(20*a, 20*b+l)] 

sage: P = polygon (Lf+Lb, rgbcolor=(0.2,0.8,0)) 

sage: Q = plot(f(x), x, a-0.5, b+0.5) 

sage: show(P+Q) 

This results in the plot shown in Figure 1.11 



2.5 - 




■l - 



Figure 1.1: Using Sage to plot y = x 2 for —1 < x < 1. 

The rough, general idea introduced in this section is the following. To compute 
the area of the "complicated" shaded region in Figure 1.11 we break it up into 
lots of "simpler" subregions, whose areas are easy to compute, then add up the 
areas of the subregions to get the total area. We shall return to this example 
later. 

The basic shape we will use is the rectangle, whose area is given by the formula 

area = base x height 

If the units for each side of the rectangle are "meters," then the area will have 
the units ( "meters" ) x ( "meters" ) = "square meters" = m 2 . Other area formulas 
needed for this section are: the area formula for triangles 

area = Wi/2 (I- 1 ) 

^eel free to try this yourself, changing a, b and x 2 to something else if you like. 
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where b is the length of the base and h is the height; and the corresponding 
formula for circles 

area = irr 2 

with r being the radius. We also assume the following familiar properties of 
area: 

• Addition Property: The total area of a region is the sum of the areas of 
the non-overlapping pieces which comprise the region (see Figure 1.2) . 

• Inclusion Property: If region B is on or inside region A, then the area of 
region B is less than or equal to the area of region A (see Figure 1.3j ) . 

• Location-Independence Property: The area of a region does not depend 
on its location (see Figure [E4] ). 











\ iW 


+ 




+ 


! V L 

























Figure 1.2: Adding areas of non-overlapping regions. 
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Figure 1.3: Estimating areas using rectangles. 

Example 1.1.1. Determine the area of the region in Figure [LW a). 

Solution. The region can easily be broken into two rectangles, Figure 1.5( b), 
with areas 35 square inches and 3 square inches respectively. So the total area 
of the original region is 38 square inches. □ 

Using the above three properties of area, we can get information about areas 
that are difficult to calculate exactly. For instance, let A be the region bounded 
by the graph of f(x) = l/x, the x-axis, and vertical lines at x = 1 and x = 
3 (refer to Figure [L6]) . Since the two rectangles in Figure [L61 are inside the region 
A and do not overlap each other, the area of the rectangles is 1/2 + 1/3 = 5/6, 
which is less than the area of region A. 
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Each area i s 6 



Figure 1.4: Independence of area under translations and rotations. 
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Figure 1.5: Calculating areas using the addition property. 



Practice 1.1.1. Build two rectangles, each with base 1 unit, outside the shaded 
region in Figure \1.6j Then use their areas to make a valid statement about the 
area of region A. 
(Ans: Outside rectangular area = 1.5.) 

Practice 1.1.2. Now use both inside and outside rectangles with base 1/2 unit. 
What can be said about the area of region A in Figure \1.6\ ? 
(Ans: The area of the region is between 0.95 and 1.2.) 

Example 1.1.2. In Figure [7T7j , there are 32 dark squares, each having a side 
length of 1 centimeter, and 31 lighter squares of the same side length. We can 
be sure that the area of the leaf is smaller than what number? 

Solution. The area of the leaf is smaller than 32 + 31 = 63 cm 2 . □ 

Practice 1.1.3. We can be sure that the area of the leaf is at least how large? 
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1.1. AREA 




Figure 1.6: The area under of y = 1/x, where 1 < x < 3. 




Figure 1.7: Approximating the area of a "leaf". 



Functions can be defined in terms of areas. For the constant function f(t) = 2, 
define A(x) to be the area of the rectangular region bounded by the graph of 
/, the i-axis, and the vertical lines at t = 1 and t = x (see Figure 1.8( a)). 
Then A(2) is the area of the shaded region in Figure 1.8( b). where A(2) = 2. 
Similarly, A(S) = 4 and A(4) = 6. In general, we have 

A(x) = base x height = (x — 1) x 2 = 2x — 2 

for any x > 1. Figure rTTBT c) shows the graph of y = A(x), whose derivative is 
A'{x) = 2 for every value of x > 1. 




Figure 1.8: The area as a function. 



1.2. SOME APPLICATIONS OF AREA 



Sometimes it is useful to move regions around. Consider the parallelogram in 
Figure 1.9 and notice the triangle on its left end. Moving this triangular region 
from the left side of the parallelogram to fill the region on the right side, we 
end up with a rectangle. So we have reduced the problem of finding the area 
of a parallelogram to a simpler, equivalent problem about finding the area of a 
rectangle. 




area of B ,' 



Figure 1.9: The area of a parallelogram. 



For complicated regions such as the shaded regions in Figure 1.10( a) , we can 
use a similar technique. At first glance, it is difficult to estimate the total area 
of the shaded regions in Figure 1.10( a). However, if we slide all of them into 
a single column as shown in Figure 1.10( b) , then it is easy to determine that 
the shaded area is less than the area of the enclosing rectangle: that is, area of 
rectangle = base x height = 1x2 = 2. 




12 3 4 
Figure 1.10: Estimating the area of an irregular region. 



1.2 Some applications of area 

One reason "areas" are so useful is that they can represent quantities other than 
simple geometric shapes. For example, if the units of the base of a rectangle 
are "hours" and the units of the height are "miles/hour", then the units of 
the "area" of the rectangle are (hours) x (miles/hour) = miles, a measure of 
distance. Similarly, if the base units are centimeters and the height units are 
grams, then the "area" units are gram x centimeters, a measure of work. 

Example 1.2.1. Distance as an "area": In Figure Yl.llj f(t) is the velocity of 
a car in "miles per hour" and t is the time in "hours. " Then the shaded "area" 
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is (base) x (height) = (3 hours) x (20 miles/hour ) = 60 miles. In other words, 
the car traveled 60 miles in the 3 hours from 1 o 'clock until 4 o 'clock. 



fit) = 20 mph 




pm pm pm pm 



clock reading 

Figure 1.11: Calculating distance using the "area" concept. 



Example 1.2.1 shows that we can use the familiar idea of "area" to calculate 
a quantity (in this case, miles) that at first seems totally unrelated to "area". 
Theorem 1.2.1 provides a general statement of the idea illustrated in Exam- 
ple ELD 

Theorem 1.2.1. ("Area" as Distance) Let f(t) be the (positive) forward ve- 
locity of an object at time t. Then the "area" bounded by the graph of f , the 
t-axis, and the vertical lines at times t = a and t = b, will be the distance that 
the object has moved forward between times a and b. 

This "area as distance" fact can make some difficult distance problems much 
easier to handle. 

Example 1.2.2. A car starts from rest (velocity = 0) and steadily speeds up so 
that 20 seconds later its speed is 88 feet per second (60 miles per hour). How 
far did the car travel during those 20 seconds? 

Solution. If "steadily speeds up" means that the velocity increases linearly, then 
we can apply the idea of "area as distance" . Figure 1.12 shows the graph of 
velocity versus time. The "area" of the triangular region in the figure represents 
the distance traveled. Using formula (1.1) for the area of a triangle, we have 

distance = - x base x height 

— - x (20 seconds) x (88 feet/second) 
= 880 feet 

In 20 seconds the car covered a total distance of 880 feet. □ 
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20 

t=me (seconds) 



Figure 1.12: Using "area" to calculate the distance a car travels. 

Practice 1.2.1. A train traveling at 45 miles per hour (66 feet/second) takes 
60 seconds to come to a complete stop. If the train slowed down at a steady rate 
(the velocity decreased linearly), how many feet did the train travel while coming 
to a stop? 

Practice 1.2.2. You and a friend start off at noon and walk in the same di- 
rection along the same path at the rates shown in Figure \1.13[ 

• Who is walking faster at 2 pm? Who is ahead at 2 pm? 

• Who is walking faster at 3 pm? Who is ahead at 3 pm? 

• When will you and your friend be together? (Answer in words.) 




pm pm pm 
clock leading 



Figure 1.13: Illustration for Practice 11.2.21 
1.2.1 Total accumulation as "area" 

In previous examples, the function represented a rate of travel (miles per hour), 
and the area represented the total distance traveled. For functions representing 
other rates, the area still represents the total amount of something. Some 
familiar rates include the production of a factory (bicycles per day) , the flow of 
water in a river (gallons per minute), traffic over a bridge (cars per minute), or 
the spread of a disease (newly sick people per week). 
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Theorem 1.2.2. ("Area" as Total Accumulation) Let f(t) be a positive rate 
(in units per time interval) at time t. Then the "area" between the graph of f, 
the t-axis, and the vertical lines at times t = a and t = b, will be the total units 
which accumulate between times a and b. 

Practice 1.2.3. Figure \1-14\ shows the number of telephone calls made per hour 
on a Tuesday. Approximately how many calls were made between 9 am and 11 
am? 




Figure 1.14: Illustration for Practice 11.2.31 



1.2.2 Problems 

1. (a) Calculate the sum of the rectangular areas in Figure [1.15( a). 

(b) From part (a) , what can we say about the area of the shaded region 
in Figure [1.15( b)? 

2. (a) Calculate the sum of the areas of the shaded regions in Figure [TTTBT c). 

(b) From part (a) , what can we say about the area of the shaded region 
in Figure [1.15( b)? 




Figure 1.15: Estimating areas. 



3. Consider the function f(x) whose graph is shown in Figure d. 16 [ Let A(x) 
represent the area bounded by the graph of f(x), the horizontal axis, and 
vertical lines at t = and t = x. Evaluate A(x) for x = 1, 2, 3, 4, and 5. 

4. Police chase: A speeder traveling 45 miles per hour (in a 25 mph zone) 
passes a stopped police car which immediately takes off after the speeder. 
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~t — f- 

3 "4- 



Figure 1.16: Computing areas. 



The police car speeds up steadily to 60 miles/hour in 20 seconds and then 
travels at a steady 60 miles/hour. Figure 1.171 shows the velocity versus 
time graph for the police car and the speeder. How long and how far 
before the police car catches the speeder who continued traveling at 45 
miles/hour? 




time (seconds} 



Figure 1.17: Computing areas to analyze the police chase problem. 

5. What are the units for the "area" of a rectangle with the given base and 
height units? 



Base units 


Height units 


"Area" units 


miles per second 


seconds 




hours 


dollars per hour 




square feet 


feet 




kilowatts 


hours 




houses 


people per house 




meals 


meals 





1.3 Sigma notation and Riemann sums 

So far, we have calculated the area of a region in terms of simpler regions. 
We cut the region into simple shapes, calculate the area of each simple shape, 
and then add these smaller areas together to get the area of the whole region. 
We will continue using that approach, but it is useful to have a notation for 
adding many values together. This is where the sigma notation (denoted with 
the capital Greek letter E) comes in handy. 
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Consider the function f(x) = 3x + 1 and for convenience we restrict x to 
take only integer values. Now consider the values of f(x) where x = 2,3,4,5 
as shown in Table 11.11 The sigma notation allows us to express the sum S = 
/(2) + /(3) + /(4) + /(5) using a convenient shorthand, namely 

5 

S = 3fc + 1 = /(2) + /(3) + /(4) + /(5) (1.2) 

k=2 



X 


/(*) = 


3x + 1 


2 


/(2) 


= 7 


3 


/(3) 


= 10 


4 


/(4) 


= 13 


5 


/(5) 


= 16 



Table 1.1: Values of f(x) = 3x + 1 for x = 2, 3, 4, 5. 

Figure 1.18 shows the various parts of this notation. The function to the right 
of the sigma E is called the summand. Below and above the sigma are numbers 
called the lower and upper limits of the summation, respectively. Notice that 
instead of using the symbol x, in equation ( |1.2j) and Figure 1.18 we replaced x 
with the variable k. The variable (typically i, j, or k) used in the summation 
is called the counter or index variable. Table 1.2 provides various ways of 
reading the E notation. Using Sage, we can obtain the value of the sum in 
Figure 1.18 as follows: 

sage: maxima (" sum (3*k + 1, k, 2, 5)") 
46 



upper limit 

^~ 3k + 1 
k = 2 

^ § summand 

variable lower 
index ii mit 

Figure 1.18: The summation or sigma notation. 
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Summation 
notation 


A way to read 
the sigma notation 


Sigma 
notation 


l 2 + 2 2 + 3 2 + 4 2 + 5 2 


the sum of k squared 
for k equals 1 to k equals 5 


C o 1 

ELifc 2 


1,1,1,1,1 

3'4~5'6'7 


the sum of 1 over k 
for k equals 3 to k equals 7 




2° + 2 1 + 2 2 + 2 3 + 2 4 + 2 5 


the sum of 2 to the j-th power 
for j equals to j equals 5 




a 2 + a 3 + a 4 + a 5 + a 6 + a 7 


the sum of a sub i 
for i equals 2 to i equals 7 


Ej=2 a * 



Table 1.2: How to read the sigma notation. 



X 


/(*) 


9(x) 


h(x) 


1 


2 


4 


3 


2 


3 


1 


3 


3 


1 


-2 


3 


4 





3 


3 


5 


3 


5 


3 



Table 1.3: Domain values and their images under the functions /, g, and h. 

Practice 1.3.1. Write the summation denoted by each of the following: 

(a) ELik 3 , (b) E; =2 M) j }. (c) El=o(2m+l). 

In practice, the sigma notation is frequently used with the standard function 
notation: 

3 

£ f(k + 2) = f(l + 2) + /(2 + 2) + /(3 + 2) = /(3) + /(4) + /(5) 
fe=l 

and 

4 

J2 f(xi) = f{xi) + f{x 2 ) + f(x 3 ) + f(x 4 ). 

k = l 

Example 1.3.1. Use the values in Table \J73\ to evaluate (i) Efc=2 2/(^);' an d 
(H)Y,U (5 + /(i-2)). 

Solution, (i) Using the values of f(k) for k = 2, 3, 4, 5, we have 

5 

^2/(fc) = 2/(2) + 2/(3) + 2/(4) + 2/(5) 

fe=2 (1.3) 
= 2(3) +2(1) + 2(0) + 2(3) 

= 14 
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On the other hand, note that in expression ( jl.3|) we can factor out the value 2 
to get an equivalent expression 



k=2 k=2 

This expression evaluates to the same value as the original expression. 

(ii) For j = 3, 4, 5 we have 

/(3-2) = /(l), /(4-2) = /(2), /(5-2)=/(3) 
Then we get 

E(5 + /0'-2)) 

= ( 5 + /(3-2)) + (5 + /(4-2)) + (5 + /(5-2)) (L4) 
= (5 + /(l)) + (5 + /(2)) + (5 + /(3)) 
= (5 + 2) + (5 + 3) + (5 + 1) 
= 21 

Now carefully look at expression ( jl.4|) again. When the index j takes on each 
of the three values 3, 4, and 5, the term 5 is added three times in evaluating the 
whole sum. Using this observation, we can write ( |1.4j ) in the equivalent form 

5 5 

£( 5 + /0*-2)) =3x5 + £/Cj'-2) 

J=3 j=3 

This observation will be formalized later on. □ 

Practice 1.3.2. Use the values of f, g and h in Table \1.3\ to evaluate the 
following: 

(a) ]T 5 (fc), (6) (c) £[/(»- !) + <#)]. 

k=2 j = l i=3 

Since the sigma notation is simply a notation for addition, it has all of the 
familiar properties of addition. Theorem 1.3.1 states a number of these familiar 
addition properties for the E notation. 

Theorem 1.3.1. (Summation Properties) Let C be a fixed constant and let 
m, n G Z be positive. Then we have the following summation properties for the 
S notation. 

1. Sum of constants: X)fe=i C ~ C + C + C + ■ ■ ■ + C = nC . 

n terms 
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2. Addition: Yl=i( a k + h) = Ylk=i a k + YJl=i b k- 

3. Subtraction: Ylk=ii a k - h) = YZ=i a k - J2l=i b k- 
4- Constant multiple: X)fc=i ^ a k = C52fc=i a k- 

5. Preserves positivity: If bk > Ofe for all k, then J2fc=i ^fc — Sfe=i a fc- ^ n 
particular, if a^ > /or aH /e i/ien £\ =1 a fc > 0. 

6. Additivity of ranges: If 1 < m < n f/ien X)fc=i a fe + Sfe=m+i a *-' = 

En 

1.3.1 Sums of areas of rectangles 

Later, we will approximate the areas under curves by building rectangles as 
high as the curve, calculating the area of each rectangle, and then adding the 
rectangular areas together. 

Example 1.3.2. Evaluate the sum of the rectangular areas in Figure [i.l9\ and 
write the sum using the sigma notation. 
Solution: We have 

sum of the rectangular areas = sum of (base) x (height) for each rectangle 

= (l)(l/3) + (l)(l/4) + (l)(l/5) = 47/60. 

Using the sigma notation, 

(l)(l/3) + (l)(l/4) + (l)(l/5) = £i 

fe=i 




12 3 4 5 



Figure 1.19: Area and summation notation. 

Practice 1.3.3. Evaluate the sum of the rectangular areas in Figure [i . 20\ and 
write the sum using the sigma notation. 
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12 3 4 5 



Figure 1.20: Area and summation notation. 



Example 1.3.3. Write the sum of the areas of the rectangles in Figure [1.21 
using the sigma notation. 
Solution: The area of each rectangle is (base)x (height). 



rectangle 


base 


height 


area 


1 


x\ - x 




(xi 


- x )f(x 1 ) 


2 


x 2 - Xi 


f(X2) 


(X2 


- xi)f(x 2 ) 


3 


x 3 - x 2 


/(*a) 


(x 3 


- x 2 )f(x 3 ) 



The area of the k-th rectangle is (xk — Xk-i).f{xk), and the total area of the 
rectangles is the sum J2k=i( Xk ~ x k-i)f{xk)- 




Figure 1.21: Area and summation notation. 
1.3.2 Area under a curve: Riemann sums 

Suppose we want to calculate the area between the graph of a positive function 
/ and the interval [a, b] on the x-axis (Fig. 7). The Riemann Sum method is to 
build several rectangles with y = f(x) bases on the interval [a, b] and sides that 
reach up to the graph of / (Fig. 8). Then the areas of the rectangles can be 
calculated and added together to get a number called a Riemann sum of / 
on [a, b]. The area of the region formed by the rectangles is an approximation 
of the area we want. 
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Example 1.3.4. Approximate the area in Figure [7. 22\ ( a) between the graph of 
f and the interval [2, 5] on the x-axis by summing the areas of the rectangles in 
Figure U7M\ (b). 

Solution: The total area of rectangles is (2)(3) + (1)(5) = 11 square units. 




In order to effectively describe this process, some new vocabulary is helpful: 
a "partition" of an interval and the mesh of the partition. A partition P 
of a closed interval [a, b] into n subintervals is a set of n + 1 points {xq = 
a, X\, X2, £3, x n -x, x n = b} in increasing order, a = x < xi < x 2 < X3 < 
... < x n — b. (A partition is a collection of points on the axis and it does not 
depend on the function in any way.) 

The points of the partition P divide the interval into n subintervals (Figure 
11.23 j ): [x ,xi], [xi,x 2 ], [0:2, X3] j ■ ■ ■ , an d [x n -i,x n ] with lengths Axi = x% — xq, 
Ax 2 = X2 — X\, . . . , Ax n = x n — x n _±. The points Xk of the partition P are 
the locations of the vertical lines for the sides of the rectangles, and the bases 
of the rectangles have lengths Axk for k = 1,2,3, ...,rt. 



first second 

subinterva] subinterual sub interval 



i 



—\ 1 1 1 1 1 1 



Figure 1.23: Partition of the interval [a, b]. 



The mesh or norm of the partition is the length of the longest of the subin- 
tervals [xk-i, Xk], or, equivalently, the maximum of the Acc/c for k = 1, 2, 3, n. 
For example, the set P = {2, 3, 4.6, 5.1, 6} is a partition of the interval [2, 6] and 
divides it into 4 subintervals with lengths Axi = 1, Ax 2 = 1.6, A.T3 = 0.5 and 
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A.T4 = 0.9. The mesh of this partition is 1.6, the maximum of the lengths of 
the subintervals. (If the mesh of a partition is "small," then the length of each 
one of the subintervals is the same or smaller.) 

A function, a partition, and a point in each subinterval determine a Rie- 
mann sum. Suppose / is a positive function on the interval [a, 6], P = {xq = 
a, x\, X2, £3, x n —i, x n = 6} is a partition of [a, b], and c k is an revalue in the k- 
th subinterval [xk-i,Xk] : Xk-i < c k < Xk- Then the area of the fc-th rectangle 
is /(c fe ) ■ (x k - x k -{) = /(cfe)Axfc. (Figure 




Figure 1.24: Part of a Ricmann sum. 

Definition 1.3.1. A summation of the form Y^l— 1 f(ck)Axk is called a Rie- 
mann sum of f for the partition P. 

This Ricmann sum is the total of the areas of the rectangular regions and is 
an approximation of the area between the graph of / and the x-axis. 

Example 1.3.5. Find the Riemann sum for f(x) = 1/x and the partition 
{1, 4, 5} using the values C\ = 2 and c 2 = 5. 

Solution: The two subintervals are [1,4] and [4,5] so Axi = 3 and Ax 2 = 1. 
Then the Riemann sum for this partition is 

n 1 

J2f(c k )Ax k = f(c 1 )Ax 1 + f(c 2 )Ax 2 = /(2)(3) + /(5)(l) = ^Hgt 1 ) = L7 ' 

k=l 

Practice 1.3.4. Calculate the Riemann sum for f{x) = l/x on the partition 
{1,4,5} using the values c\ = 3, c 2 = 4. 

Practice 1.3.5. What is the smallest value a Riemann sum for f{x) = l/x 
and the partition {1,4,5} can have? (You will need to select values for c\ and 
c 2 .) What is the largest value a Riemann sum can have for this function and 
partition? 

Here is a Sage example. 

Example 1.3.6. Using Sage, we construct the Riemann sum of the function 
y = x 2 using a partition of 6 equally spaced points, where the c k 's are taken to 
be the midpoints. 
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sage: fl(x) = x~2 

sage: f = Piecewise( [ [(-1 , 1) ,f 1] ] ) 

sage: g = f . riemann_sum(6 ,mode="midpoint" ) 

sage: P = f .plot (rgbcolor=(0 .7,0 . 1 ,0. 5) , plot_points=40) 

sage: Q = g. plot (rgbcolor= (0.7, 0.6,0. 6) , plot_points=40) 

sage: L = add( [line( [ [pf [0] [0] , 0] , [pf [0] [0] ,pf [1] (x=pf [0] [0] )] ] , \ 

rgbcolor=(0.7,0.6,0.6)) for pf in g.listO]) 

sage: show(P+Q+L) 



Here is the plot: 




-0,25 - 
-0.5 - 
-0.75 - 



Figure 1.25: Plot using Sage of a Riemann sum for y = x 2 . 

At the end of this section is a Python^ program listing for calculating Riemann 
sums of f{x) = l/i on the interval [1,5] using 100 subintervals. It can be 
modified easily to work for different functions, different endpoints, and different 
numbers of subintervals. Table 1.261 shows the results of running the program 
with different numbers of subintervals and different ways of selecting the points 
Ci in each subinterval. When the mesh of the partition is small (and the number 
of subintervals large) , all of the ways of selecting the Ci lead to approximately 
the same number for the Riemann sums. 



2 Python is a general purpose cross-platform, free and open source computer program- 
ming language. It is widely used in industry and academia, and is available for download at 
http : //www . py thon ■ org) Alternatively , you can download the mathematical software system 
Sage from |httpT / /www . sagemath . org. It comes with Python pre-installed. 
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Here is a Python program to calculate Riemann sums of f(x) = 1/x on [1,5] 
using 100 equal length subintervals, based on the "left-hand" endpoints. 



f = lambda x: 1/x 

a = 1.0 

b = 5.0 

n = 100 

Dx = (b-a)/n 

rsum = sum( [f (a+i*Dx) *Dx for i in 
print rsum 



# define the function 

# left endpoint of integral 

# right endpoint of integral 

# number of subintervals 

# width of each subinterval 
(n)]) # compute the Riemann sum 

# print the Riemann sum 



Other Riemann sums can be calculated by replacing the "rsum" line with one 
of: 

rsum = sum( [f (a+(i+0 . 5) *Dx) *Dx for i in range (n)]) "midpoint" 
rsum = sum( [f (a+(i+l)*Dx)*Dx for i in range(n)]) "right-hand" 

Written as Python "functions" , these three are written as below@ 



def rsum_lh(n) : 

f = lambda x: 1/x 

a = 1.0 

b = 5.0 

Dx = (b-a)/n 

return sum( [f (a+i*Dx) *Dx for i in range (n)]) 

def rsum_mid(n) : 

f = lambda x: 1/x 

a = 1.0 

b = 5.0 

Dx = (b-a)/n 

return sum( [f (a+(i+0 . 5) *Dx) *Dx for i in range (n)]) 

def rsum_rh(n) : 

f = lambda x: 1/x 

a = 1.0 

b = 5.0 

Dx = (b-a)/n 

return sum( [f (a+(i+l)*Dx)*Dx for i in range (n)]) 



3 If you have an electronic copy of this file, and "copy-and-pastc" 
these into Python, keep in mind indenting is very important in Python 
syntax. http: //www. python. org/doc/essays/styleguide .html 
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left-hand 


midpoint 


right-hand 


n 


Axi 


Riemann sum 


Riemann sum 


Riemann sum 


5 


0.8 


1.9779070602600015 


1.5861709609993364 


1.3379070602600014 


10 


0.4 


1.7820390106296689 


1.6032106782106783 


1.462039010629669 


20 


0.2 


1.6926248444201737 


1.6078493243021688 


1.5326248444201738 


100 


0.04 


1.6255658911511259 


1.6093739310551827 


1.5935658911511259 


1000 


0.004 


1.6110391924319691 


1.6094372724359669 


1.607839192431969 



Figure 1.26: Table for Python example. 



The command 

sizes = [5, 10, 20, 100, 1000] 

table = [[n, (b-a)/n, rsum_lh(n) ,rsum_mid(n) , rsum_rh(n)] for n in sizes] 
yields the following data: 

In fact, the exact value is log(5) = 1.609437..., so these last few lines yielded 
pretty good approximations. 

Practice 1.3.6. Replace 1/x by x 2 and [a,b] = [1,5] by [a,b] = [—1,1] in the 
Python code above and find the Riemann sum for the new function andn = 100. 
Use the midpoint approximation. (You may use Sage or Python, whichever you 
prefer.) 

Example 1.3.7. Find the Riemann sum for the function f(x) = sin(x) on 
the interval [0,7r] using the partition {0, 7r/4, 7t/2, tt} with c\ = 7r/4, C2 = 7r/2, 
c 3 = 3tt/4. 

Solution: The 3 subintervals are [0, 7t/4] ; [7t/4, 7t/2], and [tt/2,tt] so Ax\ = 
7r/4, Ax2 = 7r/4 and Axg, = n/2. The Riemann sum for this partition is 

ELi f(ck)Ax k = sin(7r/4)(7r/4) + sin(^/2)(^/4) + sin(37r/4)(7r/2) 

1 TT _|_ 1 7T _|_ 1 7T 

= 2.45148... . 

Practice 1.3.7. Find the Riemann sum for the function and partition in the 
previous example, but use C\ = 0, ci = tt/2, C3 = 7r/2. 

1.3.3 Two special Riemann sums: lower and upper sums 

Two particular Riemann sums are of special interest because they represent the 
extreme possibilities for Riemann sums for a given partition. 

Definition 1.3.2. Suppose f is a positive function on [a, b], and P is a partition 
of [a,b]. Let be the xvalue in the k-th subinterval so that f(mk) is the 
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minimum value of f in that interval, and let M k be the xvalue in the k-th 
subinterval so that f(M k ) is the maximum value of f in that interval. 

lower sum: LS = Y^ik=i /( m fc)Axfe. 

upper sum: US = £)fe=i f(M k )Ax k . 

Geometrically, the lower sum comes from building rectangles under the graph 
of / (Figure ^. 27( a)), and the lower sum (every lower sum) is less than or equal 
to the exact area A: LS < A for every partition P. The upper sum comes from 
building rectangles over the graph of / (Figure 1.27( b)). and the upper sum 
(every upper sum) is greater than or equal to the exact area A: US > A for 
every partition P. The lower and upper sums provide bounds on the size of the 
exact area: LS < A < US. 




Unfortunately, finding minimums and maximums can be a timeconsuming 
business, and it is usually not practical to determine lower and upper sums for 
"arbitrary" functions. If / is monotonic, however, then it is easy to find the 
values for m k and M/. , and sometimes we can explicitly calculate the limits of 
the lower and upper sums. 

For a monotonic bounded function we can guarantee that a Riemann sum is 
within a certain distance of the exact value of the area it is approximating. 

Theorem 1.3.2. If f is a positive, monotonically increasing, bounded function 
on [a, b] , then for any partition P and any Riemann sum for P, 

distance between the Riemann sum and the exact area < distance between the 
upper sum (US) and the lower sum (LS) < (/(£>) — /(a)) ■ (mesh of P). 

Proof: The Riemann sum and the exact area are both between the upper 
and lower sums so the distance between the Riemann sum and the exact area is 
less than or equal to the distance between the upper and lower sums. Since / is 
monotonically increasing, the areas representing the difference of the upper and 
lower sums can be slid into a rectangle whose height equals f(b) — f(a) and whose 
base equals the mesh of P. Then the total difference of the upper and lower 
sums is less than or equal to the area of the rectangle, (/(&) — f(a))- (mesh of P). 
□ 
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1.3.4 Problems 

For problems the next four problems, sketch the function and find the smallest 
possible value and the largest possible value for a Riemann sum of the given 
function and partition. 



1 

1 . 


J\ X ) ~ 


1 i t2 
= 1 + X 




(a\ P 

(a) f 


— i 1 ! A 4 j £>> 




(h) P 


— 11 9 3 4 51 

— -[1, z, o,ft, J_f 




(c\ P 


— /I 1 c. o o 4 cl 

— \L, l.O, Z, O, 'i, Oj- 


2. 


/(*) = 


= 7 -2a; 




(a)P 


= {0,2,3} 




(b)P 


= {0,1,2,3} 




(c)P 


= {0, .5, 1,1.5, 2, 3} 


3. 


m = 


- sin(a;) 




(a)P 


= {0,7r/2,7r} 




(b)P 


= {0,^/4,77/2,^} 




(c)P 


= {0,7r/4, 3tt/4, tt}. 


4. 


m = 


= a; 2 - 2x + 3 




(a)P 


= {0,2,3} 




(b)P 


= {0,1,2,3} 




(c)P 


= {0, .5, 1,2, 2.5,3} 



5. Suppose we divide the interval [1, 4] into 100 equally wide subintervals and 
calculate a Riemann sum for f(x) = 1 + x 2 by randomly selecting a point 
Ci in each subinterval. 

(a) We can be certain that the value of the Riemann sum is within what 
distance of the exact value of the area between the graph of / and the 
interval [1,4] ? 

(b) What if we take 200 equally long subintervals? 

6. If / is monotonic decreasing on [a, b] and we divide the interval [a, b] into 
n equally wide subintervals, then we can be certain that the Riemann sum 
is within what distance of the exact value of the area between f and the 
interval [a, 6]? 

7. Suppose LS = 7.362 and US = 7.402 for a positive function / and a 
partition P of the interval [1,5]. 

(a) We can be certain that every Riemann sum for the partition P is 
within what distance of the exact value of the area under the graph of / 
over the interval [1,5]? 

(b) What if LS = 7.372 and US = 7.390? 
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1.3.5 The trapezoidal rule 

This section includes several techniques for getting approximate numerical val- 
ues for definite integrals without using antiderivatives. Mathematically, exact 
answers are preferable and satisfying, but for most applications, a numerical 
answer with several digits of accuracy is just as useful. For instance, suppose 
you are a automotive or aircraft designer. You may wish to know how much 
metal is required to build your design, created using a computer-aided design 
graphics program. Due to the fluctuations in the price for metal, you only need 
the approximate cost based on a piecewise-linear approximation to your design. 
Numerical techniques such as those discussed in this section can be used for 
that. 

The methods in this section approximate the definite integral of a function / 
by building "easy" functions close to / and then exactly evaluating the definite 
integrals of the "easy" functions. If the "easy" functions are close enough to 
/, then the sum of the definite integrals of the "easy" functions will be close 
to the definite integral of /. The Left, Right and Midpoint approximations fit 
horizontal lines to / , the "easy" functions arc pieccwise constant functions, and 
the approximating regions are rectangles. The Trapezoidal Rule fits slanted lines 
to / , the "easy" functions are piecewise linear, and the approximating regions 
are trapezoids. Finally, Simpson's Rule fits parabolas to /, and the "easy" 
functions are piecewise quadratic polynomials. 

All of the methods divide the interval [a, b] into n equally long subintervals. 
Each subinterval has length h = Axt = and the points of the partition 

are xq = a, x\ = a + h, x% = a + 2h, . . . , Xi = a + ih, . . . , x„ = a + nh = 

a + n (b=*) = b. 

If the graph of f is curved, then slanted lines typically come closer to the graph 
of f than horizontal ones do, and the slanted lines lead to trapezoidal regions. 

The area of a trapezoid is (base) x (average height), so the area of the trapezoid 
with coordinates (cco,0), (£o,2/o), (^i,0), (2:1,2/1), is (see Figure fl.28j ). 



Example 1.3.8. Here is how to create and plot a piecewise linear function 
describing the trapezoidal approximation to the area under y — x 3 — 3x 2 + 2x. 



sage: x = var("x") 

sage: fl = lambda x: x~3-3*x~2+2*x 

sage: f = Piecewise ( [ [(0 , 2) ,f 1] ] ) 

sage: tf = f . trapezoid (4) 

sage: P3 = list_plot( [(1/2,0) , (0.5, f (0.5))] ,plotjoined=True,linestyle= 

sage: P4 = list_plot ( [(3/2, 0) , (1 . 5 , f (1 . 5) )] ,plot joined=True , linestyle= 

sage: show(Pl+P2+P3+P4) 

sage : f . trapezoid_integral_approximation(4) 
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0=1 >3/i ) 



foi .Sfi ) 



Figure 1.28: Trapezoid. 

sage: integrate (x~3-3*x~2+2*x, x, 0, 2) 




Here is the plot: 



0.75 - 
0.5^ 
0.25 

-1 -0.5 

-0.25 

-0.5 - 
-0.75 - 
-1 - 

Figure 1.29: Plot using Sage of a trapezoidal approximation to the integral 
J* x 3 -3x 2 + 2xdx. 

We got lucky here since both the integral and its trapezoidal approximation ac- 
tually have the same value. 

Theorem 1.3.3. ("Trapezoidal Approximation Rule") If f is integrable on 
[a,b], and [a,b] is partitioned into n subintervals of length h = then the 

Trapezoidal approximation of f(x) dx is 

Tn = \[f{xo) + 2f(xt) + 2/(2- 2 ) + ... + 2/(x n _r) + f(x n )}. 
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X 


/(*) 


1.0 


4.2 


1.5 


3.4 


2.0 


2.8 


2.5 


3.6 


3.0 


3.2 



Figure 1.30: Tabic for [LAS 

Proof: The area of the trapezoidal with coordinates {xi, 0), (xi, yi), (xj+i, 0), 
(x l+1 ,y i+ i), where yi = f(xi), is (x i+1 - x t ) Vi+ ^ +1 = | • ( yi + y i+ x). Therefore, 

the sum of the trapezoidal areas approximating f{x) dx is 

El \Vi -T y% + l h- ST^ , \ / n n \ 

(Xi+i-Xi) = - • 2_^{Vi +Vi+i) = 2 ' (^0 + 2yi + .~ + 2y n - 1 +y n ), 

i— i— 

as desired. □ 

Example 1.3.9. Calculate T 4 , the Trapezoidal approximation of fix) dx, for 
the function values tabulated in Figure \1.30[ 
Solution: The step size is h = (b — a)/n = (3 — l)/4 = 1/2. Then 

T 4 = |[/(.t ) + 2/(zi) + 2f(x 2 ) + 2f(x 3 ) + f(x 4 )} 
= | [4.2 + 2(3.4) + 2(2.8) + 2(3.6) + (3.2)] 
= (0.25)(27) = 6.75. 

Example 1.3.10. Let's see how well the trapezoidal rule approximates an in- 
tegral whose exact value we know, f^x 2 dx = Calculate T 4 , the Trapezoidal 
approximation of x 2 dx. 

Solution: As in the example above, h = 0.5 and xo = 1, X\ = 1.5, x 2 = 2, 
x 3 = 2.5, and x 4 = 3. Then 

T 4 = ^[f(x Q ) + 2f(x 1 ) + 2f(x 2 ) + 2f(x 3 ) + f(x 4 ]} 
= 0.5[/(l) + 2/(1.5) + 2/(2) + 2/(2.5) + /(3)] 
= (0.25) [1 + 2(2.25) + 2(4) + 2(6.25) + 9] = 8.75. 

Using Sage, one can show that T w = 217/25 = 8.68, T WQ = 21667/2500 = 
8.6668, and T 10 oo = 2166667/250000 = 8.666668. These trapezoidal approxima- 
tions are indeed approaching the value 8.666... of the integral. 

Example 1.3.11. Here is a Python program illustrating the Trapezoidal rule. 

f = lambda x: sin(x) 

def trapezoidal_rule (f cn,a,b,n) : 
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Does computation of the Trapezoidal rule approx of int_a~b fcn(x) dx 
using n steps. 



Deltax = (b-a)*1.0/n 
coeffs = [2]*(n-D 
coeffs = [l]+coeff s+[l] 

valsf = [f (a+Deltax*i) for i in range (n+1)] 

return (Deltax/2) *sum( [coeffs [i] *valsf [i] for i in range(n+l)]) 



Now we paste this into Sage (you may instead paste into Python if you wish) 
and see how it works. 



sage: f = lambda x: sin(x) 

sage: def trapezoidal_rule (f cn,a,b,n) : 
ii n it 

Does computation of the Trapezoidal rule approx of int_a~b fcn(x) dx 
using n steps. 



Deltax = (b-a)*1.0/n 
coeffs = [2]*(n-D 
coeffs = [l]+coeff s+[l] 

valsf = [f (a+Deltax*i) for i in range (n+1)] 

return (Deltax/2) *sum( [coeffs [i] *valsf [i] for i in range (n+1)]) 



sage: integral (f (x) ,x, , 1) 
1 - cos(l) 

sage: RR (integral (f (x) ,x,0, 1) ) 
0.459697694131860 

sage : trapezoidal_rule (f (x) ,0,1,4) 
0.457300937571502 



You see sin(x) dx = 1 — cos(l) = 0.459..., whereas the trapezoidal rule gives 
the approximation T4 = 0.457.... 

1.3.6 Simpson's rule and Sage 

If the graph of / is curved, even the slanted lines may not fit the graph of / as 
closely as we would like, and a large number of subintervals may still be needed 
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with the Trapezoidal rule to get a good approximation of the definite integral. 

1 

(^1,^ 1 ) 




Figure 1.31: Piece of a parabola. 

Theorem 1.3.4. Three points (xo,yo), (xi,yi), (2:2,2/2) are needed to determine 
the equation of a parabola, and the area under a parabolic region with evenly 
spaced Xi values (Figure \1.31\) is 

Ax 

"3- (2/o +%i +V2)- 

Theorem 1.3.5. ("Simpson's Rule") If f is integrable on [a,b], and [a, b] is 
partitioned into an even number n of subintervals of length h = Ax = . then 

the Parabolic approximation of J f(x) dx is 

At 

S n = -y\f(x ) + 4/On) + 2/(0:2) + 4/(2:3) + 2/(14) + ... + 4/(z„_ 1 ) + f(x n )}. 

Example 1.3.12. Calculate S4, the Simpson's rule approximation of f{x) dx, 
for the function values tabulated in Figure \ 1.301 
Solution: The step size is h = (b — a)/n — (3 — l)/4 = 1/2. Then 

Si ='j [f(x ) + 4/(^0 + 2/(2:2) + 4/(a;3) + f(x A )] 
= I [4.2 + 4(3.4) + 2(2.8) + 4(3.6) + (3.2)] 
= Ir = 6.83.... 

Example 1.3.13. Here is a Python program illustrating Simpson's rule. 

f = lambda x: sin(x) 

def simpsons_rule(f cn,a,b,n) : 
11 11 11 
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Does computation of the Simpson's rule approx of int_a~b fcn(x) dx 
using n steps. Here n must be an even integer. 

M II II 

Deltax = (b-a)*1.0/n 

n2 = int(n/2) 

coeffs = [4,2] *n2 

coeffs = [l]+coeffs[:n-l] + [l] 

valsf = [f (a+Deltax*i) for i in range (n+1)] 

return (Deltax/3) *sum( [coeffs [i] *valsf [i] for i in range(n+l)]) 



Now we paste this into Sage and see how it works: 



sage: f = lambda x: sin(x) 

sage: def simpsons_rule (f cn,a,b,n) : 
n n it 

....: Does computation of the Simpson's rule approx of int_a~b fcn(x) 
. . . . : using n steps. Here n must be an even integer. 

n it ii 

: Deltax = (b-a)*1.0/n 

: n2 = int(n/2) 

: coeffs = [4,2] *n2 

....: coeffs = [l]+coeffs[:n-l]+[l] 

....: valsf = [f (a+Deltax*i) for i in range(n+l)] 

....: return (Deltax/3) *sum( [coeffs [i] *valsf [i] for i in range (n+1)]) 

sage: integral (f (x) ,x, , 1) 
1 - cos(l) 

sage: RR (integral (f (x) ,x,0, 1) ) 

0.459697694131860 

sage: simpsons_rule (f (x) ,0, 1 ,4) 

(sin(l) + 4*sin(3/4) + 2*sin(l/2) + 4*sin(l/4) ) /12 
sage: RR(simpsons_rule (f (x) ,0,1,4)) 
0.459707744927311 

sage : RR(simpsons_rule (f (x) ,0 , 1 , 10) ) 
0.459697949823821 



To paste this into Python, you mus first import the sin function. 



»> from math import sin 
»> f = lambda x: sin(x) 
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»> def simpsons_rule(f cn,a,b,n) : 

M II II 

... Does computation of the Simpson's rule approx of int_a~b fcn(x) 
... using n steps. Here n must be an even integer. 

n n it 

Deltax = (b-a)*1.0/n 
n2 = int(n/2) 
coef f s = [4,2] *n2 
coeffs = [l]+coeffs[:n-l] + [l] 
... valsf = [f (a+Deltax*i) for i in range(n+l)] 

... return (Deltax/3) *sum( [coef f s [i] *valsf [i] for i in range(n+l)]) 

»> simpsons_rule (f ,0, 1 ,4) 
0.45970774492731092 



Using this and the trapezoidal approximation function, we can compare which is 
best in this example. 

sage: simpsons_rule (f (x) , , 1 ,4) 
0.459707744927311 

sage : trapezoidal_rule (f (x) ,0,1,4) 

0.457300937571502 

sage: simpsons_rule (f (x) ,0, 1 , 10) 

0.459697949823821 

sage : trapezoidal_rule (f (x) ,0,1,10) 

0.459314548857976 

sage: integral (f (x) ,x, , 1) 

1 - cos(l) 

sage: integraKf (x) ,x, , 1) *1 . 
1. 00000000000000* (1 - cos(D) 
sage: RR(integral (f (x) ,x, , 1) ) 
0.459697694131860 



We see Simpson's rule wins every time. 

1.3.7 Trapezoidal vs. Simpson: Which Method Is Best? 

The hardest and slowest part of these approximations, whether by hand or by 
computer, is the evaluation of the function at the Xi values. For n subintervals, 
all of the methods require about the same number of function evaluations. The 
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rest of this section discusses "error bounds" of the approximations so we can 
know how close our approximation is to the exact value of the integral even if 
we don't know the exact value. 

Theorem 1.3.6. (Error Bound for Trapezoidal Approximation) If the second 
derivative of f is continuous on [a,b] and M 2 > maxj-guw |/"(a;)|, then the 
"error of the T n approximation" is 

I Jj(x)dx-T n \<^M 2 . 

The "error bound" formula fe^ M 2 for the Trapezoidal approximation is a 
"guarantee" : the actual error is guaranteed to be no larger than the error bound. 
In fact, the actual error is usually much smaller than the error bound. The word 
"error" does not indicate a mistake, it means the deviation or distance from the 
exact answer. 

Example 1.3.14. How large must n be to be certain that T n is within 0.001 of 
Jq 1 sin(x) dx ? 

Solution: We want to pick n so that { j^M 2 < 1/1000. We may take M 2 = 1 , 
so < £^M 2 = < 1/1000, or n 2 > 1000/12. Taking n = 10 will work. 

Theorem 1.3.7. (Error Bound for Simpson's Rule Approximation) If the fourth 
derivative of f is continuous on [a,b] and M4 > max j£ [ a ^] |/( 4 '(x)| ; then the 
"error of the S n approximation" is 

|/ a b /(,)d,-5„|<^M, 

Example 1.3.15. How large must n be to be certain that S n is within 0.001 of 
J sin(cc) dx ? 

Solution: We want to pick n so that j^Lm 2 < 1/1000. We may take M 2 = 1 , 
so it& M 4 = Islbr ^ 1/1000, or n 4 > 1000/180. Taking n = 2 will work. 

1.4 The definite integral 

Each particular Riemann sum depends on several things: the function /, the 
interval [a, b], the partition P of the interval, and the values chosen for Cf~ in 
each subintcrval. Fortunately, for most of the functions needed for applications, 
as the approximating rectangles get thinner (as the mesh of P approaches 
and the number of subintervals gets bigger) the values of the Riemann sums 
approach the same value independently of the particular partition P and the 
points Cfe. For these functions, the limit (as the mesh approaches 0) of the 
Riemann sums is the same number no matter how the c^'s are chosen. This 
limit of the Riemann sums is the next big topic in calculus, the definite integral. 
Integrals arise throughout the rest of this book and in applications in almost 
every field that uses mathematics. 
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Definition 1.4.1. If \vc&xnesh(P)^>Q Sfe=i f( c k)Axk equals a finite number I then 
f is said to be (Riemann) integrable on the interval [a, b] . 
The number I is called the definite integral of f over [a, b] and is written 

S a /0) dx - 

The symbol J f(x)dx is read "the integral from a to b of / of x dee x"or 
"the integral from a to 6 of f(x) with respect to x." The lower limit is a, 
upper limit is b, the integrand is f(x), and x is sometimes called the dummy 
variable. Note that / f(u) du numerically means exactly the same thing, but 

with a different dummy variable. The value of a definite integral f f{%) dx 
depends only on the function / being integrated and on the endpoints a and 
b. The following integrals each represent the integral of the function / on the 
interval [a, b], and they are all equal: 

f f(x)dx= I f{t)dt= I f{u)du= f f(z)dz. 

J a J a J a J a 

Also, note that when the upper limit and the lower limit are the same then the 
integral is always 0: 




f{x)dx = 0. 



There are many other properties, as we will see later. 

Example 1.4.1. (Relation between velocity and area) 

Suppose you're reading a car magazine and there is an article about a new 
sports car that has this table in it: 



Time (seconds) 





1 


2 


3 


4 


5 


6 


Speed (mph) 





5 


15 


25 


40 


50 


60 



They claim the car drove l/8th of a mile after 6 seconds, but this just 'feels v 
wrong... Hmmm... Let's estimate the distance driven using the formula 

distance = rate x time. 

We overestimate by assuming the velocity is a constant equal to the max on each 
interval: 

195 

estimate = 5 • 1 + 15 • 1 + 25 ■ 1 + 40 • 1 + 50 • 1 + 60 ■ 1 = miles = 0.054... 

3600 

(Note: there are 3600 seconds in an hour.) But 1/8 ~ 0.125, so the article 
is inconsistent. (Doesn't this sort of thing just bug you? By learning calculus 
you'll be able to double-check things like this much more easily.) 

Insight! The formula for the estimate of distance traveled above looks exactly 
like an approximation for the area under the graph of the speed of the car! In 
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fact, if an object has velocity v(t) at time t, then the net change in position from 
time a to b is 

I v(t)dt. 

J a 

If / is a velocity, then the integrals on the intervals where / is positive measure 
the distances moved forward; the integrals on the intervals where / is negative 
measure the distances moved backward; and the integral over the whole time 
interval is the total (net) change in position, the distance moved forward minus 
the distance moved backward. 

Practice 1.4.1. A bug starts at the location x — 12 on the x-axis at 1 pm and 
walks along the axis in the positive direction with the velocity shown in Figure 
1 . 32[ How far does the bug travel between 1 pm and 3 pm, and where is the bug 
at 3 pm? 




Figure 1.32: Velocity of a bug on the x-axis. 



Practice 1.4.2. A car is driven with the velocity west shown in Figure \1 . 33\ 

(a) Between noon and 6 pm how far does the car travel? 

(b) At 6 pm, where is the car relative to its starting point (its position at 
noon)? 




Figure 1.33: Velocity of a car on the x-axis. 



Units For the Definite Integral We have already seen that the "area" under 
a graph can represent quantities whose units are not the usual geometric units 
of square meters or square feet. In general, the units for the definite integral 
f f(x)dx are (units for f(x))x (units for x). A quick check of the units can 
help avoid errors in setting up an applied problem. 

For example, if a; is a measure of time in seconds and f(x) is a velocity with 
units feet/second, then Ax has the units seconds and f(x)Ax has the units 
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(feet/second) (seconds) = feet, a measure of distance. Since each Riemann sum 
S fi x k)Axk is a sum of feet and the definite integral is the limit of the Riemann 
sums, the definite integral f(x)dx, has the same units, feet. 

If f(x) is a force in grams, and x is a distance in centimeters, then f(x)dx 
is a number with the units "gram-centimeters," a measure of work. 

1.4.1 The Fundamental Theorem of Calculus 

Example 1.4.2. For the function f(t) = 2, define A(x) to be the area of the 
region bounded by the graph of f, the t-axis, and vertical lines at t = 1 and 
t = x. 

(a) Evaluate A(l), A{2), A(3), A(4). 

(b) Find an algebraic formula for A(x), for x > 1. 

(c) Calculate ^A(x). 

(d) Describe A(x) as a definite integral. 

Solution : (a) A(l) = 0, A(2) = 2, A(3) = 4, A{A) = 6. (b) A(x) = area of a 
rectangle = (base)x (height) = (rc-l)-(2) = 2x-2. (c) £A(x) = £(2x-2) = 2. 
(d) A(x) = )*2dt. 

Practice 1.4.3. Answer the questions in the previous Example for f(t) = 3. 

A curious "coincidence" appeared in this Example and Practice problem: the 
derivative of the function defined by the integral was the same as the integrand, 
the function "inside" the integral. Stated another way, the function defined 
by the integral was an "antiderivative" of the function "inside" the integral. 
We will see that this is no coincidence: it is an important property called The 
Fundamental Theorem of Calculus. 

Let / be a continuous function on the interval [a,b]. 

Theorem 1.4.1. ("Fundamental Theorem of Calculus") If F{x) is any differ- 
entiable function on [a,b] such that F'(x) = fix), then 




f(x)dx = F(b) - F(a). 



The above theorem is incredibly useful in mathematics, physics, biology, etc. 
One reason this is amazing, is because it says that the area under the entire 
curve is completely determined by the values of a ("magic") auxiliary function 
at only 2 points. It's hard to believe. It reduces computing ( |1.4.1|) to finding 
a single function F, which one can often do algebraically, in practice. Whether 
or not one should use this theorem to evaluate an integral depends a lot on the 
application at hand, of course. One can also use a partial limit via a computer 
for certain applications (numerical integration). 
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Example 1.4.3. I've always wondered exactly what the area is under a "hump" 
of the graph o/sin. Let's figure it out, using F(x) = — cos(x). 




sin(x)dx = - cos(tt) - (- cos(O)) = -(-!)-(-!) = 2. 



In Sage, you can do this both "algebraically" and "numerically" as follows. 

sage: f = lambda x: sin(x) 
sage: integral(f (x) ,x,0,pi) 
2 

sage : numerical_integral (f (x) ,0,pi) 

(1 . 9999999999999998 , 2 . 2204460492503128e-14) 



For the "algebraic" computation, Sage knows how to integrate sin(ar) exactly, 
so can compute Jq sin(x)dx = 2 using its integral command^. On the last 
line of output, the first entry is the approximation, and the second is the error 
bound. For the "numerical" computation, Sage obtain^ the approximation 
J sin(x)dx ?s 1.99999... by taking enough terms in a Riemann sum to achieve 
a very small error. (A lot of theory of numerical integration goes into why 
numerical_integral works correctly, but that would take us too far afield to 
explain here.) 

Example 1.4.4. Let [...] denote the "greatest integer" (or "floor") function, so 
[1/2] = [0.5] = and [3/2] = [1.5] = 1. Evaluate J^[x]dx. (The function of 
y = [x] is sometimes called the "staircase function" because of the look of its 
discontinuous graph, Figure \1.S4\ ) 



Figure 1.34: Plot of the "greatest integer" function. 

Solution: f(x) = [x] is not continuous at x — 1 in the interval [1/2.3/2] so 
the Fundamental Theorem of Calculus can not be used. We can, however, use 
our understanding of the meaning of an integral as an area to get f^j^ [x] dx 

4 In fact, Sage includes Maxima (http://maxima.sf.net) and calls Maxima to compute 
this integral. 

5 In fact. Sage includes the GNU Scientific Library (http://www.gnu.org/software/gsl/) 
and calls it to approximate this integral. 
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= (area under y = between 0.5 and 1) + (area under y = 1 between 1 and 
1.5J=0 + 1/2 = 1/2. 

iVow, Zei's in/ something illegal - using the Fundamental Theorem of Calculus 
to evaluate this. Pretend for the moment that the Fundamental Theorem of 
Calculus is valid for discontinuous functions too. Let 



F(x) 



1, 1/2<x<1, 
x, 1 < x < 3/2. 



This function F is continuous and satisfies F'(x) = [x] for all x in [1/2,3/2] 
except x = 1 (where f(x) = [x] is discontinuous), so this F could be called an 
"antiderivative" of f . If we use it to evaluate the integral we get J^[x]dx = 

F{ x )\i^/2 = 3/2 — 1 = 1/2. This is correct. (Surprised?) Let's try another 
antiderivative. Let 



F[x) = { 2, 1/2<*<1, 
1 > ~ \ x, Kx< 3/2. 



This function F also satisfies F'(x) = [x] for all x in [1/2, 3/2] except x = 1. If 

3 /2 3 /2 

we use it to evaluate the integral we get [x] dx = F(x)\ 1 ^ 2 = 3/2 — 2 = —1/2. 
This doesn't even have the right sign (the integral of a non-negative function 
must be non-negative!), so it must be wrong. Moral of the story: In general, the 
Fundamental Theorem of Calculus is false for discontinuous functions. 



But does such an F as in the fundamental theorem of calculus (Theorem 1 1 .4. 1[ ) 
always exist? The surprising answer is "yes" . 

Theorem 1.4.2. Let F(x) = f* f(t)dt. Then F'(x) = f(x) for all x £ [a, b}. 

Note that a "nice formula" for F can be hard to find or even provably non- 
existent. 

The proof of Theorem 1.4.2 is somewhat complicated but is given in complete 
detail in many calculus books, and you should definitely (no pun intended) read 
and understand it. 

Proof: [Sketch of Proof] We use the definition of derivative. 

**(*) = lim F ^+ h )-n*) 

h^o h 
= lim ^jT + f(t)dt - jf f(t)<ttj /h 

= ^o{C hfm ) /h 
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Intuitively, for h sufficiently small / is essentially constant, so J x 1 f(t)dt ~ 
hf(x) (this can be made precise using the extreme value theorem). Thus 



which proves the theorem. □ 
1.4.2 Problems 

In problems 1 - 4 , rewrite the limit of each Riemann sum as a definite integral. 

1. lim mcsh (p)^o Efc=i( 2 + 3c k )Ax k on the interval [0,4]. 

2. lim mesh ( P )^ 5Zfe=i cos(5c fe )Aa; fc on the interval [0,11]. 
3- lim mcsh (p)^o YZ=i c l Ax k on the interval [2, 5]. 

4. lim mcsh(P) _ >0 YJk=i CftAarfc on the interval [2, 5]. 

5. Write as a definite integral (don't evaluate it though): The region bounded 
by y = x 3 , the rr-axis, the line x = 1, and x = 5. 

6. Write as a definite integral (don't evaluate it though): The region bounded 
by y = y/x, the x-axis, and the line x = 9. 

7. Write as a definite integral (do evaluate it, using geometry formulas): The 
region bounded by y = 2x, the x-axis, the line x = 1, and x = 3. 

8. Write as a definite integral (do evaluate it, using geometry formulas): The 
region bounded by y = \x\, the cc-axis, and the line x = — 1. 

9. For f(x) = 3 + x, partition the interval [0, 2] into n equally wide subinter- 
vals of length Ax = 2/n. 

(a) Write the lower sum for this function and partition, and calculate 
the limit of the lower sum as n — > oo. (b) Write the upper sum for this 
function and partition and find the limit of the upper sum as n — > oo. 

10. For f(x) = x 3 , partition the interval [0, 2] into n equally wide subintervals 
of length Ax = 2/n. 

(a) Write the lower sum for this function and partition, and calculate the 
limit of the lower sum as n — > oo. 

(b) Write the upper sum for this function and partition and find the limit 
of the upper sum as n — > oo. 
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1.4.3 Properties of the definite integral 

Definite integrals are defined as limits of Riemann sums, and they can be inter- 
preted as "areas" of geometric regions. This section continues to emphasize this 
geometric view of definite integrals and presents several properties of definite 
integrals. These properties are justified using the properties of summations and 
the definition of a definite integral as a Riemann sum, but they also have natural 
interpretations as properties of areas of regions. These properties are used in 
this section to help understand functions that are defined by integrals. They 
will be used in future sections to help calculate the values of definite integrals. 

Since integrals are a lot like sums (they are, after all, limits of them), their 
properties are similar too. Here is the integral analog of Theorem 1.3.11 

Theorem 1.4.3. (Integral Properties) 

• Integral of a constant function: cdx = c ■ (b — a). 

• Addition: J (f(x) + g{x)) dx = f(x) dx + g(x) dx. 

• Subtraction: f (f{x) — g(x)) dx = J f(x) dx — J g(x) dx. 

• Constant Multiple: J c ■ f{x) dx = c f{x) dx. 



Preserves positivity: If fix) > g{x) on for all x G [a,b], then 



f{x) dx > g{x) dx. 

J a 

In particular, if f(x) > on for all x G [a, b], then 

f{x)dx > 0. 

J a 

• Additivity of ranges: f{x) dx + f{x) dx = f{x) dx. 
Here are some other properties. 
Theorem 1.4.4. 

(b - a) ■ ( min f(x)) < [ f{x) dx<(b-a)-( mm f(x)). 

as£[a,6] J a x£[a,b] 

Which Functions Are Integrable? This important question was finally an- 
swered in the 1850s by Georg Riemann, a name that should be familiar by 
now. Riemann proved that a function must be badly discontinuous to not be 
integrable. 

Theorem 1.4.5. Every continuous function is integrable. If f is continuous on 
the interval [a,b], then lim mes h(p)_ + o(X)fc=i f( c k)A-Xk) is always the same finite 
number, namely, f(x)dx, so f is integrable on [a,b]. 
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Figure 1.35: Plot illustrating Theorem 11.4.4 



In fact, a function can even have any finite number of breaks and still be 
integrable. 

Theorem 1.4.6. Every bounded, piecewise continuous function is integrable. 
If f is defined and bounded ( for all x in [a,b], M < fix) < M for some 
AI > 0), and continuous except at a finite number of points in [a,b], then 
lim mcsh (p)^Q(^^ =1 /(cfe)Axfc) is always the same finite number, namely, f(x)dx, 
so f is integrable on [a, b] . 

Example 1.4.5. (A Nonintegrable Function) 

Though rarely encountered in "everyday practice", there are functions for 
which the limit of the Riemann sums does not exist, and those functions are 
not integrable. 

A nonintegrable function: The function 

. . J 1, if x is a rational number, 
1 0, if a; is an irrational number 

is not integrable on [0, 1] . 

Proof: For any partition P, suppose that you, a very rational (pun intended) 
person, always select values of Ck which are rational numbers. (Every subinterval 
contains rational numbers and irrational numbers, so you can always pick Ck to 
be a rational number.) Then /(cfc) = 1, and your Riemann sum is always 

n n 

Y P = /(c fe )A.Tfc = ^2 Ax k = x„ -x Q = l. 
fc=i fe=i 

Suppose your friend, however, always selects values of Ck which are irrational 
numbers. Then f(c.k) = 0, and your friend's Riemann sum is always 

n n 

fe=i fe=i 
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Now, take finer and finer partitions P so that mesh(P) — > 0. Keep in mind that, 
no matter how you refine P, you can always make "rational choices" for Ck and 
your friend can always make "irrational choices". We have limjneghfp^o Yp = 1 
and lim m esh(P)->o Fp = 0? so the limit of the Riemann sums doesn't have a 
unique value. Therefore the limit 



lim (Vf(c k )Ax k ) 

roesh(P)->0 ' 
k — 1 

does not exist, so f is not integrable. 



1.4.4 Problems 

Problems 1-20 refer to the graph of / in Figure 1.361 Use the graph to 
determine the values of the definite integrals. (The bold numbers represent the 
area of each region.) 




Figure 1.36: Plot for problems. 



1. J^f(x)dx 

2. fHf(x)dx 

3. J^f(x)dx 

4. $lf{x)dx 

5. J*f(x)dx 

6. J Q f(x)dx 

7. f s f(x)dx 

8. $lf{x)dx 

9. J°f(x)dx 
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10. f*f{x)dx 

11. f°f(x)dx 

12. f*2f(x)dx 

13. f*f(x) 2 dx 

14. r 3 i + /(t)* 

15. J Q 3 x + /(x) <ix 

16. J 3 5 3 + /(x)dx 

17. J 5 2 + /(x)dx 

18. J*\f(x)\dx 

19. / 7 3 l + |/(x)|dx 

For problems 21-28, sketch the graph of the integrand function and use it to 
help evaluate the integral. (|...| denotes the absolute value and [...] denotes the 
integer part.) 

21. ; 4 |x| dx, 

22. J 4 l + |x|dx, 

23. ^ x \x\dx, 

24. \l |x| - 1 dx, 

25. ff[x]dx, 

26. jf 5 [x]cfa;, 

27. Jf 2 + [x] dx, 

28. J^dx. 
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1.5 Areas, integrals, and antiderivatives 

This section explores properties of functions defined as areas and examines some 
of the connections among areas, integrals and antiderivatives. In order to focus 
on the geometric meaning and connections, all of the functions in this section are 
nonncgative, but the results are generalized in the next section and proved true 
for all continuous functions. This section also introduces examples to illustrate 
how areas, integrals and antiderivatives can be used. When / is a continuous, 
nonnegative function, then the "area function" A(x) = J f(t) dt represents the 
area between the graph of /, the i-axis, and between the vertical lines at t = a 
and t = x (Figure fl.37j ). and the derivative of A{x) represents the rate of change 
(growth) of A(x). 




Let F(x) be a differentiable function. Call F(x) an antiderivative of f(x) if 
4^F(x) = f(x). We have seen examples which showed that, at least for some 
functions /, the derivative of A(x) was equal to / so A(x) was an antiderivative 
of /. The next theorem says the result is true for every continuous, nonnegative 
function /. 

Theorem 1.5.1. ("The Area Function is an Antiderivative") If f is a contin- 
uous nonnegative function, x > a, and A(x) — J x fit) dt then ^ f(t) dt = 
■4-Aix) = f(x), so A(x) is an antiderivative of fix). 

This result relating integrals and antiderivatives is a special case (for non- 
negative functions /) of the Fundamental Theorem of Calculus. This result is 
important for two reasons: 

• it says that a large collection of functions have antiderivatives, and 

• it leads to an easy way of exactly evaluating definite integrals. 

x 

Example 1.5.1. Let G(x) = 4^ cos(t)dt. Evaluate G(x) for x = 7r/4, tt/2, 
and 37r/4. 

Solution: It is not hard to plot the graph of Aix) = J Q cos(t)dt = sin(a;) 
(Figure [L3Bl) . By the theorem, A'{x) = G{x) = cos(x) so A'{-k/4) = cos(7r/4) = 
.707..., A'(tt/2) = cos(tt/2) = 0, and A'(3tt/4) = cos(3/4) = -0.707... . 
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Here is the plot of y = A(x) and y = G{x): 

Incidentally, this was created using the following Sage commands. 



sage: P = plot(cos(x) ,x,0,2*pi,linestyle=" — ") 

sage: Q = plot(sin(x) ,x,0,2*pi) 

sage: R = text ("$y=A(x) = \sin(x) $" , (3 . 1 , 1) ) 

sage: S = text ("$y=G(x) = \cos(x)$" , (6 .8,0 .7) ) 

sage: show(P+Q+R+S) 



Theorem 1.5.2. ( "Antiderivatives and Definite Integrals") If f is a continu- 
ous, nonnegative function and F is any antiderivative of f (F'(x) = f(x)) on 
the interval [a,b], then 



area bounded between the graph 

of f and the x-axis and 
vertical lines at x = a and x = b 



= J b a f(x)dx = F(b)-F(a). 



The problem of finding the exact value of a definite integral reduces to finding 
some (any) antiderivative F of the integrand and then evaluating F(b)F(a). 
Even finding one antiderivative can be difficult, and, for now, we will stick to 
functions which have easy antidcrivatives. Later we will explore some methods 
for finding antidcrivatives of more difficult functions. 

The evaluation F(b) — F(a) is represented by the symbol F(x)\^. 

Example 1.5.2. Evaluate xdx in two ways: 

(a) By sketching the graph of y = x and geometrically finding the area. 

(b) By finding an antiderivative of F(x) of f and evaluating F(3) — F(l). 
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Solution: (a) The graph of y = x is a straight line, so the area is a triangle 
which geometrical formulas (area=^bh) tell us has area 4. 
(b) One antiderivative of x is F{x) — \x 2 (check that j^.{\x 2 ) = x), and 

F(x)\l = F(3)-F(l) = ^3 2 -il 2 =4, 

which agrees with (a). Suppose someone chose another antiderivative of x, say 
F(x) = \x 2 + 7 (check that j^{\x 2 + 7) = x), then 

F(x)\f = F(3) - F(l) = C-Z 2 + 7) - {\l 2 + 7) = 4. 
No matter which antiderivative F is chosen, F(3) — F(l) equals 4. 
Practice 1.5.1. Evaluate j^(x—l)dx in the two ways of the previous example. 

Practice 1.5.2. Find the area between the graph of y = 3x 2 and the horizontal 
axis for x between 1 and 2. 

1.5.1 Integrals, Antiderivatives, and Applications 

The antiderivative method of evaluating definite integrals can also be used when 
we need to find an "area," and it is useful for solving applied problems. 

Example 1.5.3. Suppose that t minutes after putting 1000 bacteria on a petri 
plate the rate of growth of the population is 6t bacteria per minute. 

(a) How many new bacteria are added to the population during the first 7 
minutes? 

(b) What is the total population after 7 minutes? 

(c) When will the total population be 2200 bacteria? 

Solution: (a) The number of new bacteria is the area under the rate of growth 
graph, and one antiderivative of 6i is 3t 2 (check that j^{3t 2 ) = 6t) so new 

bacteria = J^6tdt = 3t 2 \ 7 Q = 147. 

(b) The new population = (old population) + (new bacteria) = 1000 + 147 = 
1147 bacteria. 

(c) If the total population is 2200 bacteria, then there are 2200 — 1000 = 1200 
new bacteria, and we need to find the time T needed for that many new bacteria 
to occur. 1200 new bacteria = J^6tdt = 3t 2 |J = 3(T) 2 - 3(0) 2 = 3T 2 so 
T 2 = 400 and T = 20 minutes. After 20 minutes, the total bacteria population 
will be 1000 + 1200 = 2200. 

Practice 1.5.3. A robot has been programmed so that when it starts to move, 
its velocity after t seconds will be 3t 2 feet/second. 

(a) How far will the robot travel during its first 4 seconds of movement? 
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(b) How far will the robot travel during its next 4 seconds of movement? 

(c) How many seconds before the robot is 729 feet from its starting place? 
(Hint: an antiderivative of it 2 ist 3 .) 

Practice 1.5.4. The velocity of a car after t seconds is 2t feet per second. 

(a) How far does the car travel during its first 10 seconds? 

(b ) How many seconds does it take the car to travel half the distance in part 



1.5.2 Indefinite Integrals and net change 

We've seen how integrals can be interpreted using area. In this section, we will 
see how integrals can be interpreted physically as the "net change" of a quantity. 

The notation J f(x)dx = F{x) means that F'(x) = f(x) on some (usually 
specified) domain of definition of f{x). Recall, we call such an F(x) an an- 
tiderivative of fix). 

Proposition 1.5.1. Suppose f is a continuous function on an interval (a,b). 
Then any two antiderivatives differ by a constant. 

Proof: If Fx{x) and F 2 (x) are both antiderivatives of a function fix), then 



Thus Fi(x) — F2{x) = c from some constant c (since only constant functions 
have slope everywhere). Thus F\{x) — F 2 (x) + c as claimed. □ 
We thus often write 



where C is an unspecified constant. 

Note that the proposition need not be true if / is not defined on a whole 
interval. For example, f{x) — 1/x is not defined at 0. For any pair of constants 
ci, C2, the function 



satisfies F'(x) = f(x) for all i^O. We often still just write J 1/x = ln(|x|) + c 
anyways, meaning that this formula is supposed to hold only on one of the 
intervals on which l/x is defined (e.g., on (— oo,0) or (0,oo)). 

We pause to emphasize the notation difference between definite and indefinite 
integration. 



(a)? 



(F^x) - F 2 (x)Y = F[{x) - F^x) = f{x) - f{x) = 0. 





ln(|x|) + ci x < 0, 
ln(x) + C2 x > 0, 




a (family of) functions 



a specific number 
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There are no small families in the world of antiderivatives: if / has one an- 
tiderivative F (as it always does, unless / is a really unusual function), then / 
has an infinite number of antiderivatives and every one of them has the form 
F(x) + C. 

Example 1.5.4. There are many ways to write a particular indefinite inte- 
gral and some of them may look very different. You can check that F(x) = 
sin(a;) 2 , G(x) = — cos(x) 2 , and H(x) = 2sin(a;) 2 + cos(x) 2 all have the same 
derivative f(x) = 2 sin(a;) cos(x), so the indefinite integral of 2 sin(a;) cos(x), 
j 2 sin(ir) cos(a;) dx, can be written in several ways: sin(a;) 2 +C 7 or — cos(x) 2 +C, 
or 2 sin(x) 2 + cos(x) 2 + C. 

One of the main goals of this course is to help you to get really good at 
computing J f(x)dx for various functions f(x). It is useful to memorize a table 
of examples, such as the one in the appendix below, since often the trick to 
integration is to relate a given integral to one you know. Integration is like 
solving a puzzle or playing a game, and often you win by moving into a position 
where you know how to defeat your opponent, e.g., relating your integral to 
integrals that you already know how to do. If you know how to do a basic 
collection of integrals, it will be easier for you to see how to get to a known 
integral from an unknown one. 

Whenever you successfully compute F(x) = J f(x)dx 7 then you've constructed 
a mathematical gadget that allows you to very quickly compute f{x)dx for 
any a, b (in the interval of definition of f{x)). The gadget is F(b) — F(a). This 
is really powerful. 

Example 1.5.5. 

x 2 + 1 H = -dx = I x z dx + / ldx + I —= -dx 



i.T 2 + x + tan 1 (x) 



Example 1.5.6. 



Example 1.5.7. 



J ^dx = J V5x~ 1/2 dx = 2VE: 



x l ' 2 + c. 



sin(2x) , f 2 siiha;) cos(x) 

K 'dx= / \ ' v ' = / 2cos(x) = 2sin(x) + c 



sin(x) J sin(x) 

Particular Antiderivatives: You can verify the following yourself. 
• Constant Function: J k dx — kx + C 



• Powers of x: f x" dx — + C, 
n ^ -1. 



fx- 1 dx = \n(x) + C. 
Common special cases: 
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- J^dx^ |.T 3 / 2 +C. 

- f i dx = 2X 1 ' 2 + C. 



• Trigonometric Functions: 
/ cos(ax) dx = i sin(x) + C. 
J sin(ax) dx = — ~ cos(:r) + C. 
J sec(ax) 2 dx = - tan(x) + C . 
J csc(ax) 2 dx = — - cot (a) + C. 
J sec(ax) tan(x) dx = ^ sec(a;) + C. 
J csc(ax) cot (a;) dx = — - csc(x) + C. 



Common special cases: 

— / cos(cc) dx = sin (a;) + C. 

— J sin(x) dx = — cos(a;) + C . 

1.5.3 Physical Intuition 

In the previous lecture we mentioned a relation between velocity, distance, and 
the meaning of integration, which gave you a physical way of thinking about 
integration. In this section we generalize our previous observation. 
The following is a restatement of the fundamental theorem of calculus. 

Theorem 1.5.3. (Net Change Theorem) The definite integral of the rate of 
change f'{x) of some quantity f(x) is the net change in that quantity: 



For example, if p(t) is the population of your hometown at time t, then p'(t) 
is the rate of change. If p'(t) is positive then your hometown is growing. The 
net change interpretation of integration is that 



p (t) dt = p(t2)— p(t±) = change in number of residents from time ti to t 2 . 



Another very common example you'll seen in problems involves water flow into 
or out of something. If the volume of water in your bathtub is V(t) gallons at 
time t (in seconds), then the rate at which your tub is draining is V'(t) gallons 
per second. If you have the geekiest drain imaginable, it prints out the drainage 
rate V'(t). You can use that printout to determine how much water drained out 
from time t± to t-2,: 





ti 




ti 



V (t) dt = water that drained out from time t\ to t 2 
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Some problems will try to confuse you with different notions of change. A 
standard example is that if a car has velocity v(t), and you drive forward, then 
slam it in reverse and drive backward to where you start (say 10 seconds total 
elapse), then v{t) is positive some of the time and negative some of the time. 
The integral J Q v(t)dt is not the total distance registered on your odometer, 
since v(t) is partly positive and partly negative. If you want to express how far 
you actually drove going back and forth, compute J„ |w(t)|c?t. The following 
example emphasizes this distinction: 

Example 1.5.8. A bug is pacing back and forth, and has velocity v(t) = 
t 2 - It - 8. Find 

(1) the displacement of the bug from time t = 1 until time t = 6 (i.e., how 
far the bug is at time 6 from where it was at time 1), and 

(2) the total distance the bug paced from time t = 1 to t = 6. 
For (1), we compute 

-t 3 -t 2 -8t = . 

3 Ji 3 

For (2), we compute the integral of\v(t)\: 
f 6 

I \t 2 -2t-8\dt = 



(t 2 - 2t 



\dt 



- ri - r 



8f 



1 u 



-t A - t z 



8t 



= 18 
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1.5.4 Problems 

1. Two objects start from the same location and travel along the same path 
with velocities VA{t) = t + 3 and Us(t) = t 2 At + 3 meters per second. How 
far ahead is A after 3 seconds? After 5 seconds? 



2. Sketch the graph of each function and find the area between the graphs 
of f(x) = x 2 + 3, g(x) = 1 and -1 < x < 2. 

3. Sketch the graph of each function and find the area between the graphs 
of f(x) = x 2 + 3, g(x) = 1 + x and < x < 3. 

4. Sketch the graph of each function and find the area between the graphs 
of f(x) = x 2 , g(x) = x and < x < 2. 

5. Sketch the graph of each function and find the area between the graphs 
of f(x) = x + 1, g(x) = cos(x) and < x < ir/4. 

6. If f(t) denoted the velocity of a bug traveling along a line at time t, find 
the distance traveled in the first 4 seconds. 

7. If f(t) denoted the velocity of a bug traveling along a line at time t, find 
the distance traveled in the first 4 seconds. 
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Figure 1.39: Velocity of bug at time t. 




Figure 1.40: Velocity of bug at time t. 



1.6 Substitution and Symmetry 



Remarks: 



1. The total distance traveled is \v(t)\dt since \v(t)\ is the rate of 
change of F(t) = distance traveled (your speedometer displays the rate of 
change of your odometer). 



2. How to compute J \f(x)\dx. 



(a) Find the zeros of f(x) on [a, b], and use these to break the interval 
up into subintcrvals on which fix) is always > or always < 0. 

(b) On the intervals where f{x) > 0, compute the integral of /, and on 
the intervals where f(x) < 0, compute the integral of — /. 

(c) The sum of the above integrals on intervals is J \ f(x)\dx. 



This section is primarily about a powerful technique for computing definite 
and indefinite integrals. 
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1.6.1 The Substitution Rule 

In first quarter calculus you learned numerous methods for computing deriva- 
tives of functions. For example, the power rule asserts that 

(x a )' = a-x a -\ 

We can turn this into a way to compute certain integrals: 

x a dx = -^—x a+1 if a + -1. 
a + 1 r 

Just as with the power rule, many other rules and results that you already 
know yield techniques for integration. In general integration is potentially much 
trickier than differentiation, because it is often not obvious which technique to 
use, or even how to use it. Integration is a more exciting than differentiation! 

Recall the chain rule, which asserts that 

^f(g(x))=f(g(x))g'(x). 
dx 

We turn this into a technique for integration as follows: 
Proposition 1.6.1. (Substitution Rule) Let u — g(x), we have 

f{g{x))g'{x)dx = J f{u)du, 

assuming that g(x) is a function that is differentiable and whose range is an 
interval on which f is continuous. 

Proof: Since / is continuous on the range of g, Theorem 1.4.2 (the funda- 
mental theorem of Calculus) implies that there is a function F such that F' = f. 
Then 



f(g(x))g'(x)dx = I F'(g(x))g'(x)dx 
d 



lx F(g(x)))dx 
F(g(x)) + C 
F{u) + C= I F'{u)du = I f{u)du. 



□ 

If u = g(x) then du = g'(x)dx, and the substitution rule simply says if you let 
u = g(x) formally in the integral everywhere, what you naturally would hope to 
be true based on the notation actually is true. The substitution rule illustrates 
how the notation Leibniz invented for Calculus is incredibly brilliant. It is said 
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that Leibniz would often spend days just trying to find the right notation for a 
concept. He succeeded. 

As with all of Calculus, the best way to start to get your head around a new 
concept is to see severally clearly worked out examples. (And the best way to 
actually be able to use the new idea is to do lots of problems yourself!) In this 
section we present examples that illustrate how to apply the substitution rule 
to compute indefinite integrals. 

Example 1.6.1. 

r x 2 (x 3 +5fdx 



Let u — x 3 + 5. Then du = 2>x 2 dx, hence dx = du/(3x 2 ). Now substitute it all 
in: 

x 2 (x 3 + 5fdx = [ -u 9 = — U 10 = -(x 3 + 5) 10 . 
v ' J 3 30 30 v ; 

There's no point in expanding this out: "only simplify for a purpose/" 
Example 1.6.2. 

-dx 



l + e 1 

Substitute u = 1 + e x . Then du = e x dx, and the integral above becomes 



[ — =lnM =ln\l + e x \ = ln(l + e*). 
J u 



Note that the absolute values are not needed, since 1 + e x > for all x. 
Example 1.6.3. 



/ 



x 2 

-dx 



VT^x 

Keeping in mind the power rule, we make the substitution u = 1 — x. Then 
du = —dx. Noting that x = 1 — u by solving for x in u = 1 — x, we see that the 
above integral becomes 

(l-u) 2 , f l-2u + u 2 , 
-du=- I du 

u -l/2_ 2u l/2 +u 3/2 du 



2u Va _ 4 u 3/2 + 2 5/2 
3 5 

= -2(1 - x) 1 ' 2 + 4 (1 - x) 3 ' 2 - 2 (1 - xf' 2 . 
3 5 

The steps of the "change of variable" method can be summarized as 

1. set a new variable, say u , equal to some function of the original variable 
x (usually u is set equal to some part of the original integrand function) , 
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2. calculate the differential du as a function of dx, 

3. rewrite the original integral in terms of u and du, 

4. integrate the new integral to get an answer in terms of u, 

5. replace the u in the answer to get an answer in terms of the original 
variable. 

A "Rule of thumb" for changing the variable: If part of the integrand is 
a composition of functions, f(g(x)), then try setting u = g(x), the "inner" 
function. 

Example 1.6.4. elect a function for u for each integral and rewrite the integral 
in terms of u and du. 

(a) J cos(3:e)a/2 + sin(3x) dx, (b) J^^dx, (c) j e x sm(e x ) dx. 
Solution: (a) Put u = 2 + sin(3a;). Then du = 3cos(3a;) dx, and the integral 

becomes J ^^/udu. 

(b) Put u = 2 + e x . Then du = e x dx, and the integral becomes J - du. 

(c) Put u = e x . Then du = e x dx, and the integral becomes J sin(it) du. 



1.6.2 Substitution and definite integrals 

Once an antiderivative in terms of u is found, we have a choice of methods. We 
can 

(a) rewrite our antiderivative in terms of the original variable x, and then 
evaluate the antiderivative at the integration endpoints and subtract, or 

(b) change the integration endpoints to values of u, and evaluate the an- 
tiderivative in terms of u before subtracting. 

If the original integral had endpoints x = a and x — b, and we make the sub- 
stitution u — g{x) and du = g'(x)dx, then the new integral will have endpoints 
u = g(a) and u = g(b): 



ru=g(b) 

(original integrand) dx becomes / (new integrand) du. 

J u—aia) 



t^x—b pu—g{b) 

'u=g(a) 

Example 1.6.5. To evaluate 



l 

(3x — l) 4 dx, 



we can, in line with the "Rule of thumb " above, use the substitution u = 3x — 1 . 
Then du = 4^ (3a; — l)dx = "Sdx, so the indefinite integral J (3x — l) 4 dx becomes 
Jlu 4 du= ±u 5 + C. 

(a) Converting our antiderivative back to the variable x and evaluating with 
the original endpoints: 
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'0 



1 1 32 —1 11 

(3* - I) 4 dx = ( - (3a: - l) 5 + C% = - — = - = 2.2. 



(b) Converting the integration endpoints to values of u : when x = 0, then 
— 3x — 1 = 3- 0— 1 = —1, and when x = 1, then u = 3x — 1 = 3-1 — 1 = 2 so 

(3x - l) 4 dx = [ \u 4 du = (-U 5 + C)\\ = ^ = 2.2. 
o J-i 3 15 5 

Both approaches typically involve about the same amount of work and calcula- 
tion. Of course, these approaches lead to the same numerical answer, by the 
"substitution rule" ( Proposition \ 1 . 67I\ ) . 

Here's how to do this using Sage. Note that the area under the two curves 
plotted below, y = (3x — l) 4 , < x < 1, and y = x 4 /3, —1 < x < 2, are the 
same. 



sage: x,u = var("x,u") 

sage: integral ( (3*x-l) ~4,x, , 1) 

11/5 

sage: integral (u~4/3,u, -1,2) 
11/5 

sage: P = plot ( (3*x-l) ~4,x, , 1 ,rgbcolor=(0 . 7 , . 1 ,0 . 5) , plot_points=40) 

sage: Q = plot (u"4/3 ,u, -1 , 2 , linestyle=" : " ) 

sage: R = text("$y=(3x-l)~4$" , (1.4,12)) 

sage: S = text (" $y=x"4/3$ ", (2,2.5)) 

sage: plot (P+Q+R+S) 




2.5 



-2.5 f- 

Figurc 1.41: Plots of y — (3x — l) 4 and y = x 4 /3. 



1.6.3 Symmetry 

An odd function is a function f(x) (defined for all reals) such that f(—x) = 
—f(x), and an even function one for which f(—x) = f(x). If / is an odd 
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function, then for any a, 

f f(x)dx = 0. 

J —a 

If / is an even function, then for any a, 

f(x)dx = 2 [ f(x)dx. 



Both statements are clear if we view integrals as computing the signed area 
between the graph of f(x) and the x-axis. 



Example 1.6.6. An even example, 



1 f 1 

x 2 dx = 2 / x dx = 2 

l Jo 



— x 
3 



and an odd example, 



j x 3 dx = 


[K 


i 

= 









These computations are consistent with the symmetry (or "anti-symmetry") 
of the graphs and what you know about the relationship between the integral and 
area. 




-0.25 
-0.5 

-0.75 



Figure 1.42: Plots of y = x and y = x . 



1.6.4 Problems 

For the problems below, let f(x) = x 2 and g(x) = x and verify that 
1- / f(x) ■ g(x) dx ^ J f(x) dx ■ J g(x) dx. 

2. / f(x)/g(x) dx^J f{x) dx/ J g(x) dx. 

3. ^/jf(x)dx^J^f{x)dx. 
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J f(x)dx r J f( x ) 

5. f cos (3a;) dx, u = 3x. 

6. J" sin(7a;) dx, u = 7x. 

7. J e 5x dx, u = 5x. 

8. J e 3x + cos(2cc) dx, u = 3x and u = 2x. 
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Chapter 2 

Applications 



2.1 Applications of the integral to area 

The development of calculus by Newton and Leibniz was a vital step in the 
advancement of pure mathematics, but Newton also advanced the applied sci- 
ences and mathematics. Not only did he discover theoretical results, but he 
immediately used those results to answer important applied questions about 
gravity and motion. The success of these applications of mathematics to the 
physical sciences helped establish what we now take for granted: mathematics 
can and should be used to answer questions about the world. Newton applied 
mathematics to the outstanding problems of his day, problems primarily in the 
field of physics. In the intervening 300 years, thousands of people have contin- 
ued these theoretical and applied traditions and have used mathematics to help 
develop our understanding of all of the physical and biological sciences as well 
as the behavioral sciences and business. Mathematics is still used to answer 
new questions in physics and engineering, but it is also important for modeling 
ecological processes, for understanding the behavior of DNA, for determining 
how the brain works, and even for devising strategies for voting effectively. The 
mathematics you are learning now can help you become part of this tradition, 
and you might even use it to add to our understanding of different areas of life. 
It is important to understand the successful applications of integration in case 
you need to use those particular applications. It is also important that you 
understand the process of building models with integrals so you can apply it to 
new problems. Conceptually, converting an applied problem to a Ricmann sum 
is the most valuable and the most difficult step. 

2.1.1 Using integration to determine areas 

This section is about how to compute the area of fairly general regions in the 
plane. Regions are often described as the area enclosed by the graphs of sev- 
eral curves. ("My land is the plot enclosed by that river, that fence, and the 
highway.") 
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Recall that the integral J f(x)dx has a geometric interpretation as the signed 
area between the graph of f(x) and the cc-axis. We defined area by subdividing, 
adding up approximate areas (use points in the intervals) as Riemann sums, 
and taking the limit. Thus we defined area as a limit of Riemann sums. The 
fundamental theorem of calculus asserts that we can compute areas exactly 
when we can finding antidcrivatives. 



Ax 




Figure 2.1: Area between y = f(x) and y = g(x). 



Instead of considering the area between the graph of f(x) and the x-axis, 
we consider more generally two graphs, y = f(x), y = g{x), and assume for 
simplicity that f(x) > g(x) on an interval [a, b]. Again, we approximate the area 
between these two curves as before using Riemann sums. Each approximating 
rectangle has width (b — a)/n and height f(x) — g(x), so 

Area bounded by graphs ~ J~][f(cj) ~ g( c i)]Ax. 

Note that f(x) — g(x) > 0, so the area is nonnegative. From the definition of 
integral we see that the exact area is 

,6 

Area bounded by graphs = / (f(x) — g(x))dx. (2.1) 

J a 

Why did we make a big deal about approximations instead of just writing 
down (2.1) ? Because having a sense of how this area comes directly from a 
Riemann sum is very important. But, what is the point of the Riemann sum 
if all we're going to do is write down the integral? The sum embodies the 
geometric manifestation of the integral. If you have this picture in your mind, 
then the Riemann sum has done its job. If you understand this, you're more 
likely to know what integral to write down; if you don't, then you might not. 

Remark 2.1.1. By the linearity property of integration, our sought for area is 
the integral of the "top" function minus the integral of the "bottom" function, 



b 

f(x)dx — / g(x)dx, 



of two signed areas. 
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Example 2.1.1. Find the area enclosed by y = x + 1, y = 9 — x 2 , x = — 1. 

x = 2 (see Figure [EM . 




Figure 2.2: Plots of y = x + 1 and y = 9 — a; 2 . 



j4rea = 



(9-x 2 ) - (x+ 1) 



W^e /iawe reduced the problem to a computation: 
[(9 - x 2 ) - + l)]dx = J (8-x-x 2 )dx 
Here is this plot and computation in Sage: 



dx . 



8x x x* 



39 
Y' 



sage: x = var("x") 
sage: PI = plot(x+l, x, -2, 3) 

oairo' P9 = nl nt f Q - sr~0 v 



CD I * — 3 — » — j — ✓ 

age: P2 = plot (9 - x~2, x, -2, 3) 

age: Tl = text("$y = x+l$", (1,2.6)) 

,age: T2 = text("$y = 9-x~2$", (2,7)) 



— o — " ~<- J — 

sage: show(Pl+P2+Tl+T2) 

sage: integrate ( (9-x~2) 
39/2 



(x+l),x, -1, 2) 



The above example illustrates the simplest case. In practice more interesting 
situations often arise. The next example illustrates finding the boundary points 
a, b when they are not explicitly given. 

Example 2.1.2. Find area enclosed by the two parabolas y = 12 — x 2 and 
y = x 2 — 6. 

Problem: We didn't tell you what the boundary points a, b are. We have to figure 
that out. How? We must find exactly where the two curves intersect, by setting 
the two curves equal and finding the solution. We have 
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-15 - 



Figure 2.3: Plots of y = 12 - x 



2 and y = x 2 — 6. 



x 2 - 6 = 12 - x 2 , 



so = 2x 2 — 18 = 2(x 2 — 9) = 2(x — 3)(a; + 3), hence the intersect points are at 
a = — 3 and 6 = 3. FFe t/ius find the area by computing 



[12 - x 2 - {x 2 - 6)] dx = / (18 - 2.T 2 )dx = 4 / (9 - x 2 )dx = 4 • 18 = 72. 



.ffere is i/iis p/oi anrf computation in Sage: 

sage: PI = plot(12-x~2, x, -5, 5) 

sage: P2 = plot(x~2-6, x, -5, 5) 

sage: Tl = text("$y = 12-x~2$", (-3.5,-10)) 

sage: T2 = text("$y = x~2-6$", (2,-7)) 

sage: show(Pl+P2+Tl+T2) 

sage: integrate( (12-x~2) - (x~2-6),x, -3, 3) 
72 



Of course, if you had mistakenly computed J_ 3 \_(x 2 — 6) — (12 — a; 2 )] dx, then 
don't worry. You would 've gotten —72 as your answer. However, always re- 
member areas are non-negative, so the correct answer is 72. 

Example 2.1.3. A common way in which you might be tested to see if you 
really understand what is going on, is to be asked to find the area between two 
graphs x = f{y) and x = g(y). If the two graphs are vertical, subtract off 
the right-most curve. Or, just "switch x and y" everywhere (i.e., reflect about 
y = x). The area is unchanged. 
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For instance, consider the area between the two parabolas X = 12 — y 2 and 
x = y 2 — 6. 



si- 




Figure 2.4: Plots of x = 12 — y 2 and x = y 2 — 6. 

Swapping x and y amounts to reflecting the plot in Figure \2.4\ above about the 
45° line y = x. The reflected graph coincides with that in Figure \2.3\ above. 
Therefore, by Example 2. 1.2. the answer is 72. 

Example 2.1.4. Find the area (not signed area!,) enclosed by y = sin(7rx), 
y = x 2 — x, and x = 2. 



3 






\ 2 






\ 1 


y=x -x / 




-1.5 -i 0.5 / 


wJ|!L— — "% Lb / 


> 2.5 


Figure 2.5: Plots of 


y = sin(7ra;) and y = 


X — X. 



Write x 2 — x = (x— 1/2) 2 — 1 /4, so that we can obtain the graph of the parabola 
by shifting the standard graph. The area comes in two pieces, and the upper and 
lower curve switch in the middle. Technically, what we're doing is integrating 
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the absolute value of the difference. The area is 

[ silica;) — (x 2 — x)dx — [ (x 2 — x) — sm(irx)dx = — hi 

Jo Jl 7T 

Here is this plot and computation in Sage: 



sage: PI = plot (sin(pi*x) , x, -1, 2.2) 

sage: P2 = plot(x"2-x, x, -1, 2.2) 

sage: P3 = list_plot( [(2,0) , (2,2)] ,plotjoined=True,linestyle=" : ") 

sage: Tl = text("$y = \sin(\pi x)$", (-1.2,-1)) 

sage: T2 = text("$y = x~2-x$", (1.3,1.3)) 

sage: show(Pl+P2+P3+Tl+T2) 



Something to take away from this is that in order to solve this sort of problem, 
you need some facility with graphing functions. If you aren't comfortable with 
this, review. 

2.2 Computing Volumes of Surfaces of Revolu- 
tion 

The last section emphasized a geometric interpretation of definite integrals as 
"areas" in two dimensions. This section emphasizes another geometrical use 
of integration, calculating volumes of solid threedimensional objects, such as a 
volume of revolution. Our basic approach is to cut the whole solid into thin 
"slices" whose volumes can be approximated, add the volumes of these "slices" 
together (a Ricmann sum) , and finally obtain an exact answer by taking a limit 
of the sums to get a definite integral. 

Practice 2.2.1. Most people have a body density between 0.95 and 1.05 times 
the density of water which is 62.5 pounds per cubic foot. Use your weight (in 
lbs) to estimate the volume of your body (in cubic ft). (If you float in fresh 
water, your body density is less than 1.) 

First, we introduce the building blocks of this section, right solids. A right 
solid is a three-dimensional shape swept out by moving a planar region A some 
distance h along a line perpendicular to the plane of A. For instance, if A 
is a rectangle, then the "right solid" formed by moving A along the line is a 
3dimcnsional solid box B and, of course, the volume of B is 



(area of A) x (distance along the line) = (base) x (height) x (width). 

The region A is called a face of the solid. The word "right" is simply used to 
indicate that the movement is along a line perpendicular (at a right angle) to 
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the plane of A. Two parallel cuts though the shape produce a slice with two 
faces. 

Example 2.2.1. Suppose there is a fine, uniform mist in the air, and every 
cubic foot of mist contains 0.02 ounces of water droplets. If you run 50 feet in 
a straight line through this mist, how wet do you get? Assume that the front (or 
a cross section) of your body has an area of 8 square feet. 

Solution: As you run, the front of your body sweeps out a "tunnel" through 
the mist. The volume of the tunnel is the (cross sectional) area of the front of 
your body multiplied by the length of the tunnel: volume = (8 ft) (50 ft) = 400 
ft. Since each cubic foot of mist held 0.02 ounces of water which is now on you, 
you swept out a total of (400 ft )(0.02 oz/ft) = 8 ounces of water. 

A general solid can be cut into slices which are almost right solids. An indi- 
vidual slice may not be exactly a right solid since its cross sections may have 
different areas. However, if the cuts are close together, then the cross sectional 
areas will not change much within a single slice. Each slice will be almost a 
right solid, and its volume will be almost the volume of a right solid. 

Suppose an xaxis is positioned below the solid shape, and let A(x) be the 
area of the face formed when the solid is cut at x perpendicular to the xaxis. 
If P = {x = b} is a partition of [a, b], and the solid is cut at 

each Xi, then each slice of the solid is almost a right solid, and the volume of 
each slice is approximately 

(area of a face of the slice) x (thickness of the slice) = A(xi)Ax.i. 

The total volume V of the solid is approximately the sum of the volumes of the 
slices: 

V = volume of each slice = A(xi)Axi, 

i 

which is a Riemann sum. The limit, as the mesh of the partition approaches 
(taking thinner and thinner slices), of the Riemann sum is the definite integral 
oiA{x): 

V = M x i)&Xi -> / A(x) dx. 

Theorem 2.2.1. (Volume By Slices Formula) If S is a solid and A(x) is the 
area of the face formed by a cut at x and perpendicular to the x-axis, then the 
volume V of the part of S above the interval [a,b] is V = J^A(x)dx. 

Everybody knows that the volume of a solid box is 

volume = length x width x height. 

More generally, the volume of cylinder is V = irr 2 h (cross sectional area times 
height). Even more generally, if the base of a prism has area A, the volume of 
the prism is V — Ah. 
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But what if our solid object looks like a complicated blob? How would we 
compute the volume? We'll do something that by now should seem familiar, 
which is to chop the object into small pieces and take the limit of approxima- 
tions. If these small pieces are cross sections then the corresponding method of 
computing the volume of revolution is called the "disc method" . If these small 
pieces are cylindrical shells then the corresponding method of computing the 
volume of revolution is called the "shell method". We look in detail into the 
disc method first, followed by the shell method. 

2.2.1 Disc method 

Assume that we have a function 

A(x) — cross sectional area at x. 

The volume of our potentially complicated blob is approximately A(xi)Ax. 
Thus 

n 

volume of blob = lim > A{xi)Ax 

i=l 

b 

A{x)dx 

Here is the plot a picture of solid sliced vertically into a bunch of vertical thin 
solid discs: 





Figure 2.6: Plot of the cone z = |(1- 
are approximated by thin discs. 



\/ x 2 + y 2 ) sliced into thin "shells" , which 



Example 2.2.2. Find the volume of the pyramid with height H and square base 
with sides of length L. 

For convenience look at pyramid on its side, with the tip of the pyramid at the 
origin. We need to figure out the cross sectional area as a function of x, for 
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< x < H . The function that gives the distance s{x) from the x-axis to the 
edge is a line, with s(0) = and s(H) = L/2. The equation of this line is thus 



s(x) 



2H' 



Thus th 



e cross s 



ectional 



The volume is then 



A(x)dx 



A{x) = (2s(x)) 2 = 



x 2 L 2 



H x 2 L 2 
H 2 



dx 



x A L 2 



3H 2 



H 



~W 2 



= 3 HL 




Figure 2.7: How big is Pharaoh's place?. (Photo found on 

|http : //en. wikipedia . org/ wiki/Egypt ian_pyramids[ taken by Ricardo 
Libcrato.) 



When a region is revolved around a line (Figure 2.8j) a right solid is formed. 
When the face of each slice of the revolved region is a circle then the formula 
for the area of the face is easy: A(x) = area of a circle = 7r(radius), where the 
radius is often a function of the location x. Finding a formula for the changing 
radius requires care. 

Theorem 2.2.2. (Volumes of Revolved Regions by Discs) If the region formed 
between f , the horizontal line y = L, and the interval [a, b] is revolved about the 
horizontal line y = L (see Figure [278\) then the volume is V = J A(x) dx = 
7T (radius ) dx = J ir(f(x) — L) dx. 

Example 2.2.3. Find the volume of the solid obtained by rotating the following 
"flower pot" region about the x axis: the region enclosed by y = x 2 and y = x 3 
between x = and x = 1. 
The cross section is a "washer", and the area as a function of x is 

A(x) = Tr(r out (x) 2 - r m (x) 2 ) = 7r(a; 4 - x 6 ). 
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Figure 2.8: Disc method for computing a volume of revolution. 




-0,25 - 
-0.5 - 
-0.75 - 



-1 - 



Figure 2.9: Plots of y = x 2 and y = x 3 . 



The volume is thus 

J A(x)dx = J Q.t 5 - ^x 7 ^ dx = 7r 

Practice 2.2.2. Find the volumes swept out when 

(a) the region between f(x) = x and the x-axis, for < x < 2, is revolved 
about the x-axis, and 

(b) the region between f{x) = x and the line y = 2 X , for < x < 2, is revolved 
about the line y = 2. 

Example 2.2.4. One of the most important examples of a volume is the volume 
V of a sphere of radius r. Let's find it! We'll just compute the volume of a half 



= 7T. 

35 
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and multiply by 2. The cross sectional area is 
A(x) = nr(x)' 2 = ; 

Then 

1 r 

-V = / 7r(r 2 — x 2 )dx = n 

2 Jo 

Thus V = (4/3)7rr 3 . 

Example 2.2.5. Find volume of intersection of two spheres of radius r, where 
the center of each sphere lies on the edge of the other sphere. 



(Vr 2 -x 2 ) 2 = 7r(r 2 -x 2 ). 



2 3 

r x x 

3 



O 1 O 2 n 

= irr irr = — nr 

3 3 



Figure 2.10: Plot of two spheres. 
From the picture ( Figure [ 2. 10\ ) we see that the answer is 

2 f A(x), 

Jr/2 

where A(x) is exactly as in Example \2.2.4\ We have 

f r 5 
2 / 7r(r 2 - x 2 )dx = — 7rr 3 . 
Jr/2 1 ^ 12 

The previous ideas and techniques can also be used to find the volumes of 
solids with holes in them. If A{x) is the area of the face formed by a cut at x, 
then it is still true that the volume is V — J a A{x) dx. However, if the solid has 
holes, then some of the faces will also have holes and a formula for A(x) may be 
more complicated. Sometimes it is easier to work with two integrals and then 
subtract: (i) calculate the volume S of the solid without the hole, (ii) calculate 
the volume H of the hole, and (iii) subtract H from S. This is what was done 
in Example 12.2.31 

2.2.2 Shell method 

The disk method can be cumbersome if we want the volume when the region in 
the figure is revolved about the y-axis or some other vertical line. To revolve 
the region about the y-axis, the disk method requires that we represent the 
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original equation y = f(x) as a function of y: x = g(y). Sometimes that is easy: 
if y = 3x then x = y/3. But sometimes it is not easy at all: if y = x + e x , then 
we can not solve for y as an "elementary" function of x. The "shell" method 
lets us use the original equation y = f{x) to find the volume when the region 
is revolved about a vertical line. We partition the x-axis to cut the region into 
thin, almost rectangular "slices." When the thin "slice" at Xi is revolved about 
the y-axis (Figure [2.11( a)), the volume of the resulting "tube" (or cylindrical 
"shell") can be approximated by cutting the wall of the tube and laying it out 
flat (Figure [2.11( b)) to get a thin, solid rectangular box. 




radius = x. 
i 



Figure 2.11: Shell method for computing a volume of revolution. 

The volume of the tube is approximately the same as the volume of the solid 
box: 

Vol. tube = Vol. box = (length) x (height) x (thickness) 
= (27rradius) x ( height ) x (Ax») 
= 2nx i f(x l )Ax l . 

The volume swept out when the whole region is revolved is the sum of the 
volumes of these "tubes" , a Riemann sum. The limit of the Riemann sum is 



volume of rotation about the y-axis 



2nxf(x) dx. 



Theorem 2.2.3. (Volume of Revolution Using Shells) If region R is bounded 
between the functions f{x) > g(x) for < a < b (see Figure [2. 12) . then 



,6 

volume obtained when R is revolved about the y-axis = / 2nx(f(x)~g(x)) dx. 

J a 

Example 2.2.6. Find the volume when the region R in between x = 2, x = 4, 
y = x and y = x 2 is revolved about the y-axis. 
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Figure 2.12: Shell method for computing a volume of revolution. 



Solution: We can partition the interval [2,4] on the x~axis to get thin slices 
of R. When the slice at Xi is revolved around the y-axis, a tube is swept out, 
and the volume Vi of this i-th tube is 

Vi = (2ir ■ radius) x (height) x (thickness) 
= 2irxi(x 2 — Xi)Axi 
= 27r(xf - xf)Axi. 

The total approximate volume is the sum of the volumes of the tubes. As the 
partition gets finer and finer, we get 

V = 12 Vi ^ J 2^(x 3 -x 2 )da; = 27r(^-y)|4 = 27r^= 259.7... . 

2.2.3 Problems 

1 . For the solid in Figure 12.131 the face formed by a cut at x is a triangle 
with a base of 4 inches and a height of x 2 inches. Write and evaluate an 
integral for the volume of the solid for x between 1 and 2. 




Figure 2.13: Volume of a solid. 
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2. Find the volume of the squarebased pyramid in Figure 2.141 




Figure 2.14: Volume of a pyramid. 

3. Find the volume generated when the region between one arch of the sine 
curve (0 < x < ir) and the x-axis is revolved about (a) the a;-axis and (b) 
the line y = 1/2. 

4. Given that J x f(x) dx = 4 and f(x) 2 dx = 7. Represent the volumes 
of the solids (a), (b), (c) and (d) in Figure [2.151 as definite integrals and 
evaluate the integrals. 




Figure 2.15: Four volumes. 
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5. Figure |2.16j For < x < 3, each face is a circle with height (diameter) 
4 — x meters. 




Figure 2.16: Volume with circular cross-sections. 

6. Suppose A and B are solids so that every horizontal cut produces faces 
of A and B that have equal areas. What (if anything) can we conclude 
about the volumes of A and B? Justify your answer. 

7. Calculate the volume of a sphere of radius 2. 

8. Let < r < R be fixed. Revolve the circle x 2 + (y — R) 2 = r 2 about the 
x-axis. Compute the volume of this "donut" solid. 

9. (a) Find the area between f(x) = 1/x and the x-axis for 1 < x < 10, 

1 < x < 100, and 1 < x < A. What is the limit of the area for 
1 < 2: < A as A ^ oo? If ^4 = 1000000 and your think of this area 
as a long, flat wall, estimate the amount of paint (in square feet) you 
need to paint this surface. 

(b) Find the volume swept out when the region in part (a) is revolved 
about the x-axis for 1 < a; < 10, 1 < x < 100, and 1 < x < A. What 
is the limit of the volumes for 1 < x < A as A -> 00? If A = 1000000 
and your think of this volume as a room constructed by revolving 
the wall in (a) about an axis, estimate the amount of paint (in cubic 
feet) you need to completely fill the room. 

(c) Which is larger: the paint needed to paint the wall of the paint needed 
to completely fill the room? 

10. The region between y = 2x — x 2 and the xaxis for < x < 2. Sketch 
the region and calculate the volume swept out when the region is revolved 
about the yaxis. 

11. The region between y = \J\ — x 2 and the xaxis for < x < 1. Sketch 
the region and calculate the volume swept out when the region is revolved 
about the yaxis. 
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12. The region between y = ± ^ x2 and the xaxis for < x < 1. Sketch the 
region and calculate the volume swept out when the region is revolved 
about the yaxis. 

2.3 Average Values 

In this section we use Riemann sums to extend the familiar notion of an average, 
which provides yet another physical interpretation of integration. 

Recall: Suppose y±, . . . , y n are the amount of rain each day in you hometown 
so far this year. The average rainfall per day is 

_ yi h vy-a _\sr- 

2/avg / y% ■ 

i=i 

Definition 2.3.1 (Average Value of Function). Suppose f is a continuous func- 
tion on an interval [a, b]. The average value of f on [a,b] is 

1 ' 6 



/avg = 7 / f{x)dx. 

b-a J a 

Motivation: If we sample / at n points Xi, then 

1 " (h — \ " 1 " 

/av S - - E = ^ E m = ^ e /(-oa, 

since Ax = . This is a Riemann sum! 



1 If 

lim V f{xi)Ax = T r / f{x)dx. 

{b - a) n^oo ^ {b-a) J a 

This explains why we defined / avg as above. 

Example 2.3.1. What is the average value of sm(x) on the interval [0,7r]? 



TT - 



sm{x)dx 



1 



7T — 
1 



— cos(a;) 



7T 



(-l)-(-l) 



2 

7T 



Observation: If you multiply both sides by (6 — a) in Definition 12.3.11 you 
see that the average value times the length of the interval is the area, i.e., the 
average value gives you a rectangle with the same area as the area under your 
function. In particular, in Figure [2.171 the area between the x-axis and sin(x) is 
exactly the same as the area between the horizontal line of height 2/ir and the 
x-axis. 
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i 




075 



0.5 



25 



-1 



3 



-0.25 



-0.75 



-1 



Figure 2.17: What is the average value of sin(x)? 



Theorem 2.3.1 (Mean Value Theorem). Suppose f is a continuous function 
on [a, b]. Then there is a number c in [a, b] such that /(c) = /avg- 

This says that / assumes its average value. It is a used very often in under- 
standing why certain statements are true. Notice that in Example 1 2. 3. II it is just 
the assertion that the graphs of the function and the horizontal line y = / avg 
intersect. 

Proof: Let F(x) = f* f(t)dt. Then F'(x) = f(x). By the mean value 
theorem for derivatives, there is c G [a, b] such that /(c) = F'(c) = (F(b) — 
F(a))/(b — a). But by the fundamental theorem of calculus, 



In problems 1-4, use the values in Figure 12.18 to estimate the average values. 

1. Estimate the average value of / on the interval [0.5, 4.5]. 

2. Estimate the average value of / on the interval [0.5, 6.5]. 

3. Estimate the average value of / on the interval [1.5, 3.5]. 

4. Estimate the average value of / on the interval [3.5, 6.5]. 

5. Find the average value of sin(x), < x < n. 

6. Find the average value of x 2 , — 1 < x < 1. 




□ 



2.3.1 Problems 
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Figure 2.18: Table of values of the function f(x). 

2.4 Moments and centers of mass 

This section develops a method for finding the center of mass of a thin, flat 
shape the point at which the shape will balance without tilting. Centers of 
mass are important because in many applied situations an object behaves as 
though its entire mass is located at its center of mass. For example, if you are 
riding in a car with a high center of mass (such as an SUV) and you make a 
sudden sharp turn, you are more likely to tip over than if you are riding in a car 
with a high center of mass (such as a sports car) . As another example, the work 
done to pump the water in a tank to a higher point is the same as the work to 
move the center of mass of the water to the higher point (Figure |2.19j) , a much 
easier problem, if we know the mass and the center of mass of the water. Also, 
volumes and surface areas of solids of revolution can be easy to calculate, if we 
know the center of mass of the region being revolved. 




woik= (force) (di stance) 
Figure 2.19: Work depends on the center of mass. 

Before looking for the centers of mass of complicated regions, we consider 
point masses and systems of point masses, first in one dimension and then in 
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two dimensions. 

2.4.1 Point Masses 

First we discuss point masses along a line. 

Two people with different masses can position themselves on a seesaw so that 
the seesaw balances (Figure [2. 20 j) . The person on the right causes the seesaw to 
"want to turn" clockwise about the fulcrum, and the person on the left causes 
it to "want to turn" counterclockwise. If these two "tendencies" are equal, the 
seesaw will balance. A measure of this tendency to turn about the fulcrum is 
called the moment about the fulcrum of the system, and its magnitude is the 
mass multiplied by the distance from fulcrum. 



Q 

t* 

Figure 2.20: Balance on a see-saw depends on the center of mass. 

In general, the moment about the origin, Mo, produced by a mass m at a 
location x is mx, the product of the mass and the "signed distance" of the mass 
from the origin. For a system of masses mj , m 2 , . . . , m n at locations x\, x%, 
. . . , x n , respectively, 

M = total mass of the system 

and 



Mq = moment about the origin = x\m\ + 2:2^2 H + x n m n — Xiirii. 

i=i 

If the moment about the origin is positive then the system tends to rotate 
clockwise about the origin. If the moment about the origin is negative then the 
system tends to rotate counterclockwise about the origin. If the moment about 
the origin is zero, then the system does not tend to rotate in either direction 
about the origin; it balances on a fulcrum at the origin. The moment about the 
point p, M v , produced by a mass m at the location x is the signed distance of x 
from p times the mass m: (x — p) ■ m. The moment about the point p produced 
by masses mi, 7712, . . . , m n at locations x\, X2, ■ ■ ■ , x n , respectively, is 




n 

»=l 



73 



2.4. MOMENTS AND CENTERS OF MASS 



Mp = moment about p = (xi~p)mi+{x2—p)m2+- ■ ■+{x n —p)m n = p) 



i=l 



The point at which the system balances is called the center of mass of the 
system and is written x (pronounced "a; bar"). Since the system balances at 
x, the moment about p = x must be zero. Using this fact and properties of 
summation, we can find a formula for x: 



2 — 1 i— 1 2 — 1 

so 

En 
i=l x i m i 

E i= i m { 

This is summarized as follows. 

Theorem 2.4.1. The center of mass of a system of masses mi, mi, ■ ■ ■ , m n 

at locations x\, xi, . . . , x n , is given by 

moment about the oriqin 

x = — y — = M M. 

total mass 

Now wc discuss point masses in the plane. 

The ideas of moments and centers of mass extend nicely from one dimension 
to a system of masses located at points in the plane. 
For a system of masses m., located at the points (xt,yi), 



M = total mass of particles = Y^h=i m ii 
My = moment about the y-axis = X)"=i m i x ii 
M x = moment about the x-axis = X)"=i m i2/i- 

Theorem 2.4.2. The center of mass of a system of masses mi, mi, 
at locations (x^yi), (x 2 ,y 2 ), ■ • (x n ,y n ), is given by (x,y), where 



* = = M v/ M 



moment about the y-axis 

total mass 
moment about the x-axis 

M x /M. 



total mass 

Example 2.4.1. Consider a regular hexagon in the plane centered at the origin^] 
Suppose that all the vertices of the hexagon have equal mass 1. The center of 
mass of this hexagon is the same as the average value of the vertices. Therefore, 
by construction, (x,y) = (0,0). 



x To draw this, simply draw a circle and slice it up fairly into 6 equal "pie pieces". The 
points on the "crust" where your slices start are the vertices of the hexagon. 
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2.4.2 Center of mass of a region in the plane 

When we move from discrete point masses to whole, continuous regions in the 
plane, we move from finite sums and arithmetic to limits of Ricmann sums, 
definite integrals, and calculus. The following material extends the ideas and 
calculations from point masses to uniformly thin, flat plates that have a constant 
density given as mass per area. The center of mass of one of these plates is the 
point (x, y) at which the plate balances without tilting. It turns out that the 
center of mass (x, y) of such a plate depends only on the region of the plane 
covered by the plate and not on its density. In this situation, the point (x,y) is 
also called the centroid of the region. In the following discussion, you should 
notice that each finite sum that appeared in the discussion of point masses has 
a counterpart for these thin plates in terms of integrals. 

The rectangle is the basic shape used to extend the point mass ideas to regions. 
The total mass of a rectangular plate is the area of the plate multiplied by the 
density constant: mass M = area x density. We assume that the center of 
mass of a thin, rectangular plate is located half way up and half way across the 
rectangle, at the point where the diagonals of the rectangle cross. Then the 
moments of the rectangle can be found by treating the rectangle as a point with 
mass M located at the center of mass of the rectangle. 

To find the moments and center of mass of a plate made up of several rectangu- 
lar regions, just treat each of the rectangular pieces as a point mass concentrated 
at its center of mass. Then the plate is treated as a system of discrete point 
masses. 



Example 2.4.2. The plate in Figure \2. 21 can be divided into two rectangular 
plates, one with mass 24 g and center of mass (1,4) , and one with mass 12 g 
and center of mass (3, 3). 




Figure 2.21: Centroid of two rectangles. 



The total mass of the pair is M = 36 g, and the moments about the axes are 
M x = (24 g) (4 cm) + (12 g) (3 cm) = 132 gem, and M y = (24 g)(l cm) + 
(12 g) (3 cm) — 60 gem. 

Thenx = M y /M = (60 gem)/ (36 g) = 5/3 cm andy = M x /M = (132 gem)/ (36 g) 
11/3 cm so the center of mass of the plate is at (x,y) = (5/3,11/3). 
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To find the center of mass of a thin plate, we will "slice" the plate into narrow 
rectangular plates and treat the collection of rectangular plates as a system of 
point masses located at the centers of mass of the rectangles. The total mass 
and moments about the axes for the system of point masses will be Riemann 
sums. Then, by taking limits as the widths of the rectangles approach 0, we 
will obtain exact values for the mass and moments as definite integrals 

2.4.3 x For A Region 

Suppose f(x) > g{x) on [a, b] and R is a plate on the region between the graphs 
of / and g for a < x < b (Figure [2. 22) ) . 




Figure 2.22: Centroid of a region. 



If the interval [a, b] is partitioned into subintervals [xj_i,a;j] and the point Cj is 
the midpoint of each subinterval, then the slice between vertical cuts at Xi—\ 
and Xi is approximately rectangular and has mass approximately equal to 

(area) x (density) = (height) x (width) x (density) 
= (/(cj) - g(c i ))(xi_ 1 x i )k 

where fc denotes the density. The mass of the whole plate is approximately 

M = ^2(f( c i) - 9^))^ -> k / (f(x)g(x)) dx 
= k ■ (area of the region between / and g). 
The moment about the y-axis of each rectangular piece is 

My = (distance from y-axis to center of mass of piece) x (mass) 
= Ci(f(ci) - g(a))(Axi)k, 
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so 



M v = J2 c *(/( c ») - g(c t ))(Ax t )k x(f(x) - g(x)) dx. 

i Ja 

The x-coordinate of the center of mass of the plate is 

w= My = J*x(f(x)-g(x))dx 
M ' f b a m-g(x)dx ' 

since the common factor of fc on top and bottom cancel. 

Practice 2.4.1. Find the x-coordinate of the center of mass of the region be- 
tween f(x) — x 2 and the x-axis for < x < 2. (In this case, g(x) = 0.) 

2.4.4 y For a Region 

Again, suppose f(x) > g(x) on [a, b] and R is a plate on the region between the 
graphs of / and g for a < x < b (Figure [2.22 j ). 



i i i j 



c . 
1 



Figure 2.23: Finding the y-coordinate of the centroid of a region. 

To find y, the y-coordinate of the center of mass of the plate R, we need to find 
M x , the moment of the plate about the cc-axis. When R is partitioned vertically 
(Figure [2.23j h the moment of each (very narrow) strip about the x-axis, M x , is 



(signed distance from x-axis to the center of mass of strip) x (mass of strip). 

Since each thin strip is approximately rectangular, the y-coordinate of the center 
of mass of each strip is approximately half way up the strip: y i = (/(cj) + 
g{d))/2. Then 
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M x for the strip = (signed distance from the x-axis 

to the center of mass of the strip) x (mass of strip) 
= (signed distance from xaxis) x ( height of strip) 
x (width of strip) x (density constant) 
f{ci) t 9(ci) (f(ci)-9(ci))(A Xi )k. 

The moment about the x~axis of each rectangular piece is 



(distance from the x-axis to the center of mass of the piece) x (mass) 
= /(cO+ g (cO (/(c . ) _ g(c . ))(Aa; . )fc) 



so 



M x = /(Ci T (Cij (/fa) - g( Cl ))(A Xl )k 

- kj b a iM+iM {f(x) _ g(x)) dx = | f a f{x f _ g[x f dx _ 

Practice 2.4.2. Show that the centroid of a triangular region with vertices 
(0,0), (0,h) and (6,0) is (x,y) = (6/3, h/3). 

Example 2.4.3. Find the centroid of the region bounded between the graphs of 
y = x and y = x 2 , for < x < 1. 

Solution: M = k f (% — x 2 ) dx = k/6, M y = k f Q x(x — x 2 ) dx = k/12 and 
M x = | $l(x 2 - x 4 ) dx = k/15. Then x = M y /M = 1/2 and y = M x /M = 2/5 



2.4.5 Theorems of Pappus 

When location of the center of mass of an object is known, the theorems of 
Pappus make some volume and surface area calculations very easy. 

Volume of Revolution: If a plane region with area A and centroid (x, y) is re- 
volved around a line in the plane which does not go through the region (touching 
the boundary is alright), then the volume swept out by one revolution is the 
area of the region times the distance traveled by the centroid (Figure T2.24j ): 

Theorem 2.4.3. (Pappas' theorem for volume) 

Volume about line L = A ■ 2tt distance of (x, y) from the line L , 
so in particular, 

Volume about x-axis = A ■ 2ny, 

and 

Volume about y-axis = A ■ 2ttx. 
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area sweeps out a volume 




Figure 2.24: Pappas' theorem for a volume of revolution. 

Surface Area of Revolution If a plane region with perimeter P and centroid of 
the edge (x, y) is revolved around a line in the plane which does not go through 
the region (touching the boundary is alright), then the surface area swept out 
by one revolution is the perimeter of the region times the distance traveled by 
the centroid (Figure |2.25j ) : 

perimeter sweeps out an area (surface area) 




Figure 2.25: Pappas' theorem for a surface area of revolution. 

Theorem 2.4.4. (Pappas' theorem for surface area) 

Surface area about line L = P ■ 2tt distance of (x, y) from the line L , 
so in particular, 

Surface area about x-axis = P ■ 2?:y, 
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and 

Surface area about y-axis = P ■ 2ttx. 

Example 2.4.4. The center of a square region with 2 foot sides is at the point 
(3,4). Use the Theorems of Pappus to find the volume and surface area swept 
out when the square is rotated (a) about the x-axis, (b) about the y-axis, and 
(c) about the horizontal line y = 6. 

Solution: (a) Volume about x-axis = A-2ny = 32tt, Surface area about x-axis 
= P ■ 2ny = 64tt. 

(b) Volume about y-axis = A-2nx = 247r, Surface area about y-axis = P-2ixx = 
48tt. 

(c) Volume about line y = 6 = A -2t: (distance o/(3, 4) to the line y = 6) =l6ir, 
Surface area about line y = 6 = P -2tt (distance of (3,4) to the line y = 6) =327r. 

2.5 Arc lengths 

This section introduces another geometric applications of integration: finding 
the length of a curve, i.e., the total distance you travel if you are moving along 
a curve. The general strategy is the same as before: partition the problem 
into small pieces, approximate the solution on each small piece, add the small 
solutions together in the form of a Riemann sum, and finally, take the limit of 
the Riemann sum to get a definite integral. 

2.5.1 2-D Arc length 

Suppose C is a curve, and we pick some points (xi,j/i) along C and connect 
the points with straight line segments. Then the sum of the lengths of the line 
segments will approximate the length of C . We can think of this as pinning 
a string to the curve at the selected points, and then measuring the length 
of the string as an approximation of the length of the curve. Of course, if 
we only pick a few points, then the total length approximation will probably 
be rather poor, so eventually we want lots of points (xi,yi), close together all 
along C. Suppose the points are labeled so (xo;2/o) 1S one endpoint of C and 
{%m Un) is the other endpoint and that the subscripts increase as we move along 
C. Then the distance between the successive points (xj_i, yi—i) and (xi,yi) is 
•*/ (Axi) 2 + (Aj/i) 2 , and the total length of the line segments is simply the sum 
of the successive lengths. This is an important approximation of the length of 
C, and all of the integral representations for the length of C come from it. The 
length of the curve C is approximately 

J2 V(^i) 2 + (Ay,) 2 = E + (A yi /Ax 4 ) 2 Ax 4 . 

i i 

This is a Riemann sum. We could have factored out a Ayi instead: the arc 
length of C is approximately 
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Theorem 2.5.1. • The length of the curve C described by the graph of the 
function y — f{x), a < x < b, is given by 

b 

y/l + f'(x)*dx. 

• The length of the curve C described by the graph of the function x = g{y), 
O < U 5= b, is given by 

,6 

/ V^+dWdy. 

J a 

Example 2.5.1. Use the points (0,0), (1,1), and (3,9) to approximate the 
length of y = x 2 , for < x < 3. 

Solution: The lengths of the pieces are y(T— 0) 2 + (1 — 0) 2 = -\/2 and 
-\/(3 — l) 2 + (9 — l) 2 = \/68, so i/ie totoZ length is approximately y/2 + %/68 = 
9.66... . 

W^e can wse Sage to compute some more Riemann sum approximations and 
also this arc length exactly. 



sage: fl = lambda x: sqrt (l+4*x~2) 
sage: f = PiecewiseC [[(0,3) ,f 1]] ) 

sage: n = 10; RR(f . riemann_sum_integral_approximation(n) ) 
8.99946939777166 

sage: n = 50; RR(f . riemann_sum_integral_approximation(n) ) 
9.59519771936512 

sage: n = 100; RR(f .riemann_sum_integral_approximation(ii) ) 
9.67099527976211 

sage: n = 200; RR(f .riemann_sum_integral_approximation(n) ) 
9.70900502940468 

sage: integral(sqrt(l+(2*x) ~2) ,x,0,3) 

(arcsinh(6) + 6*sqrt (37) ) /4 

sage: RR (integral (sqrt (l+(2*x) ~2) ,x,0,3)) 

9.74708875860856 



In other words, 



sinh -1 6 + 6\/37 



v/l + (2.x) 2 dx = ^— — — = 9.74... . 

Note that if we reflect this curve about the y = x line then the resulting part of 
the curve must have the same arc length. Agreed? Do you see that the reflected 
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curve is y = *Jx, < x < 9? In this case the reflected points are (0,0),. (1, 1), 
and (9,3). The lengths of the reflected pieces are yJY— 0) 2 + (1 — 0) 2 = \[2 
and \/(9 — l) 2 + (3 — l) 2 = -s/68, so the total length is (still) approximately 
\/2 + \/68 = 9.66... . The integral describing the arc length is 



1 + (^x- 1 / 2 ) 2 dx= [ Jl + -x- 1 dx. 
o v 2 ,y V 4 

This is a harder integral to computed but we can use Sage to compute this 
reflected arc length fairly accurately. 



sage: fl = lambda x: sqrt(l+l/(4*x)) 

sage: mimerical_integral(f 1 , 0, 9, max_points=100) 

(9 . 7470886680795221 , 7 . 9546440984616276e-06) 



The output is a pair, the first coordinate is the approximate numerical value of 
the integral and the second is an upper bound for the error term. This is in 
agreement with the above answer. 

Parametric equations: When the curve C is described by pairs (x, y), where 
x and y are functions of t, x = x(t) and y = y(t), for a < t < f3, we can factor 
(Ati) from inside the radical and simplify: 

J2 V(^) 2 + (A^) 2 = J2 Vi^/^ti) 2 + (AyJAuyAU. 

i i 

This is a Riemann sum. Taking limits, we get the following formula. 

Theorem 2.5.2. The length of the curve C described by the graph of the para- 
metric equations x = x(t) and y = y(t), for a < t < (3, is given by 

f ^x'(t) 2 + y'(t) 2 dt. 

J OL 

Example 2.5.2. Represent the length of each curve as a definite integral. 

(a) The length of y = e x between (0, 1) and (1, e). 

(b) The length of the parametric curve x{t) = cos(t) and y(t) = sm(t) for 
< t < 2tt . 

Solution: (a) J Q v/T-j- e 2x dx. This looks complicated (and it is) but amazingly 
enough Sage has no problem with it: 



sage: fl = lambda x: sqrt (l+exp(2*x) ) 
sage: integral (fl (x) , x, 0, 1) 

-arcsinh(e~-l) + arcsinh(l) + sqrt(e~2 + 1) - sqrt (2) 

2 Of course, there is no need to, since we already know its value! 
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sage: RR ( integral (f l(x) , x, 0, 1)) 
2.00349711162735 



(b) V(-sin(*)) 2 + (cos(i)) 2 d* = J^Vldt = 2tt. 
2.5.2 3-D Arc length 

The parametric equation form of arc length extends very nicely to 3 dimensions. 
If a curve C in 3-dimensions (Figure [2.26|) is given parametrically by x = x(t), 
y = y(t), and z = z(t) for a < t < b, then the distance between the successive 
points and (xuy%,Zi) is ^Aa; 2 + Ayf + Azf. 




Figure 2.26: The arc length or a space curve. 

We can, as before, factor (Ati) 2 from each term under the radical, sum the 
pieces to get a Ricmann sum, and take a limit of the Riemann sum to get a 
definite integral representing the length of the curve C. 

yfch + Azj? + Azf = y/iAxi/Atif + (Aj/i/Ati) 2 + (Az t /At t ^At t 

i 

b 

V^W + y'(t) 2 + z'{t) 2 dt. 

Theorem 2.5.3. If a curve C in 3- dimensions is given parametrically by x = 
x(t), y = y(t), and z = z(t) for a < t < b, then the arc length is 

b 

\/x'{t) 2 + y'(t) 2 + z'{t) 2 dt. 
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Example 2.5.3. Find the arc length of the helix x = cos(t), y = sin(i), z = t 
forO <t< 4?r. 
Solution: We have 



yjx'{t) 2 + y'{ty + z'(t) 2 alt = J Q 47r y/sin(t) 2 + cos(t) 2 + 1 alt 
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Polar coordinates and 
trigonometric integrals 

The rectangular coordinate system is immensely useful, but it is not the only 
way to assign an address to a point in the plane and sometimes it is not the 
most useful. In applications to physical problems where there is some "cylin- 
drical symmetry" , such as a vibrating drum or water moving along a pipe, the 
most natural coordinates are often polar coordinates rather than rectangular 
coordinates. 

In many experimental situations, our location is fixed and we, using sonar or 
radar, take readings in different directions (Figure HO) ; this information can be 
graphed using rectangular coordinates (e.g., with the angle on the horizontal 
axis and the measurement on the vertical axis). 




Taking Measurements Rectangular Coordinate Graph of Data 



Figure 3.1: Sonar and radar use polar coordinates. 
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Sometimes, however, it is more useful to plot the information in a way similar 
to the way in which it was collected, as magnitudes along radial lines (Figure 
13. 2) . This system is called the Polar Coordinate System. 



ss:' 




270" 

Polar C oordi nate Graph of Data 



Figure 3.2: Polar coordinates. 

Example 3.0.4. SOS! You've just received a distress signal from a ship located 
at A on your radar screen (Fiqure \3.3\ ). Describe its location to your captain so 
your vessel can speed to the rescue. 




Figure 3.3: Polar coordinate figure for Example 13.0.41 

Solution: You could convert the relative location of the other ship to rectan- 
gular coordinates and then tell your captain to go due east for 7.5 miles and 
north for 13 miles, but that certainly is not the quickest way to reach the other 
ship. It is better to tell the captain to sail for 15 miles in the direction of 60° . 
// the distressed ship was at B on the radar screen, your vessel should sail for 
10 miles in the direction 150°. (Real radar screens have 0° at the top of the 
screen, but the convention in mathematics is to put 0° in the direction of the 
positive x-axis and to measure positive angles counterclockwise from there. And 
of course a real sailor speaks of "bearing" and "range" instead of direction and 
magnitude.) 
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Practice 3.0.1. Describe the locations of the ships at C and D in Figure \373\ 
by giving a distance and a direction to those ships from your current position at 
the center of the radar screen. 

Points in Polar Coordinates: To construct a polar coordinate system we need 
a starting point (called the origin or pole) for the magnitude measurements and 
a starting direction (called the polar axis) for the angle measurements. A polar 
coordinate pair for a point P in the plane is an ordered pair (r, 9), where r is 
the directed distance along a radial line from O to P, and 9 is the angle formed 
by the polar axis and the segment OP. The angle 8 is positive when the angle 
of the radial line OP is measured counterclockwise from the polar axis, and 9 
is negative when measured clockwise. 

Degree or Radian Measure for 9?: Either degree or radian measure can be 
used for the angle in the polar coordinate system, but when we differentiate 
and integrate trigonometric functions of 9 we will always want all of the angles 
to be given in radian measure. From now on, we will primarily use radian 
measure. You should assume that all angles are given in radian measure unless 
the units "°" ("degrees") are shown. 

In the rectangular coordinate system, the derivative dy/dx measured both the 
rate of change of y with respect to x and the slope of the tangent line. In the 
polar coordinate system dr/d9 measures the rate of change of r with respect to 
9. The sign of dr/d9 tells us whether r is increasing or decreasing as 9 increases. 



r(6 > 

As & increases 

so does r: 

d B 



Figure 3.4: r changing as a function of 9. 



3.1 Polar Coordinates 

Rectangular coordinates allow us to describe a point (x, y) in the plane in a 
different way, namely 

(x,y) <-> (r,8), 

where r is any real number and 9 is an angle. 

Polar coordinates arc extremely useful, especially when thinking about com- 
plex numbers. Note, however, that the (r, 9) representation of a point is very 
non-unique. 
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First, 9 is not determined by the point. You could add 2n to it and get the 
same point: 

(^)=( 2 '?H'i+— )-(>■=?) 

Also that r can be negative introduces further non-uniqueness: 




Think about this as follows: facing in the direction 37r/2 and backing up 1 meter 
gets you to the same point as looking in the direction 7r/2 and walking forward 
1 meter. 

We can convert back and forth between Cartesian and polar coordinates using 
that 



x = r cos(0) (3-1) 

y = ran(e), (3.2) 

and in the other direction 

r 2 =x 2 +y 2 (3.3) 

tan(0) = - (3.4) 

x 

(Thus r = ±\/x 2 + y 2 and 8 — tan -1 (?/ /x).) 




x = r-cos{ 8 ) 
y = r- sin( 6 > 




x 2 + v 2 =r 



VAil ( > '— 



Figure 3.5: Rectangular to polar coordinate conversion. 



Example 3.1.1. Sketch r = sin(#) ; which is a circle sitting on top the x axis. 

We plug in points for one period of the function we are graphing — in this case 
[0,2tt]: 
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Figure 3.6: Plot of (a) r = 1, < 6 < tt/6, and (b) r = sin(6»), < 9 < 2n. 






sin(O) = 


tt/6 


sin(vr/6) = 1/2 


tt/4 


sin(7r/4) = & 


tt/2 


sin(7r/2) = 1 


3tt/4 


sin(37r/4) = ^ 


7T 


sin(7r) = 


7T + TT/6 


sin(7r + 7r/6) = -1/2 



Notice it is nice to allow r to be negative, so we don't have to restrict the input. 
BUT it is really painful to draw this graph by hand. 

To more accurately draw the graph, let's try converting the equation to one 
involving polar coordinates. This is easier if we multiply both sides by r: 

r 2 = r sin(9). 

Note that the new equation has the extra solution (r = 0, 6 = anything), so we 
have to be careful not to include this point. Now convert to Cartesian coordinates 
using [3J\ ) to obtain (373\ ): 

x 2 +y 2 =y. (3.5) 

The graph of (3. 5\) is the same as that of r = sm(8). To confirm this we complete 
the square: 

2 i 2 

x + y = y 
x 2 +y 2 -y = 
x 2 + (y-l/2) 2 = l/4 

Thus the graph of (3. 5\ ) is a circle of radius 1/2 centered at (0, 1/2). 

Actually any polar graph of the form r = a sm(9) + b cos(6*) is a circle (exercise 
for the interested reader). 
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3.2 Areas in Polar Coordinates 

The previous section introduced the polar coordinate system and discussed how 
to plot points, how to create graphs of functions (from data, a rectangular 
graph, or a formula), and how to convert back and forth between the polar 
and rectangular coordinate systems. This section examines calculus in polar 
coordinates: rates of changes, slopes of tangent lines, areas, and lengths of 
curves. The results we obtain may look different, but they all follow from the 
approaches used in the rectangular coordinate system. 

We know how to compute the area of a sector, i.e., piece of a circle with angle 
9. [[draw picture]]. This is the basic polar region. The area is 

A = (fraction of the circle) • (area of circle) = (^^j ' 717,2 = 2 r2 ^' 

We now imitate what we did before with Riemann sums. We chop up, ap- 
proximate, and take a limit. Break the interval of angles from a to b into n 
subintervals. Choose 9* in each interval. The area of each slice is approxi- 
mately (l/2)/(6»*) 2 6f. Thus 



n 

A = Area of the shaded region - ^ - f(9*) 2 A{9). 
Taking the limit, we see that 



2" 

i=i 



A = l^E^m'fAW = - ■ \ f(9) 2 d9. (3.6) 

Amazing! By understanding the definition of Riemann sum, we've derived a 
formula for areas swept out by a polar graph. But does it work in practice? 

Example 3.2.1. Find the area enclosed by one leaf of the four-leaved rose 
r = cos(26»). 




Figure 3.7: Graph of y = cos(2x) and r = 008(26*). 
This was plotted in Sage using these commands: 
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sage: PI = polar_plot (lambda x:cos(2*x), 0, 2*pi, rgbcolor=(0 ,0 , 1) ) 

sage: P2 = plot (lambda x:cos(2*x), 0, 2*pi, rgbcolor=( 1,0,0) ,linestyle=" : ") 

sage: show(Pl+P2) 



To find the area using the methods we know so far, we would need to find a 
function y = f(x) that gives the "height" of the leaf. 
Multiplying both sides of the equation r = cos (2$) by r yields 

r 2 = rcos(26>) = r(cos 2 9 - sin 2 9) = -((rcosfl) 2 - (rsmO) 2 ). 

r 

Because r 2 = x 2 + y 2 and x — r cos((9) and y = rsm(9), we have 

2,2 1 / 2 2\ 

x +y = i 9 , A x -v )• 
v x + v 

Solving for y is a crazy mess, and then integrating? It seems impossible! 

But it isn't... if we remember the basic idea of integral calculus: integral equals 
area. 

We need the boundaries of integration to determine the area. Start at 9 = 
— 7r/4 and go to 9 — 7r/4. As a check, note that cos((— 7r/4) • 2) = = cos((7r/4) • 
2). We evaluate 



7r/4 /■t'74 

2ja / „„„tr,n\'2 



cos(2#) d9 = cos(2#) d9 (even function) 

tt/4 JQ 

tt/4 

(1 + cos(40))d0 



1 

2 

1 

2 

7T 







\ ■ sin(40) 



tt/4 




We used that 

cos 2 {x) = (1 + cos(2a;))/2 and sin 2 (x) = (1 - sin(2a;))/2, (3.7) 
which follow from 

cos(2a;) = cos 2 (a;) - sin 2 (a;) = 2 cos 2 (a;) -1 = 1- 2 sin 2 (a;). 
Therefore, by (3.6\) , the area is A = ^. 

Example 3.2.2. Find area of region inside the curve r = 3cos(0) and outside 
the cardiod curve r = 1 + cos(#). 
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1.5 



-1 



-1.5 



Figure 3.8: Graph of r = 3cos(a;) and r = 1 + 008(6*). 



Figure \ 3.8\ was plotted in Sage using these commands: 



sage: PI = polar_plot (lambda x:3*cos(x), 0, 2*pi, rgbcolor= (0,0,1)) 

sage: P2 = polar_plot (lambda x:l+cos(x), 0, 2*pi, rgbcolor=( 1,0,0) ,linestyle=" : ") 

sage: show(Pl+P2) 



Solution: This is the same as before. It's the difference of two areas. Figure out 
the limits, which are where the curves intersect, i.e., the 9 such that 



3cos(6>) = 1 + cos(0). 



Solving, 2cos(6 ) ) = 1, so cos(9) = 1/2, hence 9 = 7r/3 and 9 



7r/3. Thus the 



area is 




since the integrand is an even function. Now, expand this out algebraically and 
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integrate term-by-term: 

tt/3 



r/3 



(3cos(0)) 2 - (1 + cos(6)yd6 = / (8cos 2 (<9) - 2cos(6>) - l)dO 

Jo 

tt/3 / 1 

I 8 ■ - (1 + cos(20)) - 2cos(0) - 1 ) d0 

tt/3 



/ 3 + 4cos(26»)-2cos(6i)d6i 



3^ + 2sin(20) -2sin(6>) 



tt/3 



Practice 3.2.1. T7ie area of the shaded region in Figure \37Uj 




3-H/2 

Figure 3.9: Graph of r 



3.3 Complex Numbers 

A complex number is an expression of the form a + bi 7 where a and b are real 
numbers, and i' 2 = — 1. We add and multiply complex numbers as follows: 



(a + bi) + (c + di) = (a + c) + (b + d)i 
(a + bi) ■ (c + di) = (ac — bd) + (ad + bc)i 

The complex conjugate of a complex number is 

a + bi = a — bi. 
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Note that 



(a + bi)(a + bi) = a 2 + b 2 



is a real number (has no complex part). 
If c + di ^ 0, then 



a + bi (a + bi)(c—di) 1 



((ac + bd) + (be - ad)i). 



c + di c? + d 2 c 2 + d 2 



Example 3.3.1. (l-2i)(8-3i) = 2— 19z and l/(l+i) = (1-i)/ 2 = l/2-(l/2)i. 

Complex numbers are incredibly useful in providing better ways to under- 
stand ideas in calculus, and more generally in many applications (e.g., electrical 
engineering, quantum mechanics, fractals, etc.). For example, 

• Every polynomial f(x) factors as a product of linear factors (x — a), if we 
allow the a's in the factorization to be complex numbers. For example, 



This will provide an easier to use variant of the "partial fractions" inte- 
gration technique, which we will see later. 

• Complex numbers are in correspondence with points in the plane via 
(x, y) <-» x + iy. Via this correspondence we obtain a way to add and 
multiply points in the plane. 

• Similarly, points in polar coordinates correspond to complex numbers: 



• Complex numbers provide a very nice way to remember and understand 
trig identities. 

3.3.1 Polar Form 

The polar form of a complex number x + iy is r(cos(#) + isin(0)) where (r, 6) 
are any choice of polar coordinates that represent the point (x, y) in rectangular 
coordinates. Recall that you can find the polar form of a point using that 



NOTE: Historically, the "existence" of complex numbers wasn't generally ac- 
cepted until people got used to a geometric interpretation of them. 

Example 3.3.2. Find the polar form of 1 + i. 



f(x) = x 2 + 1 = (x - i)(x + i). 



(r, 6) <-> r(cos(0) + i sin(0)). 



\/ x 2 + y 2 and 9 = tan 1 (y/x). 



Solution. We have r = \/2, 



so 
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Example 3.3.3. Find the polar form of ' y/3 — i. 
Solution. We have r = \/3 + 1 = 2, so 

V3-i = 2 [^Y+ i ^f] = 2(cos(-tt/6) +«sin(-7r/6)) 









jf u 

/ 1 

' 0.5 








-1.5 -1 -0.5 

\ -1 
\^ 1.5 


\M 1 1.5 


2.5 







Figure 3.10: Plot of \fi — i, as a vector. 
This was plotted in Sage using these commands: 

sage: PI = circle ( (0 , 0) , 2) 

sage: P2 = arrow((0,0) , (sqrt(3) ,-1)) 

sage: P3 = text("$\sqrt{3}-i$" , (2,-1)) 

sage: P4 = text ("$-i$" , (0 . 2 , -1) ) 

sage: P5 = text ("$i$" , (0 . 2, 1) ) 

sage: show(Pl+P2+P3+P4+P5) 

Finding the polar form of a complex number is exactly the same problem as 
finding polar coordinates of a point in rectangular coordinates. The only hard 
part is figuring out what 9 is. 

If we write complex numbers in rectangular form, their sum is easy to compute: 

(a + bi) + (c + di) = (a + c) + (b + d)i 

The beauty of polar coordinates is that if we write two complex numbers in 
polar form, then their product is very easy to compute: 

r 1 (cos(^i)+isin(^i))-r 2 (cos(6' 2 )+ism(6'2)) = (r 1 r 2 )(cos{9 1 +9 2 )+ism(9 1 +9 2 )). 

The magnitudes multiply and the angles add. The above formula is true because 
of the double angle identities for sin and cos: 
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(cos(0i) +isin(0i)) • (cos(0 2 ) + isin(0 2 )) 

= (cos (0i ) cos (0 2 ) - sin(0i) sin(0 2 )) + i(sin(0i) cos(0 2 ) + cos(0i) sin(0 2 )) 
= cos(0i + 2 ) + i sin(0x + 2 )). 

For example, the power of a singular complex number in polar form is easy to 
compute; just power the r and multiply the angle. 

Theorem 3.3.1 (De Moivrc's). For any integer n we have 

(r(cos(0) + isin(0))) n = r n (cos(n0) + isin(n0)). 

Example 3.3.4. Compute (1 + i) 2006 . 
Solution. PFe /lave 

(1 + i) 2006 = (V2 (cos(7r/4) + » sin(7r/4))) 2006 

/ — 2006 . , . . , . . „ 

= V2 (cos(2006tt/4) + i sin(20067r/4))) 
= 2 1003 (cos(37r/2) + isin(37r/2))) 
= -2 1003 i 

To get cos(20067r/4) = cos(37r/2) we use that 2006/4 = 501.5, so by periodicity 
of cosine, we have 

cos(2006tt/4) = cos((501.5)tt - 250(2tt)) = cos(1.5tt) = cos(3tt/2). 

Here's a quick summary of what we've just learned: Given a point (x, y) in the 
plane, we can also view it as x + iy or in polar form as r(cos(0) + zsin(0)). Polar 
form is great since it's good for multiplication, powering, and for extracting 
roots: 

ri(cos(0i) + isin(0i))r 2 (cos(0 2 )+isin(0 2 )) = (nr 2 )(cos(0i +0 2 ) +isin(0i +0 2 )). 

(If you divide, you subtract the angle.) The point is that the polar form works 
better with multiplication than the rectangular form. For any integer n, we have 

(r(cos(0) + * sin(0)))" = r"(cos(?i0) + * sin(n0)). 

Since we know how to raise a complex number in polar form to the nth power, 
we can find all numbers with a given power, hence find the nth roots of a complex 
number. 

Proposition 3.3.1 (nth roots). A complex number z = r(cos(0) + i sin(0)) has 
n distinct nth roots: 

rl (^cos^— — ]+*sm(— — jj, 
for k = 0, 1, . . . , n — 1. Here r x l n is the real positive n-th root of r. 
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As a double-check, note that by De Moivre, each number listed in the propo- 
sition has nth power equal to z. 

An application of De Moivre is to computing sin(n0) and cos(n0) in terms of 
sin(0) and cos((9). For example, 

cos(30) + isin(30) = (cos(0) + i sin(0)) 3 

= (cos(0) 3 - 3 cos (0) sin(0) 2 ) + i(3 cos(0) 2 sin(0) - sin(0) 3 ) 

Equate real and imaginary parts to get formulas for cos(30) and sin(30). In the 
next section we will discuss going in the other direction, i.e., writing powers of 
sin and cos in terms of sin and cosine. 

Example 3.3.5. Find the cube roots of 2. 
Solution. Write 2 in polar form as 

2 = 2(cos(0) -Hsin(O)). 

Then the three cube roots of 2 are 

2 1/3 (cos(27rfc/3) +isin(27rfc/3)), 

fork = 0,1,2. I.e., 

2 i/3 j 2 1 / 3 (-l/2 + z^/2), 2^ 3 (-l/2-iVS/2). 

3.4 Complex Exponentials and Trigonometric Iden- 
tities 

Recall that 

n (cos(0i) + isin(0 1 ))r 2 (cos(0 2 )-(-isin(0 2 )) = (rir 2 )(cos(0 1 + 2 )+isin(0 1 +02)). 

(3.8) 

The angles add. You've seen something similar before: 

a b a+b 

e e = a . 

This connection between exponentiation and ( j3.8|) gives us an idea! 
If z = x + iy is a complex number, define 

e z = e x (cos(y) + ism(y)). 

We have just written polar coordinates in another form. It's a shorthand for 
the polar form of a complex number: 

r(cos(0) +isin(0)) = re w . 
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Theorem 3.4.1. If z\, z 2 are two complex numbers, then 

e Zl e Z2 = e Zl+Z2 

Proof. 

e Zl e Z2 = e ai (cos(6i) + isin(&i)) • e Q2 (cos(6 2 ) + isin(6 2 )) 
= e ai+a ' 2 (cob(6i + 63) + i sui(&i +b 2 )) 
= e Zl+Z2 . 
Here we have just used ( ]3.8[ ) . □ 

The following theorem is amazing, since it involves calculus. 
Theorem 3.4.2. If w is a complex number, then 

^ wx wx 

— e = we , 
ax 

for x real. In fact, this is even true for x a complex variable (but we haven't 
defined differentiation for complex variables yet). 

Proof. Write w = a + bi. 

^ ^ax-\-bix 

dx 

-^(e ax {cos(bx) +isin(bx))) 
dx 

^-{e ax cos(frx) + ie ax sin(fcc)) 
dx 

^-{e ax cos(bx))+i-^(e ax sm{bx)) 
dx dx 

Now we use the product rule to get 

^-{e ax cos(bx))+i-^(e ax sm(bx)) 
dx dx 

— ae ax cos(fcc) — be ax sin(bx) + i(ae ax sin(&x) + be ax cos(bx)) 
= e ax (acos(bx) — bsm(bx) + i(asm(bx) + bcos(bx)) 

On the other hand, 

we wx = (a + bi)e ax+bxi 

= (a + bi)e ax (cos{bx) + i sin(6sc)) 
= e ax (a + bi)(cos(bx) + i sin(6x)) 

= e ax ((a cos(bx) — bsin(bx)) + i(asin(bx)) + bcos(bx)) 
Wow!! We did it! □ 
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That Theorem |3.4.2| is true is pretty amazing. It's what really gets complex 
analysis going. 

Example 3.4.1. Here's another amusing fact (if only for its obfuscating effect): 
1 = -e M . 

Solution. By definition, have e l7r — cos(7r) + isin(7r) = — 1 + iO = —1. 

3.4.1 Trigonometry and Complex Exponentials 

Amazingly, trig functions can also be expressed back in terms of the complex 
exponential. Then everything involving trig functions can be transformed into 
something involving the exponential function. This is very surprising. 

In order to easily obtain trig identities like cos(x) 2 + sin(x) 2 = 1, let's write 
cos(x) and sin(a;) as complex exponentials. From the definitions we have 

e lx = cos(x) + i sin(a;), 

so 

e~ ix = cos(— x) + isin(— x) = cos(x) — zsin(x). 

Adding these two equations and dividing by 2 yields a formula for cos(a;), and 
subtracting and dividing by 2i gives a formula for sin(;r): 

cos(a;) = — sin(x) = — . (3.9) 

We can now derive trig identities. For example, 

sin(2a;) = 

v ' 2i 

{e lx - e~ lx )(e lx + e- 



2i 

2 — = 2 sin(x) cos(a;) . 

I'm unimpressed, given that you can get this much more directly using 

(cos(2x) + i sin(2x)) = (cos(x) + i sin(a;)) 2 = cos 2 (cc) — sin 2 (x) + i2 cos(a;) sin(x), 

and equating imaginary parts. But there are more interesting examples. 
Next we verify that ( |3.9|) implies that cos(x) 2 + sin(a;) 2 = 1. We have 

4(cos(x) 2 + sin(a;) 2 ) = (e tx + e~ lx f + 

= e 2lx + 2 + e~ 2ix — (e 2lx — 2 + e~ 2vx ) = 4. 
The equality just appears as a follow- your-nose algebraic calculation. 
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1 




□.75 




0.5 




0.25 








-0.25 




0.5 




-0.75 




-1 






Figure 3.11: Plot of y — sin(x) 3 . 



Example 3.4.2. Compute sin(a;) 3 as a sum of sines and cosines with no powers. 
Solution. We use [3M : 

sin(x) 3 = 



2i 

3 



ix g — ix\3 



(e ix - e~ ix )(e ix - e~ ix )(e ix - e~ ix ) 



(e lx -e~ lx )(e Mx -2 + e~ MX ) 



3ix I „ix 



{e 6lx - 2e 



e lx + 2e~ lx - e~ Alx ) 



({e 6lx - e~ Mx ) - 3(e™ - e~ lx )) 



„ix —ix 



2i 



- 3 



2/ 



3 sin(a;) — sin(3x) 



You can also do this in Sage: 



sage: y = sin(x)~3 

sage: maxima (y ) .trigreduceO 

(3*sin(x)-sin(3*x))/4 



100 



3.5. INTEGRALS OF TRIGONOMETRIC FUNCTIONS 



3.5 Integrals of Trigonometric Functions 

There are an overwhelming number of combinations of trigonometric functions 
which appear in integrals, but fortunately they fall into a few patterns and most 
of their integrals can be found using reduction formulas and tables of integrals. 
This section examines some of the patterns of these combinations and illustrates 
how some of their integrals can be derived. 



/ sin(ax) sin(bx) dx, J cos(ax) cos(bx) dx, J sin(acc) cos(bx) dx. 
Products of Sine and Cosine: 

All of these integrals are handled by referring to the trigonometric identities 
for sine and cosine of sums and differences: 



sin (A + B) = sin(A) cos(-B) + cos(A) sin(JB) 
sm.{A - B) = sin(A) cos(B) - cos(A) sin(B) 
cos(A + B) = cos(A) cos(B) - sm(A) sin(B) 
cos(,4 - B) = cos(A) cos(S) + sm(A) sin(B). 

By adding or subtracting the appropriate pairs of identities, we can write the 
various products such as sin(aa;) cos(bx) as a sum or difference of single sines or 
cosines. For example, by adding the first two identities we get 2 sin(A) cos(-B) = 
sin(A + B) + s'm(A - B) so sin(^4) cos(-B) = 2sm(A + B) + sin(A - B). Using 
this last identity, the integral of sin(ax) cos{bx) for a ^ b is relatively easy: 

sin(aa;) cos{bx) dx = J — [sin((a + b)x) + sin((a — b)x)] dx 

1 .— cos((a — b)x) — cos((a + b)x) 
~ 2 a-b + a+b ' 

The other integrals of products of sine and cosine follow in a similar manner. 
For a 7^ b: 



J sin(ax) sin(6x) dx = 


1 

2 


.sin((a — b)x) 
a-b 


sin((o + Mx), 

— + 

a + b 1 


c, 


J 003(^)003(6,)^ = 


1 


r sin((a — b)x) 
— + 


sin((a + b)x) 


c, 


2 


a — b 


a + b J + 


J sin(ax) cos(bx) dx 




1 cos((a — b)x) 
2^ a-b 


cos((a + b)x) 

+ 

a + b 


+ C. 



This is confirmed by Sage: 
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sage: a,b = var("a,b") 

sage: integral (sin(a*x) *cos (b*x) ,x) 

-((b - a)*cos((b + a)*x) + (-b - a)*cos((b - a)*x))/(2*b~2 - 2*a~2) 
sage: integral (sin (a*x) *sin(b*x) ,x) 

-((b - a)*sin((b + a)*x) + (-b - a)*sin((b - a)*x))/(2*b~2 - 2*a~2) 
sage: integral (cos (a*x) *cos (b*x) ,x) 

((b - a)*sin((b + a)*x) + (b + a)*sin((b - a) *x) ) / (2*b~2 - 2*a~2) 



For a = b: 



/sin(a*)^ = 


x sin(aa;) cos(ax) 
2 2a 




J cos(ax) 2 dx = 


x sin(aa;) cos(ax) 
2 + 2a 




J sin(aa;) cos(ax) dx = 


s'm(ax) 2 „ 

\ +c. 

la 





The first and second of these integral formulas follow from the identities sin(ax) 2 = 
i-cos(2aa;) an( ^ cos(ax) 2 — 1 + co H 2az ) ; anc ] the third can be derived by a substi- 
tution using the variable to u = sin(aa;). These formulas too are confirmed by 
Sage: 



sage: a = var("a") 

sage: integral (cos (a*x) ~2,x) 

(sin(2*a*x) + 2*a*x)/(4*a) 

sage: integral (sin(a*x) ~2,x) 

-(sin(2*a*x) - 2*a*x)/(4*a) 

sage: integral (sin(a*x) *cos (a*x) ,x) 

-cos(a*x) ~2/(2*a) 



Remark 3.5.1. Note that Sage tells us that J sin(ax) cos(ax) dx = — c °2° + 

C but the table tells us that J sin(ax) cos(ax) dx = s ' n ^ t 2 ^ — h C . Aside from the 
ambiguity in the notation "+C", these are the same since sin(aa;) 2 = — cos(ax) + 
1. In other words, if you keep in mind that "+C" in one equation is not the 
same as "+C" in another, these formulas are the same. 
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Example 3.5.1. Compute J sin 3 (x)dx. 

We use trig, identities and compute the integral directly as follows: 



sin (x)dx = / sin (x) sin(x)dx 



[1 — cos 2 (cc)] sin(x)dx 

= — cos(.t) + - cos 3 (a;) + c (substitution u = cos(a;) ) 
o 

This idea always works for odd powers of sm{x). 

Example 3.5.2. What about even powers'? '! Compute J sm 4: (x)dx. We have 

sin 4 (a;) = [sin 2 (x)] 2 

1 — cos(2x)" 1 ' 



Thus 



sin i (x)dx = 



[l - 2cos(2.t) + cos 2 (2z)] 
1-2 cos(2a;) + - + - cos(4x] 



3 1 ,„ s 1 fl s 

cos(2a;) H — cos(4a;) 

8 2 v ; 8 y ' 



dx 



3 1 . ,„ \ 1 • r a \ 

= 8 X ~ 4 sm ( ) + 32 sm ^ 4x ) + c - 

Key Trick: Realize that we should write sin 4 (x) as (sin 2 (a;)) 2 . The rest is 
straightforward. 

Patterns for J sin(a;) m cos(x) n dx: If the exponent of sine is odd, we can split 
off one factor sin(x) and use the trig identity sin(a:) 2 = 1 — cos(.t) 2 to rewrite 
the remaining even power of sine in terms of cosine. Then the change of variable 
u = cos (a;) makes all of the integrals straightforward. If the exponent of cosine 
is odd, we can split off one factor cos(a;) and use the trig identity cos(a;) = 
1 — sin(a;) 2 to rewrite the remaining even power of sine in terms of cosine. Then 
the change of variable u — sin(x) makes all of the integrals straightforward. 
If both exponents are even, we can use the identities sin(x) 2 = 1 c °s( 2x ) an( j 
cos(a:) 2 = 1+co ^- 2x } to rewrite the integral in terms of powers of cos(2x) and 
then proceed with integrating even powers of cosine. 

Example 3.5.3. This example illustrates a method for computing integrals of 
trig functions that doesn't require knowing any trig identities at all or any tricks. 
It is very tedious though. We compute $ siv?(x)dx using complex exponentials. 
We have 
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cos(a;) 



~IX i „ — IX 



sin(x) 



2% 



sin 3 (x)dx 



dx 



1 

Yi 
1 

Yi 
1 

i 

1 

8> 



2i 

( e « _ e - zx )(e rx - e~ lx )(e lx - e~ lx )dx 
(e 2ix - 2 + e- 2lx )(e ix - e- ix )& 
e 3lx - e ix - 2e" + 2e 



ux _ e ~3ix + 3e - 



3e lx dx 



8i V 3i 
1 /l 



-3/ 



4 V 3 



cos(3x) — 3cos(a;) I + c 



1 3 

= — cos(3ai) — - cos(a;) + C. 

The answer looks totally different, but is in fact the same function. 

Example 3.5.4. The complex exponentials method used in the previous example 
also works for powers of different trig functions. For instance, 



sm 3 (x)cos 2 (x)dx = IC—^- ) 3 )(^^ f)dx 

-ie 5ix /32 + ie 3zx /32 + ie ix /16 - ie~ vx /16 - ie~ 3 "/32 + ie~ 5ix /32 

= -e 5tx /160 + e 3ix /96 + e"/16 + e~ ix /16 + e- 3ix /96 - e _5i! 7:L60 + C 

= (-e 5i 7l60 - e- 5ix /160) + (e 3 "/96 + e" 3ix /96) + (e ix /16 + e~ ix /W) + C 

= - cos(5a;) /80 + cos(3x) /32 + cos(.t) /8 + C. 

Here are some more identities that we'll use in illustrating some tricks below. 



dx 



tan(x) = sec 2 (x) 



(3.10) 



and 
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— sec(x) = sec(a;) tan(x). (3-11) 
dx 



Also, 



1 + tan 2 (:r) = sec 2 (x). (3.12) 



Example 3.5.5. Compute J t&n 3 (x)dx. We have 



tan(x) [sec 2 (2) — l] dx 
tan(x) sec 2 (x)dx — / ta.n(x)dx 



= - tan 2 (a;) — In | sec(x)| + c 

Here we used the substitution u — tan(a;) ; so du = sec 2 (x)dx, so 

f 1 1 

tan(x) scc 2 (x)dx = / udu = —u + c = — tan 2 (2;) + c. 

Also, with the substitution u = cos(x) and du = — sin(x)dx we get 
f f sin (x*) f 1 

/ t'&n(x)dx = / --—dx = — I —du = — In lul + c — — In I scc(a;)| + c. 

J J cos{x) J u 

Key trick: Write tan 3 (x) as t&n(x) tan 2 (a;). 

Example 3.5.6. Here's one that combines trig identities with the funnest vari- 
ant of integration by parts. Compute J sec 3 (x)dx. 
We have 



sec (x)dx = / sec(x)sec (x)dx. 



Let 's use integration by parts. 



u = scc(x) v = tan(x) 

du = sec(a;) t&n(x)dx dv = sec 2 (x)dx 
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The above integral becomes 

'sec(,)sec 2 (^ = sec(x)tan(,)-/sec(,)tan 2 (,)^ 

= B ec(x)taa(*)-/aec(x)[ fl ec>(x)-l]dx 
= sec(,)tan(,)-/sec3(,) + /sec(,)d, 
= s C c(,)tan(,)-/scc3(,) + ln| S ec(,) + tan(,)| 



This is familiar. Solve for J sec 3 (a;). We get 

sec 3 (x)dx = - sec(x) tan(x) + In | sec(x) + tan(x)| 

3.5.1 Some Remarks on Using Complex- Valued Functions 

Consider functions of the form 

f(x)+ig(x), (3.13) 
where x is a real variable and /, g are real- valued functions. For example, 

e lx = cos(x) + i sin(x). 

We observed before that 

4-e wx = we wx 
ax 

hence 

e wx dx = -e wx + c. 



w 

For example, writing it e lx as in (|3.13| , we have 



e lx dx = J cos(x)dx + i J sin(x)dx 

= sin(a;) — i cos(a;) + c 

= — i(cos(a;) + i sin(a;)) + c 

= l e lx . 

i 

f 1 

Example 3.5.7. Let's compute / dx. Wouldn't it be nice if we could just 

J x + i 

write ln(x + i) + c? This is useless for us though, since we haven't even defined 
hi(x + i)! However, we can "rationalize the denominator" by writing 
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-dx = 



-dx 



X + l X — I 

x — i 



x 2 + 1 

x 



.v- 



1 



dx 

dx — i 



-dx 



In \x 2 + 1| - ttan 1 (x) + c. 



This informs how we would define ln(z) for z complex (which you'll do if you 
take a course in complex analysis). 
Key trick: Get the i in the numerator. 

The next example illustrates an alternative to the method of Section 13.51 

Example 3.5.8. 



sin(5x) cos(3x)dx = 



1 

4i 
4i [lii 



2i 



e iSx _ e -iSx 



dx 



e t2x - e~ l2x ) dx 



8/ 



2i 



2/ 



cos(8a:) + cos(2a;) 



This is more tedious than the method in \3.5[ But it is completely straightfor- 
ward. You don't need any trig formulas or anything else. You just multiply it 
out, integrate, etc., and remember that i 2 = — 1. 
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Chapter 4 

Integration techniques 



4.1 Trigonometric Substitutions 

The first homework problem is to compute 

f 2 1 

Your first idea might be to do some sort of substitution, e.g., u = x — 1, 
but du — 2xdx is nowhere to be seen and this simply doesn't work. Likewise, 
integration by parts gets us nowhere. However, a technique called "inverse trig 
substitutions" and a trig identity easily dispenses with the above integral and 
several similar ones! Here's the crucial table: 



Expression 


Inverse Substitution 


Relevant Trig Identity 


V a 2 — x 2 


x = asin(0),-§ < 6 < f 


1 -sin 2 (0) =cos 2 (0) 


\/ a 2 + x 2 


x = atan(0),-f < 9 < f 


1 +tan 2 (0) =sec 2 (6») 


\J x 2 — a 2 


x = asec(0),O < 0<for7r<0<^ 


sec 2 (6>) - 1 =tan 2 (0) 



Inverse substitution works as follows. If we write x — g(t), then 



J f(x)dx = J f(g(t))g'(t)dt. 

This is not the same as substitution. You can just apply inverse substitution 
to any integral directly — usually you get something even worse, but for the 
integrals in this section using a substitution can vastly improve the situation. 

If g is a 1 — 1 function, then you can even use inverse substitution for a definite 
integral. The limits of integration are obtained as follows. 

,6 rg-Hb) 

/ f(x)dx= / f{g{t))g\t)dt. (4.2) 

J a J g 1 (a) 
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To help you understand this, note that as t varies from g^ 1 {a) to 5 _1 (6), the 
function g(t) varies from a = g{g^ 1 {a) to b = g{g~ 1 {b)), so / is being integrated 
over exactly the same values. Note also that ( ]4.2j ) once again illustrates Leibniz's 
brilliance in designing the notation for calculus. 
Let's give it a shot with ( 14. lj ). From the table we use the inverse substitution 



x = scc{9). 



We get 



y/2 x 3 \fx 2 — 1 



zdx 



3 1 



sec(6») 
s 1 



5 scc(0) 

cos ( 6)d6 
1 f f 



v/sec 2 (6») - lsec((9) tan((9)d(9 
tan(6») sec(6») tan(6>)d6> 



1 + cos(28)d8 



1 



sin(20) 



7T 1 

24 + ~8~ ~ 4 



Wow! That was like magic. This is really an amazing technique. Let's use it 
again to find the area of an ellipse. 

Example 4.1.1. Consider an ellipse with radii a and b, so it goes through 
(0,±6) and (±a, 0). An equation for the part of an ellipse in the first quadrant 
is 

V = b\ l—^— = - \J a 2 - x 2 . 
V a z a 

Thus the area of the entire ellipse is 

A = i[ - \J a? - x 2 dx. 
Jo a 

The 4 is because the integral computes l/Ath of the area of the whole ellipse. So 
we need to compute 

\J a 2 — x 2 dx 

Obvious substitution with u = a 2 — x 2 ...? nope. Integration by parts...? nope. 
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Let's try inverse substitution. The table above suggests using x = osin(^), so 
dx = acos(8)d8. We get 



a 2 - a 2 sin 2 (8)d6 = a 2 / cos 2 (9)d6 
o 



1 + cos{2d)d9 



2 Jo 



2 r 



Thus th 



e area is 



b ira 2 

A 

a 4 



2 2 



•nab. 



-sin(20) 



(4.3) 
(4.4) 
(4.5) 
(4.6) 



Consistency Check: If the ellipse is a circle, i.e., a = b = r , this is irr 2 , 
which is a well-known formula for the area of a circle. 

Remark 4.1.1. Trigonometric substitution is useful for functions that involve 
V 'a 2 — x 2 , V x 2 + a 2 , \J x 2 — a, but not all at once!. See the above table for how 
to do each. 



One other important technique is to use completing the square. 



Example 4.1.2. Compute J y/5 + Ax — x 2 dx. We complete the square: 
5 + Ax - x 2 = 5 - (x - 2) 2 + 4 = 9 - (x - 2) 2 . 



^5 + Ax - x 2 dx = / v/9 - (x - 2) 2 dx. 



Thus 



We do a usual substitution to get rid of the x — 2. Let u = x — 2, so du = dx. 
Then 



a/9 - (x - 2) 2 dx = / ^9 - V 2 dy 



Now we have an integral that we can do; it's almost identical to the previous 
example, but with a = 9 (and this is an indefinite integral). Let y = 3sin(#), so 
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dy = 3cos(6)d9. Then 

J V9 - (x - 2) 2 dx = J V9 - y 2 dy 



J ^3 2 - 3 2 sin 2 (6<)3 cos{8)d9 
9 J cos 2 (6>) d9 



1 + cos(28)d8 



Of course, we must transform back into a function in x, and that's a little tricky. 
Use that 

x — 2 = y = 3sin(0), 

so that 



sin 



x-2 



y/9 - {X - 2) 2 dx = ■ ■ ■ 



+ - sin(20) + C 



X 



sin 



3 

x-2 



sin(6*) cos(8) 
x - 2 



+ C 



2) 2 



C. 



Here we use that sin(20) = 2 sin(0) cos(0). Also, to compute cos(sin 1 (^y^)), 
we draw a right triangle with side lengths x — 2 and \/9 — [x — 2) 2 , and hy- 
potenuse 3. 



Example 4.1.3. Compute 

f 1 

J Vt 2 - 6i + 13 

To compute this, we complete the square, etc. 



dt 



Vi 2 -64+13 i v/(i-3) 2 + 4 
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(You may want to visualize a triangle with sides 2 and t — 3 and hypotenuse 
v/(<-3) 2 +4.J Then 

t-3 = 2 tan(0) 
y/(t- 3) 2 +4 = 2sec((9) 



cos(6>) 



(it = 2 sec 2 (0)d0 



5acfc to i/ie integral, we have 
1 



V / (i"3) 2 + 4 



di = 



2see 2 (fl) 
2sec(0) 

scc(6»)d6» 

ln|sec(0) + tan(0)| +C 

t-3 



In 



>/(* - 3) 2 + 42 + 



C. 



4.2 Integration by Parts 

The product rule is that 



dx 



[f(x)g(x)] = f(x)g'(x) + f'(x)g(x). 



Integrating both sides leads to a new fundamental technique for integration: 

f(x)g(x) = J f(x)g'(x)dx + J g(x)f'(x)dx. (4.7) 
Now rewrite ( |4.7j) as 

f{x)g'(x)dx = f(x)g(x) - / g(x)f'(x)dx. 



Shorthand notation: 



u = f(x) 
v = g(x) 



du = f'(x)dx 
dv = g ' (x)dx 



Then have 



udv = uv — I vdu. 



So what! But what's the big deal? Integration by parts is a fundamental 
technique of integration. It is also a key step in the proof of many theorems in 
calculus. 
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Example 4.2.1. J xcos(x)dx. 

u = x v = sin(x) 

du = dx dv = cos(x)dx 

We get 

xcos(x)dx = xsin(x) — I sin(x)c?a; = xsin(x) + cos(a;) + c. 



"Did this do anything for us?" Indeed, it did. 

Wait a minute — how did we know to pick u = x and v = sin(cc) ? We could 
have picked them other way around and still written down true statements. Let 's 
try that: 

u = cos(x) v = -x 2 

v ; 2 

du = — sm(x)dx dv = xdx 

1 , , fin 



x cos(x)dx = -x cos(x) + J -x sin(x)dx. 

Did this help!? NO. Integrating x 2 sin(x) is harder than integrating xcos(x). 
This formula is completely correct, but is hampered by being useless in this case. 
So how do you pick them? 

Choose the u so that when you differentiate it you get something sim- 
pler; when you pick dv, try to choose something whose antiderivative 
is simpler. 

Sometimes you have to try more than once. But with a good eraser nobody will 
know that it took you two tries. 

Question If integration by parts once is good, then sometimes twice is even 
better? Yes, in some examples (see Example \4-2.4 1 ) - But in the above example, 
you just undo what you did and basically end up where you started, or you get 
something even worse. 

i 

Example 4.2.2. Compute / sin _1 (a;)da;. Two points: 

Jo 

1. It's a definite integral. 

2. There is only one function; would you think to do integration by parts? 
But it is a product; it just doesn't look like it at first glance. 

Your choice is made for you, since we 'd be back where we started if we put 
dv = sin _1 (a;)c?2;. 

u = sin _1 (.x) v = x 

du = ; dv = dx 

VT^x* 
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We get 



o 



sin (x)dx = [a;sin _1 (a;)]g 



o \/l - x* 



zdx. 



Now we use substitution with w ~ 1 — x 2 , dw 

zdx = — — f w~*dw = — i/) 3 + c 



o Vl 



Ixdx, hence xdx 
■\/l-x 2 -\ 



\dw. 



Hence 



1 r , 

sm~ l (x)dx = [a;sin _1 (x)]g + v 1 — x 



7T V3 

h — 

12 2 



1 



But shouldn't we change the limits because we did a substitution? (No, since 
we computed the indefinite integral and put it back; this time we did the other 
option.) 

Is there another way to do this? / don't know. But for any integral, there 
might be several different techniques. If you can think of any other way to guess 
an antiderivative, do it; you can always differentiate as a check. 

Note: Integration by parts is tailored toward doing indefinite integrals. 

Example 4.2.3. This example illustrates how to use integration by parts twice. 
We compute 

x 2 e~ 2x dx 



1 



u — x 
du = 2xdx 



v = — -e 
dv = e 



2x 

- 2x dx 



We have 



c dx 



-x 2 e~ 2x 



"dx. 



Did this help? It helped, but it did not finish the integral off. However, we can 
deal with the remaining integral, again using integration by parts. If you do it 
twice, you what to keep going in the same direction. Do not switch your choice, 
or you '11 undo what you just did. 



-2.T 



dx 



U = X 

du = dx 
1 



1 

v = --e 

-2.v 



-xe 



-2.T 



1 

2~~ ' 2 
Now putting this above, we have 



' 2x dx 



dv = e 
1 



dx 



-2x 



-2x 



I- 



■ 2 e- 2x dx = -lx 2 e- 2x -\xe- 2x 



\e- 2x + c= -\e- 2x (2x 2 + 2x + 1) + c. 
4 4 
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Do you think you might have to do integration by parts three times? What if 
it were J x i e~ 2x dxl Grrr - you'd have to do it three times. 

Example 4.2.4. Compute J e x cos(x)dx. Which should be u and which should 

be v? Taking the derivatives of each type of function does not change the type. 
As a practical matter, it doesn't matter. Which would you prefer to find the 
antiderivative of? (Both choices work, as long as you keep going in the same 
direction when you do the second step.) 

u = cos(x) v = e x 

du = — sin(x)dx dv = e x dx 

We get 

Je* C os( X )d X = e*cos{ X ) + je*Mx)dx. 

We have to do it again. This time we choose (going in the same direction,) : 

u = sin(x) v = e x 

du = cos(x)dx dv = e x dx 

We get 



e x cos(x)dx = e x cos(x) + e x sin(a;) — / e x cos(x)dx. 

Did we get anywhere? Yes! No! First impression: all this work, and we're back 
where we started from! Yuck. Clearly we don't want to integrate by parts yet 
again. BUT. Notice the minus sign in front of J e x cos(x)dx; You can add the 
integral to both sides and get 

e x cos(dx) = e x cos(x) + e x sin(x) + c. 

Hence 

1 



e x cos(dx) = -e x (cos(x) + sin(x)) + c. 

4.2.1 More General Uses of Integration By Parts 

The Integration By Parts Formula is also used to derive many of the entries 
in the Table of Integrals. For some integrands such as x n ln(a;), the result is 
simply a function, an antiderivative of the integrand. For some integrands such 
as sin(x)™, the result is a reduction formula, a formula which still contains an 
integral, but the new integrand is the sine function raised to a smaller power, 
sin(a;) ra_2 . By repeatedly applying the reduction formula, we can evaluate the 
integral of sine raised to any positive integer power. 

Practice 4.2.1. Let n —1 be an integer. Evaluate J x n hi(x)dx using u = 
ln(x) and dv = x n dx. 

Practice 4.2.2. Let n =/= — 1 be an integer. Apply integration by parts to 
J x n e x dx using u = x n and dv — e x dx. 
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4.3 Factoring Polynomials 

How do you compute something like 

x 2 +2 



(x-l)(x + 2)(x + 3) dx7 

So far you have no method for doing this. The trick (which is called partial 
fraction decomposition), is to write 

+ 2 t f 1 2 11 , 

-dx = I — + — — dx (4.8) 



3 +4x 2 +x-6 J 4(i-l) x + 2 4(x + 3) 

The integral on the right is then easy to do (the answer involves ln's). 

But how on earth do you right the rational function on the left hand side as 
a sum of the nice terms of the right hand side? Doing this is called "partial 
fraction decomposition", and it is a fundamental idea in mathematics. It relies 
on our ability to factor polynomials and solve linear equations. As a first hint, 
notice that 

x 3 + Ax 2 + x - 6 = (x - 1) ■ {x + 2) ■ (x + 3), 

so the denominators in the decomposition correspond to the factors of the de- 
nominator. 

Before describing the secret behind ( |4.8|) , we'll discuss some background about 
how polynomials and rational functions work. 

Theorem 4.3.1 (Fundamental Theorem of Algebra). If f(x) — a n x n + - ■ ■ a±x+ 
cto is a polynomial, then there are complex numbers c, cti,...a n such that 

f(x) = c(x - ai){x - a 2 ) ■ ■ ■ (x - a„). 

Example 4.3.1. For example, 

3x 2 + 2x - 1 = 3 • (x - ij • (x + 1). 

And 

(x 2 + l) = {x + i) 2 -{x-i) 2 . 
If f(x) is a polynomial, the roots a of / correspond to the factors of /. Thus 

if 

f(x) = c(x - ai)(x -a 2 )---(x- a n ), 
then f(ai) = for each i (and nowhere else). 

Definition 4.3.1 (Multiplicity of Zero). The multiplicity of a zero a of f(x) is 
the number of times that {x — a) appears as a factor of f . 

For example, if f(x) = 7(x — 2)" • (x + 17) 5 • (x — it) 2 , then 2 is a zero with 
multiplicity 99, tt is a zero with multiplicity 2, and —1 is a "zero multiplicity 
0". 
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Definition 4.3.2 (Rational Function). A rational function is a quotient 

where g(x) and h{x) are polynomials. 
For example, 



^ ^ (x — i) 2 (x + ir)(x — 3) 3 ^ 



is a rational function. 



Definition 4.3.3 (Pole). A pole of a rational function f(x) is a complex number 
a such that \f(x)\ is unbounded as x — > a. 

For example, for ( ]4.9j) the poles are at i, tt, and 3. They have multiplicity 2, 
1, and 3, respectively. 



4.4 Partial Fractions 

Rational functions (polynomials divided by polynomials) and their integrals are 
important in mathematics and applications, but if you look through a table of 
integral formulas, you will find very few formulas for their integrals. Partly that 
is because the general formulas are rather complicated and have many special 
cases, and partly it is because they can all be reduced to just a few cases using the 
algebraic technique discussed in this section, Partial Fraction Decomposition. 
In algebra you learned to add rational functions to get a single rational function. 
Partial Fraction Decomposition is a technique for reversing that procedure to 
"decompose" a single rational function into a sum of simpler rational functions. 
Then the integral of the single rational function can be evaluated as the sum of 
the integrals of the simpler functions. 

Example 4.4.1. Use the algebraic decomposition ^J^-tx ~ ^ + 2^-5 ^° eva ^ ua t e 

f 17a35 J„ 
J 2x 2 -5x " X - 

Solution: The decomposition allows us to exchange the original integral for 
two much easier ones: 

f 17x35 , f 7 , f 3 , 

= 71n>| + |ln|2a: — 5| + C. 
When Sage computes this integral, it implicitly assumes that x > 0: 



sage: integral((17*x-35)/(2*x~2-5*x) ,x) 
3*log(2*x - 5)/2 + 7*log(x) 
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Note that Sage also leaves off the constant of integration, but this is more of 
an abbreviation than a matter of precision. 

Practice 4.4.1. Use the algebraic decomposition ^^Z^x-?, = 3x+i ^° 
evaluate J ^J^-^xS dx. 

The Example illustrates how to use a "decomposed" fraction with integrals, 
but it does not show how to achieve the decomposition. The algebraic basis 
for the Partial Fraction Decomposition technique is that every polynomial can 
be factored into a product of linear factors ax + b and irreducible quadratic 
factors ax 2 + bx + c (with b 2 — 4<xc < 0) . These factors may not be easy to find, 
and they will typically be more complicated than the examples in this section, 
but every polynomial has such factors. Before we apply the Partial Fraction 
Decomposition technique, the fraction must have the following form. 

(a) (the degree of the numerator) < (degree of the denominator); 

(b) the denominator has been factored into a product of linear factors and 
irreducible quadratic factors. 

If assumption (a) is not true, we can use polynomial division until we get a 
remainder which has a smaller degree than the denominator. If assumption (b) 
is not true, we simply cannot use the Partial Fraction Decomposition technique. 
Distinct Linear Factors 

If the denominator can be factored into a product of distinct linear factors, 
then the original fraction can be written as the sum of fractions of the form 
imca"factor • Our j ^ ^ s to nn d the values of the numbers in the numerators, and 
that typically requires solving a system of equations. 

Example 4.4.2. Find constants A and B such that ,} 7x3 Jl = — + . 

(2x — b)x x 2x — 5 

Solution: First, note the roots of the denominator are {0, 5/2}. Cross multiply: 

17x35 = A(2x - 5) + Bx. 
Now "kill B and solve for A using the first root x = 0: 

-35 = 17 • 035 = A(2 ■ - 5) + B ■ = -5A 
This gives A = 7. Next, "kill A and solve for B using the second root x = 5/2; 

17 • 5/235 = A(2 ■ 5/2 - 5) + B5/2 = 55/2. 

This gives B = 3, so 9 17 2 x3 - 5 = 7 - + 
In Sage this is very easy: 

sage: f = (17*x-35) / (2*x~2-5*x) 
sage: f .partial_f ractionO 
3/(2*x - 5) + 7/x 



119 



4.4. PARTIAL FRACTIONS 



Practice 4.4.2. Find values of A and B so -, — , 6 ,fw 7 „, = -4^ H — 

J [x+3){x — 2) x+3 x — 2 

In general, there is one unknown coefficient for each distinct linear factor of 
the denominator. However, if the number of distinct linear factors is large, we 
would need to solve a large system of equations for the unknowns. 

Practice 4.4.3. Using partial fractions, solve J ^+i)(x+3) ■ 

Other possible cases are listed as follows. (Of course, a rational function can 
involve more than one case as well.) 

• Distinct Irreducible Quadratic Factors 

If the factored denominator includes a distinct irreducible quadratic fac- 
tor, then the Partial Fraction Decomposition sum contains a fraction of 
the form of a linear polynomial with unknown coefficients divided by the 
irreducible quadratic factor: 



ax + b 
cx 2 + dx + e 

• Repeated Factors 

If the factored denominator contains a linear factor raised to a power 
(greater than one), then we need to start the decomposition with several 
terms. There should be one term with one unknown coefficient for each 
power of the linear factor. For example, 



ax + b 
(cx + d) s ' 

Here is the general procedure: 
Partial fraction decomposition of N(x) / 'D(x): Let N(x) be a polynomial 
of lower degree than another polynomial D(x). 

1. Factor D(x) into irreducible factors having real coefficients. Now D(x) is 
a product of distinct terms of the form (ax + b) r or (ax 2 + bx + c) s , for some 
integers r > 0, s > 0. For each term (ax + b) r the partial fraction decomposition 
of N(x) / D(x) contains a sum of terms of the form 

A\ A r 

1 H 

(ax + b) (ax + b) r 

for some constants Ai, and for each term (ax 2 + bx + c) r the PFD of N(x)/D(x) 
contains a sum of terms of the form 

Bix + Ci | | B s x + C s 
(ax 2 + bx + c) (ax 2 + bx + c) s 

for some constants Bi, Cj . \ is the sum of all these simpler rational functions. 
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2. Now you have an expression for which is a sum of simpler rational 
functions. The next step is to solve for these constants A's, B's, C"s occurring 
in the numerators. Cross multiply both sides by D(x) and expand out the 
resulting polynomial identity for N(x) in terms of the A's, B's, C's. Equating 
coefficients of powers of x on both sides gives rise to a linear system of equations 
for the A's, B's, C's which you can solve. 

Practice 4.4.4. Logistic Growth: The growth rate of many different populations 
depends not only on the number of individuals (leading to exponential growth ) but 
also on a "carrying capacity" of the environment. If x is the population at time 
t and the growth rate of x is proportional to the product of the population and 
the carrying capacity M minus the population, then the growth rate is described 
by the differential equation 



where k and M are constants for a given species in a given environment. 
Let k = 1 and M = 100, and assume the initial population is x(0) = 5. 

(a) Solve the differential equation ^ = kx(M — x), for x. 

(b) Graph the population x(t) for < t < 20. 

(c) When will the population be 20? 50? 90? 100? 

(d) What is the population after a "long" time? (Find the limit, as t becomes 
arbitrarily large, of x.) 

(e) Explain the shape of the graph in (a) in terms of a population of bacteria. 

(f) When is the growth rate largest? (Maximize dx/dt.) 

(g) What is the population when the growth rate is largest? 

Practice 4.4.5. Chemical Reaction: In some chemical reactions, a new ma- 
terial X is formed from materials A and B, and the rate at which X forms is 
proportional to the product of the amount of A and the amount of B remaining 
in the solution. Let x represent the amount of material X present at time t, and 
assume that the reaction begins with a grams of A, b grams of B, and no mate- 
rial X (x(0) = 0). Then the rate of formation of material X can be described 
by the differential equation 



Solve the differential equation for x if k = 1 and the reaction begins with (i) 7 
grams of A and 5 grams of B, and (ii) 6 grams of A and 6 grams of B. 



dx 



kx(M - x) 



dx 



k(a — x)(b — x). 
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4.5 Integration of Rational Functions Using Par- 
tial Fractions 

Our goal today is to compute integrals of the form 

P(x) 



■dx 



by decomposing / = q\^)) ■ This is called partial fraction expansion. 

Theorem 4.5.1 (Fundamental Theorem of Algebra over the Real Numbers). 
A real polynomial of degree n > 1 can be factored as a constant times a product 
of linear factors x — a and irreducible quadratic factors x 2 + bx + c. 

Note that x 2 + bx + c = (x — a) {x — a), where a — z + iw, a = z — iw are 
complex conjugates. 

Types of rational functions f(x) = quI . To do a partial fraction expansion, 
first make sure deg(P(x)) < deg(Q(x)) using long division. Then there are four 
possible situation, each of increasing generality (and difficulty) : 

1. Q{x) is a product of distinct linear factors; 

2. Q(x) is a product of linear factors, some of which arc repeated; 

3. Q(x) is a product of distinct irreducible quadratic factors, along with 
linear factors some of which may be repeated; and, 

4. Q{x) is has repeated irreducible quadratic factors, along with possibly 
some linear factors which may be repeated. 

The general partial fraction expansion theorem is beyond the scope of this 
course. However, you might find the following special case and its proof inter- 
esting. 

Theorem 4.5.2. Suppose p, q\ and q2 are polynomials that are relatively prime 
(have no factor in common). Then there exists polynomials ci\ and a.2 such that 

p ai Ot2 



<7i<72 Qi <72 

Proof. Since q± and q2 are relatively prime, using the Euclidean algorithm (long 
division), we can find polynomials si and S2 such that 

1 = siQi + s 2 q2- 
Dividing both sides by q±q2 and multiplying by p yields 

P _ Ol 02 

qi 92 qi <?2 ' 

which completes the proof. □ 
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Example 4.5.1. Compute 

x 3 - Ax - 10 , 

~~2 ^ dx - 

x z — x — o 

First do long division. Get quotient of x + 1 and remainder of 3x — A. This 
means that 

x 3 - Ax - 10 3x - A 

— — = x + 1 + — -. 

x z — x — o x z — x — o 

Since we have distinct linear factors, we know that we can write 



3x - 4 A B 

x 2 — x — 6 x — 3 x + 2 ' 



/or reaZ numbers A,B. A clever way to find A, B is to substitute appropriate 
values in, as follows. We have 

f{x)(x - 3) = ^4 =A + B- 

M A ; x + 2 x + 2 

Setting x = 3 on 6o£/i sides we /iaue (taking a limit): 

A = m = 3 -^ = 5 - = l. 
JK ' 3 + 2 5 



Likewise, we have 



B = f{ -2) = 3 -^- 4 = 2. 

J v 7 -2-3 



Thus 

x 3 - 4x - 10 , /" 1 
a!x = / x + 1 + 



x 2 — x — 6 7 x — 3 x + 2 

a; 2 + 2a; 



2 



21og|x + 2|+log|x-3|+C. 



Example 4.5.2. Compute the partial fraction expansion of ( - a ._ 3 ^ a ._|_ 2 ) 2 • ^ e 
partial fraction theorem, there are constants A,B,C such that 

x 2 ABC 



(x-3)(x + 2) 2 x - 3 x + 2 (x + 2) 2 

./Voie </iai there's no possible way this could work without the (x + 2) 2 term, since 
otherwise the common denominator would be (x — 3)(x + 2). We have 



A = [/(x)(x - 3)] x=3 = lx=3 



C=[f{x){x + 2f] x= 



x 2 _9_ 
(x + 2) 2 '^ 3 ~ 25' 
4 



x=-2 



5 
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This method will not get us B! For example, 

r 2 t 4- 2 d 

f(x)(x + 2) = -. ^ = A- 



(sc-3)(a: + 2) a: -3 x + 2' 

While true this is useless. 

Instead, we use that we know A and C , and evaluate at another value of x, 
say 0. 

9 B 4 



f (0) = = — H 1 ^ 

A ' -3 2 (2) 2 



so B = i| . TTius finally, 



_9_ 16 _4 

25 i 25 i 5 



(x-3)(x + 2) 2 7 ar-3 x + 2 (x + 2) 2, 

9 16 - 

= — In \x - 3| + — In \x + 2| + — ^— + constant. 
2o 25 x -|- 2 

Example 4.5.3. Let's compute J x -} +1 dx. Notice that x + 1 is a factor, since 
— 1 is a root. We have 

x a + 1 = (x + 1) (x 2 - x + 1) . 
There exist constants A,B,C such that 

1 A + C 

X 3 + 1 X + 1 x 2 — X + 1 

TTien 

A = /(*)(* + 1)1*=-! = i. 

You could find B,C by factoring the quadratic over the complex numbers and 
getting complex number answers. Instead, we evaluate x at a couple of values. 
For example, at x — we get 



so C = |. Aferf, use x = 1 to get B. 



/(I) 



I B(l) + 1 



so 



1 3 + 1 (1) + 1 (l) 2 - (1) + 1 
11 2 
2 = 6 +S+ 3' 



B=i_l_4 = _l 
6 6 6 3 
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Finally, 



1 



x 3 + 1 



dx 



x + 1 x 2 — x — 1 a; 2 — x — 1 
x - 2 



(7.r 



iln|* + l|-|/ 



x 2 — x + 1 



-dx 



remains to compute 



First, complete the square to get 



x-2 



x 2 — x + 1 



-dx. 



X — X + 1 = ( x — 



Let u=(x- \), so du ~ dx and x = u + h. Then 



,2.3 



udu 

,2,3 



= _ In 



3 2 _i / 2it 
- • — p= tan — = 
2 V3 V\/3 



- In |.x 2 - x + 1 1 - \/3 tan 1 ( 2X - 1 
2 1 1 V 73 



Finally, we put it all together and get 



x-2 



1 1 , 1 

— dx = — In x + 1 , 

x 3 + l 3 1 1 3 7 x 2 -x + l 

= - In |x + 1| In |x 2 — x+ l| 

3 6 1 1 



dx 



V3 t f2x- 1 

— tan" 1 =- 

3 V 73 



is 



Problem: Compute J cos 2 (x)e 3K dx using complex exponentials. The answer 

11 3 

--e~ 3:E + — e~ 3x sin(2x) - — e~ 3x cos(2x) + c. 
6 13 26 



Here's how to get it. 

cos 2 (x)e~ 3x dx 



2 + e" 



3 (2i-3)a: 2 



e" 3a; dx 



2i-3 3 6 



— 3x 



?(-2i-3)x 



= e 

6 



— 3x 



— 2i — 3 
2i - 3 ~ 2i + 3 
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Simplify the inside part requires some imagination 



= ^-{-2ie 2lx - 3e 2tx + 2ie- 2lx - 3e~ 2 ") 



2i - 3 2i + 3 13 

= — (4sin(2x) - 6cos(2x)) 



4.6 Improper Integrals 

Example 4.6.1. Make sense of J Q e~ x dx. The integrals 



t 

e~ x dx 

o 



make sense for each real number t. So consider 



lim / e x dx = lim [— e x ]l = 1. 

£ — >oo Jq t — >oo 



Geometrically the area under the whole curve is the limit of the areas for finite 
values oft. 

Example 4.6.2. Consider —p==dx. 

Problem: The denominator of the integrand tends to as x approaches the 
upper endpoint. Define 

1 1 f* 1 

dx = lim / dx 

o VI - x 2 t->i- Jo VI - x 2 

= lim (sin _1 (i) - sin" 1 ^)) = sin _1 (l) = - 



Here t — > 1 means the limit as t tends to 1 from the left. 

Example 4.6.3. There can be multiple points at which the integral is improper. 
For example, consider 

zdx. 



-oo^+X 2 

A crucial point is that we take the limit for the left and right endpoints inde- 
pendently. We use the point (for convenience only!) to break the integral in 
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half. 



1 ,o r~ 1 



1 + x 2 J_ 00 1+x 2 J 1 + x 2 

f° 1 /■* 1 

= lim / dx + lim / dx 

s->-ooJ g 1 + x 2 t-*ooJ Q I + X 2 

= lim (tan -1 (0) - tan _1 (s)) + lim (tan _1 (t) - tan _1 (0)) 

s — ► — oo t— »oc 

= lim (-tan _1 (s)) + lim (tan _1 (t)) 

s — > — oo t — >oo 

— 7T 7T 

= h - = 7T. 

2 2 

Example 4.6.4. Consider xdx. Notice that 

r° ■ C 

xdx = lim / xdx + lim / xdx. 

S — > — oo / t — >oo / n 

J s J 

TTiis diverges since each factor diverges independently. But notice that 

lim / xdx = 0. 

t->ooJ_ t 

This is not what f°° xdx means (in this course - in a later course it could be 

J — oo 1 

interpreted this way)! This illustrates the importance of treating each bad point 
separately ( since Example \4--6.3) doesn't. 



Example 4.6.5. Consider J j -h^dx. We 



have 



1 1 f s f 1 

—j=dx= lim / x~ 5 dx + lim / x _s dx 
.! s/x s^o- J_! t^o+ J t 

,. /3 2 3\ ,. /3 3 2 \ 

= hm -S3 + lim 1* = 0. 

s^o-\2 2 J t^o+ \2 2 

This illustrates how to be careful and break the function up into two pieces when 
there is a discontinuity. 

Example 4.6.6. Compute J"_ -^rz^dx. A few weeks ago you might have done 
this: 



3 1 



_i x- 2 



dx = [In |x - 2\] 3 _! = ln(3) - ln(l) (totally wrong!) 



This is not valid because the function we are integrating has a pole at x = 2. 
The integral is improper, and is only defined if both the following limits exists: 



4 i . . .. r 3 i 



lim / dx and lim / dx. 

s->2- J_i x — 2 t^2+ J t x — 2 
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However, the limits diverge, e.g., 

,3 j 

lim / dx = lim (In 111 - In It — 21) = - lim In It - 21 = -oo. 

t-»2+ J t X — 2 t^2+ t— >2+ 

TTius J j -^rz^dx is divergent. 

4.6.1 Convergence, Divergence, and Comparison 

In this section we discuss using comparison to determine if an improper integrals 
converges or diverges. Recall that if / and g are continuous functions on an 
interval [a, b] and g(x) < f(x), then 

b pb 

g(x)dx < I f(x)dx. 



This observation can be incredibly useful in determining whether or not an 
improper integral converges. 

Not only does this technique help in determining whether integrals converge, 
but it also gives you some information about their values, which is often much 
easier to obtain than computing the exact integral. 

Theorem 4.6.1 (Comparison Theorem (special case)). Let f and g be contin- 
uous functions with < g(x) < f{x) for x > a. 

1. If J°° f{x)dx converges, then g{x)dx converges. 

2. If J a °° g(x)dx diverges then f{x)dx diverges. 
Proof. Since g{x) > for all x, the function 

G{t) = [ g{x)dx 

J a 

is a non-decreasing function. If J a °° f(x)dx converges to some value B, then for 
any t > a we have 

ft r t 



G(t) = [ g(x)dx < [ f{x)dx < B. 

J a J a 



Thus in this case G(t) is a non-decreasing function bounded above, hence the 
limit linif^oo G(t) exists. This proves the first statement. 
Likewise, the function 

F(t)= f f(x)dx 

J a 

is also a non-decreasing function. If J a °° g(x)dx diverges then the function G(t) 
defined above is still non-decreasing and limbec G{t) does not exist, so G{t) is 
not bounded. Since g(x) < f(x) we have G(t) < F(t) for all > a, hence F(t) is 
also unbounded, which proves the second statement. □ 
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The theorem is very intuitive if you think about areas under a graph. "If the 
bigger integral converges then so does the smaller one, and if the smaller one 
diverges so does the bigger ones." 

Example 4.6.7. Does f™ ^f^-dx converge? Answer: YES. 
Since < cos 2 (a;) <l, we really do have 

< cos 2 (x) < 1 
~ 1 + x 2 ~ 1 + x 2 ' 



Thus 



f°° 1 _ 7T 

/ ^dx = lim tan _1 (t) = — , 

J 1 + x 2 t^oo w 2 



roc cob 2 (x) 7 

so J t A ai converges. 

But why did we use 1 _^ x2 ? It 's a guess that turned out to work. You could 
have used something else, e.g., -% for some constant c. This is an illustration 
of how in mathematics sometimes you have to use your imagination or guess 
and see what happens. Don't get anxious — instead, relax, take a deep breath and 
explore. 

For example, alternatively we could have done the following: 

cos 2 (X) , f°° 1 , 

K -^-dx< / — dx = l, 

1 + £ J i x z 

and this works just as well, since c ° s +x 2 - dx converges (as c °^^ is continu- 
ous). 

Example 4.6.8. Consider z , e 1 _ 2 x dx. Does it converge or diverge? For 
large values of x, the term e~ 2x very quickly goes to 0, so we expect this to 
diverge, since \dx diverges. For x > , we have e~ 2x < 1, so for all x we 
have 

> (verify by cross multiplying). 



1 



But 



f 1 

/ dx = lim [Info; + l)]i = oo 

J l X + 1 t^oo 

Thus J °° x i dx must also diverge. 

Note that there is a natural analogue of Theorem 14.6.11 for integrals of functions 
that "blow up" at a point, but we will not state it formally. 

Example 4.6.9. Consider 

fl e -x pi e -x 

—f=dx = lim / — -j=dx. 
o V x t^o+J t ^x 

We have 

e~ x < 1 
129 



4.6. IMPROPER INTEGRALS 



( Coming up with this comparison might take some work, imagination, and trial 
and error.) We have 

f 1 e~ x f 1 1 

/ -^dx < / — dx = lim 2V1 - 2Vt = 2. 
Jo V x Jo V x f ^ 0+ 

thus J Q ^j^dx converges, even though we haven't figured out its value. We just 
know that it is < 2. (In fact, it is 1.493648265 ....) 

What if we found a function that is bigger than and its integral diverges?? 
So what! This does nothing for you. Bzzzt. Try again. 

Example 4.6.10. Consider the integral 




This is an improper integral since f{x) = ^— has a pole at x = 0. Does it 
converge? NO. 

On the interval [0, 1] we have e~ x > e _1 . Thus 

f 1 e~ x . f 1 e _1 
lim / dx > lim / dx 

i^0+ J t X t^0+ J t X 

= e • lim / —dx 

*->o+ Jt x 

= eT 1 ■ lim ln(l) - ln(t) = +oo 

Thus Jg 1 ^-dx diverges. 
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Chapter 5 

Sequences and Series 



Our main goal in this chapter is to gain a working knowledge of power series and 
Taylor series of function with just enough discussion of the details of convergence 
to get by. 



5.1 Sequences 

What is 



lim — ? 

n — >oc 2 n 



You may have encountered sequences long ago in earlier courses and they 
seemed very difficult. You know much more mathematics now, so they will 
probably seem easier. On the other hand, we're going to go very quickly. 

A sequence is an ordered list of numbers. These numbers may be real, complex, 
etc., etc., but in this book we will focus entirely on sequences of real numbers. 
For example, 

1111111 1 

2'4'8' 16' 32' 64' 128' " ' ' 2"' " ' 
Since the sequence is ordered, we can view it as a function with domain the 
natural numbers = 1, 2, 3, . . .. 

Definition 5.1.1 (Sequence). A sequence {a n } is a function a : N — ► R that 
takes a natural number n to a n = a(n). The number a n is the nth term. 

For example, 

a(n) = a n = — , 
which we write as {^}- Here's another example: 



(bn) 



oo 



1 2 3 



i;„ =1 2'3'4' 



131 



5.2. SERIES 



Example 5.1.1. The Fibonacci sequence (F n )°° =1 is defined recursively as fol- 
lows: 

F 1 = 1, F 2 = 1, F n = F n - 2 + F„_! /or n > 3. 

Let's return to the sequence (^r) We write lining ^ = 0, since the 
terms get arbitrarily small. 

Definition 5.1.2 (Limit of sequence). // (a n )'^' =1 is a sequence then that se- 
quence converges to L, written linin^oo a n = L, if a n gets arbitrarily close to 
L as n get sufficiently large. SECRET RIGOROUS DEFINITION: For every e > 
there exists B such that for n > B we have \a n — L\ < e. 

This is exactly like what we did in the previous course when we considered 
limits of functions. If f{x) is a function, the meaning of lim^^oo f{x) = L is 
essentially the same. In fact, we have the following fact. 

Proposition 5.1.1. // / is a function with lim.j.^00 /(a;) = L and (<!„) ■, is 
the sequence given by a n = f(n), then limn—^ a„ = L. 

As a corollary, note that this implies that all the facts about limits that you 
know from functions also apply to sequences! 



Example 5.1.2. 



lim = lim = 1 

n— >oo n + 1 x^oo x + 1 



Example 5.1.3. The converse of Proposition \ 5.1A\ is false in general, i.e., 
knowing the limit of the sequence converges doesn't imply that the limit of the 
function converges. We have \\m n -too cos (27m) = 1, but lim :1 ,_ +00 cos(27rx) di- 
verges. The converse is OK if the limit involving the function converges. 

Example 5.1.4. Compute lim 5 . Answer: i. 

rwoo Yin 6 — 2006n + 15 ' 

5.2 Series 

What is 
What is 



What is 



1 


1 


1 


1 


1 


2 




f 8- 


-16' 


"32- 




1 


1 


1 


1 


; + 


9 + 


27" 


"ST 


" 243 


1 


1 


1 


1 


1 


1 


f 4" 


f 9- 


-16- 


"25- 



Consider the following sequence of partial sums 

N 



^ 1 11 
a N = 2-.^= + 2 



n=l 



1 

2' 1 ~~ 2 ' 4 ' " + W ' 
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Can we compute 



oo 

Y- 



? 

2" 



n=l 

These partial sums look as follows: 

1 3 1023 1048575 

ai ~r 02 "4' ai °- 1024' 020 " 1048576 

oo ^ 

It looks very likely that — = 1, if it makes any sense. But does it? 

n=l 

In a moment we will define 



N 



— = lim — = lim om. 

2" JV->oo ^ 2™ JV->ao 
n—1 n — 1 

A little later W6 will sllOW tllclt CljsJ — — TTff -, hence indeed Y^Li g^r = !• 

Definition 5.2.1 (Sum of scries). If '(a n )^° =1 is a sequence, then the sum of the 

series is 

oo N 

Ea n = lim > a n = lim sm 

n—1 n—1 

provided the limit exists. Otherwise we say that Yl^L-f a n diverges. 

Example 5.2.1 (Geometric series). Consider the geometric series X^Li ar ™ _1 
for a^O. Then 



sn — > or 



1 — r 



E 1 



To see this, multiply both sides by 1 — r and notice that all the terms in the 

a(l-r N 
->oo 1—r 



middle cancel out. For what values of r does limjv-joo _ r — - converge? If 



\r\ < 1, then limjv^oo r N = and 

,. a(l-r w ) a 
lira = . 

Ar->oo 1 — r 1—r 

// |r| > 1, then liniAr^oo r N diverges, so X^^Li ar ™ _1 diverges. If r = ±1. 
it's clear since a ^ i/ia£ i/ie series also diverges (since the partial sums are 
s N = ±Na). 
For example, if a = 1 and r = 1 , we get 



oo 

i r ,-. 



as claimed earlier. 

133 



5.3. THE INTEGRAL AND COMPARISON TESTS 



5.3 The Integral and Comparison Tests 

What is E^=! sr? What is E,T=i 

Recall that Section [572] began by asking for the sum of several series. We found 
the first two sums (which were geometric series) by finding an exact formula for 
the sum of the first ./V terms. The third series was 

°° 1 1111 1 
A = £^ = T + 4 + 9 + 16 + 25 + ---- (5 ' 1} 

n—l 

It is difficult to find a nice formula for the sum of the first n terms of this series 
(i.e., I don't know how to do it). 

Remark 5.3.1. Since I'm a number theorist, I can't help but make some further 
remarks about sums of the form (5. 1\) . In general, for any s > 1 one can consider 
the sum 

oo 

coo = E-r- 

n=l 

The number A that we are interested in above is thus C(2)- The function Q{s) is 
called the Riemann zeta function. There is a natural (but complicated) way of 
extending £(s) to a (differentiable) function on all complex numbers with a pole 
at s = 1. The Riemann Hypothesis asserts that if s is a complex number and 
C(s) = then either s is an even negative integer or s = \ + bi for some real 
number b. This is probably the most famous unsolved problems in mathematics 
(e.g., it's one of the Clay Math Institute million dollar prize problems) . Another 
famous open problem is to show that £(3) is not a root of any polynomial with 
integer coefficients (it is a theorem of Apeery that zeta(3) is not a fraction). 

The function £(s) is incredibly important in mathematics because it governs 
the properties of prime numbers. The Euler product representation of C(s) gives 
a hint as to why this is the case: 

«*>=£^= n 

n— 1 primes p 

To see that this product equality holds when s is real with Re(s) > 1, use Exam- 
ple \5.2.1\ with r = p~ s and a = 1 above. We have 



= 1 + p~ s + p~ 2s + ■■■ . 
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Thus 

1 1 



n u== - ii 



oo 1 
n=l 



primes p 

11 \ / 1 1 

1 1 1 

= |1+ 2^ + 3^ 



where the last line uses the distributive law and that integers factor uniquely as 
a product of primes. 

of \C( S )\ f or complex s. 



This section is how to leverage what you've learned so far in this book to say 
something about sums that are hard (or even "impossibly difficult" ) to evaluate 
exactly. For example, notice (by considering a graph of a step function) that if 
fix) = l/x 2 , then for positive integer t we have 



Thus 



x - 1 1 

^— ' n 2 

n=l 



l 2 



' 1 



rdx. 



v — < — r —dx 

2— 1 n 2 ~ l 2 A x 2 

n=l J 1 



= 1+ lim 



1 + lim 

t— +OC 



1 






rdx 


X 1 


r 


t 


X 


1 


1 


1 






t 


+ T 



= 2 



We conclude that JZ^Li converges, since the sequence of partial sums is getting 
bigger and bigger and is always < 2. And of course we also know something 
about Y^Li ^2 even though we do not know the exact value: J^^Li ^? — 2- 
Using a computer we find that 
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t 




1 


i 


2 


I = 1-25 


5 


m = 1.4636T 


10 


™* = 1.54976773117 


100 


1.63498390018 


1000 


1.64393456668 


10000 


1.64483407185 


100000 


1.6449240669 



The table is consistent with the fact that V°° , \ converges to a number < 2. 
In fact Euler was the first to compute J2™=i exactly; he found that the exact 
value is 



7i~~ 

~6 



1.644934066848226436472415166646025189218949901206798437735557. 



There are many proofs of this fact, but they don't belong in this book; you 
can find them on the internet, and are likely to see one if you take more math 
classes. 

We next consider the harmonic series 



oo 

E 

71=1 



(5.2) 



Does it converge? Again by inspecting a graph and viewing an infinite sum as 
the area under a step function, we have 



E 

n=l 



1 



> 



1 



-dx 



= lim [ln(x)][ 

= lim ln(£) - = +00. 

t — >oo 



Thus the infinite sum (|5.2|) must also diverge. 

We formalize the above two examples as a general test for convergence or 
divergence of an infinite sum. 

Theorem 5.3.1 (Integral Test and Bound). Suppose f(x) is a continuous, pos- 
itive, decreasing function on [1, oo) and let a n = f(n) for integers n > 1. Then 
the series 2r»=i a ™ converges if and only if the integral J, f{x)dx converges. 
More generally, for any positive integer k, 



f 

Jk 



f{x)dx 



< 



E< 

n=k 



< 



f{x)dx. 



(5.3) 
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5.3. THE INTEGRAL AND COMPARISON TESTS 



The proposition means that you can determine convergence of an infinite series 
by determining convergence of a corresponding integral. Thus you can apply the 
powerful tools you know already for integrals to understanding infinite sums. 
Also, you can use integration along with computation of the first few terms of 
a series to approximate a series very precisely. 

Remark 5.3.2. Sometimes the first few terms of a series are "funny" or the 
series doesn't even start at n = I, e.g., 



^ (n-3) 3 ' 

n=4 v ' 



In this case use (5. 3\) with any specific k > 1. 

Proposition 5.3.1 (Comparison Test). Suppose ^ a„ and^2b n are two series 
with positive terms. If J^6 n converges and a n < b n for all n. then ^a„ 
converges. Likewise, if b n diverges and a n > b n for all n. then ^ a n must 
also diverge. 

Example 5.3.1. Does J^^Li ^ converge? No. We have 

f 00 1 

> / —dx = lim (2Vt - 2VT) = +oo 
Ji Vx t^oo 



oc 1 



Example 5.3.2. Does X^^Li n ^+i converge? Let's apply the comparison test: 
we have ?t 2 +1 < for every n, so 



^ n 2 + 1 ^ n 2 ' 

n—l n—1 



Alternatively, we can use the integral test, which also gives as a bonus an upper 
and lower bound on the sum. Let f(x) = 1/(1 + x 2 ). We have 

1 /"* 1 

dx — lim / Trdx 



1 + x 2 t->oo J l 1 + x 2 



1 7T 7T 7T 7T 

hm tan (t) = = — 

t^oo w 4 2 4 4 



Thus the sum converges. Moreover, taking k = 1 in Theorem \5.3.1\ we have 

7T \ 1 1 7T 

- < > < - + -. 

4-^n 2 + l~2 4 

n— 1 

the actual sum is 1.07 . . which is much different than ^ ^2 = 1-64 .... 

We could prove the following proposition using methods similar to those illus- 
trated in the examples above. 

Proposition 5.3.2. The series Y^=i * s convergent if p > 1 and divergent 
ifp < I- 
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5.3.1 Estimating the Sum of a Series 

Suppose ^2 a n is a convergent sequence of positive integers. Let 



Or, 
n—m+1 



Rm — ^ CLn — ^ a n — ^ 
n— 1 n— 1 n— m+ 

which is the error if you approximate ^ a„ using the first n terms. From 
Theorem 5.3.1 we get the following. 

Proposition 5.3.3 (Remainder Bound). Suppose f is a continuous, positive, 
decreasing function on [m, oo) and ^ a n is convergent. Then 

/>oo 

f(x)dx < R m < / f(x)dx. 

m+l m 



Proof. In Theorem 5.3.1 set k = m + l. That gives 

/oo °° />oo 

a n < flm+i + / f(x)dx. 
.1+1 n=m+l "' m + 1 

But 

/•OO /• OO 

a m +i + / f(x)dx < / f{x)dx 

J m+l J m 

since / is decreasing and f{m + 1) = a m +\. □ 

Example 5.3.3. Estimate C,('i) = YlnLi ^ using the first 10 terms of the 
series. We have 

n— 1 

TTie proposition above with m = 10 ieZZs ms £/iai 

f 00 1 10 f°° 1 11 

0.00413223140495867 . . . = / -^dx < C(3l- V < / — dx = ^ = = 0.005. 

in z 3 " W Ao ^ 3 2-10 2 200 

in /act, 

C(3) = 1.202056903159594285399738161511449990. . . 

and we have 

10 

(•(3) - ^ = 0.0045249174854010 . . . , 

n=l 

so i/*e integral error bound was really good in this case. 

Example 5.3.4. Determine if ii7n^+°Ain+3 converges or diverges. An- 

swer: It converges, since 

2006 2006 _ 2006 1 

117n 2 + 41n + 3 ~ 117n 2 ~ TT7 ' n 5 ' 

and ^2 ^2 converges. 
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5.4 Tests for Convergence 

5.4.1 The Comparison Test 

Theorem 5.4.1 (The Comparison Test). Suppose ^a„ and Y] b n are series 
with all a n and b n positive and a n < b n for each n. 

1. If^2b n converges, then so does Y2 a n- 

2. If ^2 a n diverges, then so does b n ■ 

Proof Sketch. The condition of the theorem implies that for any k, 

k k 
a„ < ^ b n , 

n—1 n—1 

from which each claim follows. □ 

Example 5.4.1. Consider the series Vl 00 ^ Q 2 7 , . For each n we have 

7 7 1 

3n 2 + 2n - 3 ' r? ' 

Since X^^Li converges, Theorem \5.4-l\ implies that Yl^—i z n ^+2n a ^ so con ~ 
verges. 

Example 5.4.2. Consider the series Y2n=i '"i"^ • ^ diverges since for each 
n > 3 we have 

ln(n) 1 
n — n ' 

tm<i Y2n=3 n diverges. 

5.4.2 Absolute and Conditional Convergence 

Definition 5.4.1 (Converges Absolutely). We say that YL^Li a n converges ab- 
solutely if Y2n=i l a «J converges. 

For example, 

n—1 

converges, but does noi converge absolutely (it converges "conditionally", though 
we will not explain why in this book). 
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5.4. TESTS FOR CONVERGENCE 



5.4.3 The Ratio Test 



Recall that X^Li a » * s a geometric series if and only if a n = ar n ~ x for some 
fixed a and r. Here we call r the common ratio. Notice that the ratio of any 
two successive terms is r: 

a n+1 ar n 



r-n— 1 



Moreover, we have Y^Li ar ™ 1 converges (to jzr) if and only if |r| < 1 (and, 
of course it diverges if |r| > 1). 

Example 5.4.3. For example, 5Z^=i 3 (§)" 1 converges to = 9. However, 
Er=i 3 (l) n_1 diverges. 

Theorem 5.4.2 (Ratio Test). Consider a sum Y^=i a n- Then 



1. 7/lim„^ c 

2. 7/lim„^ c 

3. 7/lim„^ c 



O-n + 1 



gn+1 



= _L < 1 i/ien E^Li °" * s absolutely convergent. 

= L > 1 i/ien a « diverges. 

= L = 1 i/ien we may conclude nothing from this! 



Proof. We will only prove 1. Assume that we have lim,- 



L < 1. Let 



_ i+i 



and notice that L < r < 1 (since < L < 1, so 1 < L + 1 < 2, so 



1/2 < r < 1, and also r-L= (L+ l)/2 - L = (1 - L)/2 > 0). 



Since lim»- 



gra+1 



= L, there is an N such that for all n > N we have 



a,: 



< r, so |a n+ i| < |a„| ■ r. 



Then we have 

OO CO 

53 l<»n| < \a N +i\-J2 rn - 

n—N+1 n=0 

Here the common ratio for the second one is r < 1, hence thus the right-hand 
series converges, so the left-hand series converges. □ 



Example 5.4.4. Consider 



(-10) 



n+l 



(71+1)! 

(-10)" 



(-10)" 



10 n+l 



TTie ratio of successive terms is 



10 



(n + l)n! 10" n + l 



< 1. 



TTius this series converges absolutely. Note, the minus sign is missing above 
since in the ratio test we take the limit of the absolute values. 
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5.4. TESTS FOR CONVERGENCE 



Example 5.4.5. Consider 1+3tt - We have 



(n + l) n+1 
3 • (27)" +1 



il+3n 



n=l 

(n + l)(n + l)" 27 ,; 
27 • 27" ' n™ 



n + 1 
27 



77 + 1 



- 00 



Thus our series diverges. (Note here that we use that (H±l)" — ► 
Example 5.4.6. Lei's appZy £/ie ratio test to „• ^ e have 

1 

1 77, 77, 

1. 



lim 

n — >oc 



77+1 



1 

n 



n + 1 1 77+1 



This tells us nothing. If this happens... do something else! E.g., in this case, 
use the integral test. 

5.4.4 The Root Test 

Since e and In are inverses, we have x = e ln ' x '. This implies the very useful fact 
that 

x a 



As a sample application, notice that for any nonzero c, 



lim c r 

n — >oo 



lim e 

n — »oc 



1 log(c) = e 



1. 



Similarly, 



where we've used that lim„_> 
rule. 



lim 77- = lim e" los(n) = e° = 1 

n — >oo n — >oo 

log(») 



n 



0, which we could prove using L'Hopital's 



Theorem 5.4.3 (Root Test). Consider the sum Y^=i a «- 

1. J/liixin^oo |a„|" = L < 1, then X^^Li a n converges absolutely. 

2. i71im„_> 00 |a„|« = Z > 1, £/ie77 Y^n=i a « diverges. 

3. If L = 1, then we may conclude nothing from this! 

Proof. We apply the comparison test (Theorem 1 5.4.1}) . First suppose linin^oo \a n \ 
L < 1. Then there is a such that for n > N we have |a„|« < k < 1. Thus 
for such n we have \a n \ < k n < 1. The geometric series Y^Ln ^ converges, so 
J^i^N l a »l a ^ S0 does, by Theorem 15.4.11 If \a n \ n > 1 for 77 > TV, then we see 
that J2i1n l a «l diverges by comparing with J2i^N 1=1 
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5.4. TESTS FOR CONVERGENCE 



Example 5.4.7. Let's apply the root test to 

oo oo 
n—l n—1 

We have 



lim |r n |s 

n — >oc 



Thus the root test tells us exactly what we already know about convergence of 
the geometry series (except when \r\ = 1). 



Example 5.4.8. The sum Y^Li 
have 

n 2 + 1 



2k 2+1 J is a candidate for the root test. We 



lim 

2n 2 + 1 
Thus the series converges. 



1 



= lim — — = lim ^ - 

n->oo 2n 2 + 1 «^oo 2 + X 2 



Example 5.4.9. The sum J2^Li 
have 

'2n 2 + 1 



oo / 2n-* + l 
n 2 + l 



is a candidate for the root test. We 



lim 

n — >oo 



n 2 + 1 



lim 



2n 2 + 1 



2 + \ 
lim = 2. 



n— >oo n 2 + 1 n— >oo 1 -| \ 



hence the series diverges! 

Example 5.4.10. Consider Y^n=i ^ • We have 



lim 

n — >oo 



1. 



so we conclude nothing! 

Example 5.4.11. Consider Y^Li 3-{27 n ) ■ To apply the root test, we compute 

'1\ " n 



lim 



3 • (27™) 



= lim 



27 



= +00. 



Again, the limit diverges, as in Example \5.4-5\ 
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5.5. POWER SERIES 



5.5 Power Series 

Recall that a polynomial is a function of the form 

f(x) = c + ax + c 2 x 2 H h c k x k 



Polynomial^ar^^as^J 



They are easy to integrate, differentiate, etc. 

' k \ k 



d 

dx 



Vn— / n— 1 



/ fc k ™n+l 
Z = c + Z c » irrr 
n n 



n=0 n=0 

Definition 5.5.1 (Power Series). A power series is a series of the form 



f(x) = y~] c„a;™ = co + cix + c 2 x 2 



where x is a variable and the c n are coefficients. 

A power series is a function of x for those x for which it converges. 
Example 5.5.1. Consider 



f(x) = ]T x n = 1 + x 
When \x\ < 1, i.e., —1 < i < 1, i«c /iai>e 

1 



/(*) = 



1-x 

But what good could this possibly be? Why is writing the simple function 
as the complicated series X^^Lo a; ™ °^ an y va l ue ? 

f . Power series are relatively easy to work with. They are "almost" polyno- 
mials. E.g., 

, oo oo oo 

da Z = nxn ~ x = 1 + 2x + 3x2 + • • ■ = H ( m + 1 ) a;m ' 

n— n— 1 m— 

where in the last step we "re-indexed" the series. Power series are only 
"almost" polynomials, since they don't stop; they can go on forever. More 
precisely, a power series is a limit of polynomials. But in many cases we 
can treat them like a polynomial. On the other hand, notice that 

1 \ 1 00 

Y^(m + l)x m . 



dx \1 — x J (1 — a;) 2 

v ' y ' m=0 
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5.5. POWER SERIES 



2. For many functions, a power series is the best explicit representation avail- 
able. 

Example 5.5.2. Consider Jq(x), the Bessel function of order 0. It arises 
as a solution to the differential equation x 2 y" + xy' + x 2 y = 0, and has 
the following power series expansion: 

(-l) n x 2n 

2~ 



M') = £ 2 2n (n!) 

n=l v ' 



= 1 - -x 2 + —x 4 — x 6 + — 1 — .t 8 1 x 10 H . 

4 64 2304 147456 14745600 

This series is nice since it converges for all x (one can prove this using 
the ratio test). It is also one of the most explicit forms of Jq(x). 

5.5.1 Shift the Origin 

It is often useful to shift the origin of a power series, i.e., consider a power series 
expanded about a different point. 

Definition 5.5.2. The series 

oo 

£ c n (x — a) n = c + c\{x — a) + c 2 (x — a) 2 + ■ ■ ■ 
n=0 

is called a power series centered at x = a, or "a power series about x = a". 
Example 5.5.3. Consider 

oo 

Y^{x- 3)" = 1 + (x - 3) + (x - 3) 2 + ■ ■ ■ 

n=0 

1 



1- (x- 3) 
1 



equality valid when \x — 3| < 1 



4 - x 



Here conceptually we are treating 3 like we treated before. 

Power series can be written in different ways, which have different advantages 
and disadvantages. For example, 

1 1 1 



4 — x 4 1 - x/A 

oo 

1 /X\ n 

— ■ 2_. ( t ) converges for all \x\ < A. 



4 ^ V4 

n=0 



Notice that the second series converges for \x\ < A, whereas the first converges 
only for \x — 3| < 1, which isn't nearly as good. 
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5.5. POWER SERIES 



5.5.2 Convergence of Power Series 

Theorem 5.5.1. Given a power series X^^Lo c n(x — &)" . there are exactly three 
possibilities: 

1. The series converges only when x = a. 

2. The series converges for all x. 

3. There is an R > ( called the "radius of convergence") such that ^2^ =0 c n (x— 
a) n converges for \x — a\ < R and diverges for \x — a\ > R. 

Example 5.5.4. For the power series X^Lo a; ™' ^ e ra dius R of convergence is 
1. 

Definition 5.5.3 (Radius of Convergence). As mentioned in the theorem, R is 
called the radius of convergence. 

If the series converges only at x = a, we say R = 0, and if the series converges 
everywhere we say that R = oo. 

The interval of convergence is the set of x for which the series converges. It 
will be one of the following: 

(a-R,a + R), {a-R,a + R), (a — R,a + R], [a-R,a + R] 

The point being that the statement of the theorem only asserts something about 
convergence of the series on the open interval (a — R, a + R). What happens 
at the endpoints of the interval is not specified by the theorem; you can only 
figure it out by looking explicitly at a given series. 

Theorem 5.5.2. If Y^Lo c n( x ~ a )™ ^as radius of convergence R > 0, then 
f( x ) = S^Lo c n( x ~ a ) n * s differentiable on (a — R, a + R), and 



and both the derivative and integral have the same radius of convergence as f . 

Example 5.5.5. Find a power series representation for f{x) = tan _1 (x). No- 
tice that 



which has radius of convergence R = 1, since the above series is valid when 
| — x | < 1, i.e., \x\ < 1. Next integrating, we find that 



1. f'(x) = ■ c n (x - a) 



n=l 
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5.6. TAYLOR SERIES 



for some constant c. To find the constant, compute c = /(0) = tan (0) = 0. 
We conclude that 



tan- 1 (x) = £(-l)" 



2n + 1 

71=0 

Example 5.5.6. We will see later that the function fix) = e~~ x has power 
series 

e~* 2 = l-x 2 + 1 -x i - 1 -x e + --- . 
2 6 

Hence 

-t 2 i 1 1 *i 1 7 

This despite the fact that the antiderivative of e~ x is not an elementary func- 
tion. 



5.6 Taylor Series 

Example 5.6.1. Suppose we have a degree-3 (cubic) polynomial p and we know 
that p(0) = 4, p'(0) = 3, p"{0) = 4, and p"'(0) = 6. Can we determine p? 
Answer: Yes! We have 

p(x) = a + bx + ex 2 + dx 3 

p'{x) = b + 2cx + 3dx 2 
p (x) = 2c + 6dx 
p"'(x) = 6d 

From what we mentioned above, we have: 

a = p(0) = 4 
6 = p'(0) = 3 
p"(0) 



2 

= 1 
6 

Thus 



2 

p"'(0) 



p[x) = 4 + 3x + 2.t 2 + x 3 



Amazingly, we can use the idea of Example 5.6.1 to compute power series 
expansions of functions. E.g., we will show below that 

< v ,. 

n=0 
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5.6. TAYLOR SERIES 



Consider a general power series 

oo 

f(x) = ^2 °n( x ~ a Y l = c + ci(x - a) + c 2 {x — a) 2 H 



n=0 

We have 



CO 


= /(«) 


Cl 


= /» 




_ /"(«) 


C2 


2 



C " ~ n! ' 
where for the last equality we use that 

f^\x)=n\c n + {x-a){- ■• + ■■■) 

Remark 5.6.1. The definition of 0! is 1 (i£'s i/ie empty product). The empty 
sum is and the empty product is 1. 

Theorem 5.6.1 (Taylor Series). If f(x) is a function that equals a power series 
centered about a, then that power series expansion is 

/w-f^(.-.r 



n=0 



= f(a) + f'(a)(x -a) + ^(x - a) 2 + ■ ■ ■ 

Remark 5.6.2. WARNING: There are functions that have all derivatives de- 
fined, but do not equal their Taylor expansion. E.g., f(x) = e~ x l x forx ^ and 
/(0) = 0. It's Taylor expansion is the series (which converges everywhere), 
but it is not the function. 

Definition 5.6.1 (Maclaurin Scries). A Maclaurin series is just a Taylor series 
with a = 0. I will not use the term "Maclaurin series" ever again (it's common 
in textbooks). 

Example 5.6.2. Find the Taylor series for f{x) = e x about a = 0. We have 
fW(x) = e x . Thus / (n) (0) = 1 for all n. Hence 

—x n = l + x+ — + — + ■■■ 
n! 2 6 
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5.6. TAYLOR SERIES 



What is the radius of convergence? Use the ratio test: 



lim 



1 x n+l 



("+!) 

1 



lim 7 — TTTT^I 



lim 



= 0, for any fixed x. 



n^oo n + 1 

TTms i/ie radius of convergence is oo . 

Example 5.6.3. Find the Taylor series of fix) = sin(a;) about x = %c\ We 



m = E 



n=0 



To do i/iis we have to puzzle out a pattern: 

f(x) = sm(x) 
f'(x) = cos(x) 
f"(x) = -sm(x) 
f"'(x) = -co S (x) 
f^\x) = sin{x) 

First notice how the signs behave. For n = 2m even, 

/(»)(!) = f 2m \x) = (-1)"/ 2 sm(x) 
and for n = 2m + 1 odd, 

/(»)(a:) = f( 2m+1 \x) = (-l) m cos(.T) = (_i)(«-D/ 2 CO s(a;) 
For n = 2m even we have 

/< n) ^/2) = / (2m) (!) =(-l) m . 
and /or n = 2??i + 1 odd we have 

/W(tt/2) = /( 2m+1 ) (|) = (-l) m cos(7r/2) = 0. 

Finally, 



sin 



n=0 



E 



m=0 



m 
(-1)' 



7TN 2m 

- (2m)! K X ~ 2, 



1 Evidently this expansion was first found in India by Madhava of Sangamagrama (1350- 
1425). 
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Next we use the ratio test to compute the radius of convergence. We have 

(2m)! 



lim 

m — >oo 



(_l)m+l 


{'-I 


2(m+l) 


(2(m + l))! 




(-l) m / 
(2m)! ' 




2m 





lim 

m^oo (2m + 2)! 



lim 



(*-!)■ 



m^oo (2m + 2) (2m + 1) 
which converges for each x. Hence R = oo. 

Example 5.6.4. find £/ie Taylor series for cos(a;) about a = 0. M^e /iave 
cos(x) = sin fx + f ) • XTws /rom Example 5. 6. 3\ (with infinite radius of conver- 
gence) and that the Taylor expansion is unique, we have 



s(x) = sin (x + |J 



E 

oo 

E 



(-i)" 

(2n)! 



(2n)! 
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5.7 Applications of Taylor Series 



This section is about an example in the theory of relativity. Let to be the 
(relativistic) mass of an object and too be the mass at rest (rest mass) of the 
object. Let v be the velocity of the object relative to the observer, and let c be 
the speed of light. These three quantities are related as follows: 



to 



(relativistic) mass 



1 - 



The total energy of the object is mc 2 : 



E = mc 2 \\ 



In relativity we define the kinetic energy to be 

K = mc 2 — tooC 2 . 

What? Isn't the kinetic energy ^mov 2 ? 
Notice that 



(5.4) 



mc 2 — TOqC 



TTIqC 



TOqC 2 = TOqC 2 



1 

r 



v 2 

l -V 2 



Let 



f(x) = (l-x) 5 -l 
Let's compute the Taylor series of /. We have 

f{x) = {l-x)- 1 * - I 

f(x) = \{i-x)-i 



f"(x) 



1 3 

2 ' 2 



(1-x)- 



Thus 



Hence 



/ (n) (0) 



1 • 3 • 5 • • • (2n - 1) 



2" 



/(*) = E 



/ ( ' l) (o)^„ 



~ 1.3-5---(2n-l) 
2" • n\ 

n=l 

1 3 2 5 , 35 
-x H — x H a; H x 

2 8 16 128 



150 
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We now use this to analyze the kinetic energy ( [5.4) : 



mc 2 — jtiqc 2 = toqc 2 • / ( 



m c 



( 1 w 2 3 ^ 
2 ' <? + 8 ' c? 



1 2 2 /3« 2 



And we can ignore the higher order terms if \ is small. But how small is 

2 ... 

; 'small" enough, given that \ appears in an infinite sum? 



5.7.1 Estimation of Taylor Series 

Suppose 



Write 



We call 



^(x) :=/(*)- V ^— P(z-a)" 



n=0 



JV 



= 2^ j — - a ) 



n=0 



the TVth degree Taylor polynomial. Notice that 



lim T N (x) = f(x) 



if and only if 



N 



lim Rn{x) = 0. 



We would like to estimate f(x) with Tn(x). We need an estimate for R^(x). 
Theorem 5.7.1 (Taylor's theorem). If |/ (w+1) 0)l < M for \x - a\ < d, then 

\R N (x)\ < ^^y \x - a\ N+1 for \x a\ < d. 

For example, if N = 0, this says that 

\R(x)\ = \f(x)-f(a)\<M\x-a\, 



< M, 
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which should look familiar from a previous class (Mean Value Theorem). 
Applications: 

1. One can use Theorem |5.7.1| to prove that functions converge to their Taylor 
series. 

2. Returning to the relativity example above, we apply Taylor's theorem 
with N = 1 and a = 0. With x = —v 2 /c 2 and M any number such that 
\f"(x)\ < M, we have 

\Ri(*)\ < y J- 
For example, if we assume that \v\ < lOOm/s we use 

|/"0)| < ^(1 - 100 2 /c 2 )~ 5/2 = M. 

Using c = 3x 10 8 m/s, we get 

\Ri{x)\ <4.17-10~ 10 -m . 

Thus for v < lOOm/s ~ 225mph, then the error in throwing away rela- 
tivistic factors is 10 _10 too- This is like 200 feet out of the distance to the 
sun (93 million miles). So relativistic and Newtonian kinetic energies are 
almost the same for reasonable speeds. 
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Chapter 6 

Some Differential Equations 



This chapter is an introduction to differential equations, a major field in applied 
and theoretical mathematics and a very useful one for engineers, scientists, 
and others who study changing phenomena. The physical laws of motion and 
heat and electricity can be written as differential equations. The growth of a 
population, the changing gene frequencies in that population, and the spread 
of a disease can be described by differential equations. Economic and social 
models use differential equations, and the earliest examples of "chaos" came 
from studying differential equations used for modeling atmospheric behavior. 
Some scientists even say that the main purpose of a calculus course should be 
to teach people to understand and solve differential equations. 
Differential Equations 

Algebraic equations contain constants and variables, and the solutions of an 
algebraic equation are typically numbers. For example, x = 3 and x = —2 
are solutions of the algebraic equation x 2 = x + 6. Differential equations con- 
tain derivatives or differentials of functions. Solutions of differential equations 
are functions. The differential equation y' = 3x 2 has infinitely many solutions, 
and two of those solutions are the functions y = x 3 + 2 and y = x 3 — 4. You 
have already solved lots of differential equations: every time you found an an- 
tiderivative of a function f(x), you solved the differential equation y' = f(x) 
to get a solution y. You have also used differential equations in applications. 
Areas, volumes, work, and motion problems all involved integration and finding 
antiderivatives so they all used differential equations. The differential equation 
y' = f(x), however, is just the beginning. Other applications generate different 
differential equations. 

Checking Solutions of Differential Equations 

Whether a differential equation is easy or difficult to solve, it is important to be 
able to check that a possible solution really satisfies the differential equation. A 
possible solution of an algebraic equation can be checked by putting the solution 
into the equation to see if the result is true: x — 3 is a solution of 5x + 1 = 16 
since 5(3) + 1 = 16 is true. Similarly, a solution of a differential equation can 
be checked by substituting the function and the appropriate derivatives into the 
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6.1. SEPARABLE EQUATIONS 



equation to see if the result is true: y = x 2 is a solution of xy' = 2y since y' = 2x 
and x(2x) = 2(x 2 ) is true. 
2 

Example 6.0.1. For every value of C, the function y = Cx 2 is a solution of 
xy' = 2y. Find the value of C so that y(5) = 50. 

Solution: Substituting the initial condition x — 5 and y = 50 into the solution 
y = Cx 2 , we have that 50 = C5 2 so C = 50/25 = 2. The function y = 2x 2 
satisfies both the differential equation and the initial condition. 

6.1 Separable Equations 

A separable differential equation is a first order differential equation that can be 
written in the form 

= f( x ) 
dx h(y) 

These can be solved by integration, by noting that 

Hv)dy = f(x)dx, 

hence 

J h(y)dy = J f(x)dx + C. 

This latter equation defines y implicitly as a function of x (we have added 
a "+C" onto one side just to emphasize that you only need one constant of 
integration for the solution), and in some cases it is possible to explicitly solve 
for y as a function of x. 

6.2 Logistic Equation 

The logistics equation is a differential equation that models population growth. 
Often in practice a differential equation models some physical situation, and 
you should "read it" as doing so. 
Exponential growth: 

1 dP - k 

This says that the "relative (percentage) growth rate" is constant. As we saw 
before, the solutions are 

P(t) = Po ■ e kt . 

Note that this model only works for a little while. In everyday life the growth 
couldn't actually continue at this rate indefinitely. This exponential growth 
model ignores limitations on resources, disease, etc. Perhaps there is a better 
model? 

Over time we expect the growth rate should level off, i.e., decrease to 0. What 
about 
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6.2. LOGISTIC EQUATION 



— = *(l-£), (6-1) 
P dt \ KJ' y ' 

where K is some large constant called the carrying capacity, which is much big- 
ger than P = P(t) at time 0. The carrying capacity is the maximum population 
that the environment can support. Note that if P > K, then dP/dt < so the 
population declines. The differential equation ( j6.1J ) is called the logistic model 
(or logistic differential equation). There are, of course, other models one could 
use, e.g., the Gompertz equation. 

First question: are there any equilibrium solutions to ( |6.1|) , i.e., solutions with 
dP/dt = 0, i.e., constant solutions? In order that dP/dt = then = k (l — ^) , 
so the two equilibrium solutions are P(t) = and P(t) = K. 

The logistic differential equation ( ]6.1| ) is separable, so you can separate the 
variables with one variable on one side of the equality and one on the other. 
This means we can easily solve the equation by integrating. We rewrite the 
equation as 

dP KP(P-K). 



Now separate: 



and integrate both sides 



dt K 

KdP 
P{P - K) 



KdP 



P{P - K) 



= —k ■ dt, 



-k-dt = -kt + C. 



On the left side we get 



KdP f ( 1 1 

P(P - K) ~ J \P-K ~ P 



dP = \n\P-K\ - lnlPl 



Thus 



\n\K -P\ -ln|P| = -kt + c, 

so 

In | (If - P)/P| = -kt + c. 
Now exponentiate both sides: 

{K-P)/P = e~ kt+c = Aer kt , where A = e c 



Thus 



so 



K = P(l + Ae- kt ) 



K 
1+Ae- 
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6.2. LOGISTIC EQUATION 



Note that A = also makes sense and gives an equilibrium solution. In general 
we have lim^oo P(t) = K . In any particular case we can determine A as a 
function of P = P(0) by using that 

K K K-P Q 

P(0 = so A= 1 = -. 

1 + ^4 P P 
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Chapter 7 

Appendix: Integral tables 



Trigonometry 

Trigonometric Functions 

Tl. sin 2 x + cos 2 x = 1 

T2. tan 2 x + 1 = sec 2 x 

T3. cot 2 x + l = esc 2 X 

T4. sin(x ± j/) = sin x cos y ± cos x sin y 

T5. cos(x ± y) = cos a; cosy =F sin a; siny 

tan a; ± tan y 

To. tanfx ±y) — 

1 =p tan x tan y 

„ . x . sin a; 

T7. tan -) = 

y 2' 1 + cos a; 

T8. sin(2x) = 2 sin a; cos a; 

T9. cos(2a;) = cos 2 x — sin 2 x 
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sin 2 x = ^(l — cos(2x)) 
cos 2 x = V2(l + cos(2x)) 
sin a; siny = Y2(cos(x — y) — cos(x + y)) 
cosx cosy = l /i(cas{x — y) + cos(x + y)) 
sinx cos y = y2(sin(x — y) + sin(x + y)) 
c\ cos(ajt) + C2 sin(u>t) = A sin(wt + 0), 



where A = yl c\ + c 2 , , 4> — 2 arctan 



Hyperbolic Functions 



e x + e x 

coshx = 



c 2 + A 



2 

sinh x 



cosh 2 x — sinh 2 x = 1 



tanh 2 x + sech 2 x = 1 



coth 2 x — csch 2 x = 1 



sinh(x ± y) = sinhx coshy ± coshx sinhy 
cosh(x ± y) = cosh x cosh y ± sinh x sinh y 
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tanh x ± tanh y 

tanh(ai ± y) 



1 ± tanh x tanh y 
sinh(2a;) = 2 sinh a; cosh a; 

cosh ( 2a; ) = cosh 2 x + sinh 2 x 

sinh a; sinh y = 1 /2(cosh(a; + y) ~ cosh(x — y)) 

cosh a; coshy = 1 /2{cos\\{x + y) + cosh(x — y)) 

sinh a: coshy = 1 /2(sinh(a: + y) + sinh(x — y)) 

Power Series 



00 n 2 3 

x sr^ x , x x 

^ n\ 2! 3! 



n=0 



sina; = E(- 1 ) n ( 2n + l)l = X ~ 3! + 5! " 



n=0 



cosx — y^(— 1)" — — - 7 = 1 — — + — — ■•■ , — OO < X < OO 



(2n)\ 2! 4! 

n=0 y 1 



,:■'< 



tlmX = X+ Y + T5 x5 + 315 x7 + ---> -2 <X< 2 



7T 7T 



1 + X + x 2 + X 3 H , -1< X < 1 



n=0 

00 a; 2rl+1 x 3 x 5 

Sinh - = E(^TT)! = - T+ 3! + 5! + '"' -°°<*<°° 



n=0 

n=0 ^ ; ' 

X 3 2 - 17 7 7T 7T 

tanh a: = a; 1 x x + • • • , < x < — 

3 15 315 2 2 
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P9. f(x) = f(a) + f(a)(x - a ) + f -^{x - af + f -^-{x - af 



P10. Taylor Scries with remainder: 



m ^Zto^ix-ar + RN+iix), where 

Rn+i(x) = ^ N+1 y ' > {x — a) N+1 for some £ between a and 



Table of Integrals 

A constant of integration should be added to each formula. The letters a, 
b, m, and n denote constants; u and v denote functions of an independent 
variable such as x. 

Standard Integrals 



11. / u n du = , n± -1 

J n + 1 

12. / — = In « 

J u 

13. J e u du = e" 

14. / a u du = —, a > 

J m a 

15. J cosudu = sinu 

16. / sin udu = — cos u 

17. y sec 2 udu = tan u 

18. / esc 2 udu = — cot u 
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19. 



/ 



sec u tan udu = sec i 



110. 



esc u cot udu = — esc \ 



111. 



tan u du = — In I cos \ 



112. 



cot u du = In I sin i 



113. 



114. 



115. 



116. 



sec u du = In I sec w + tan u| 



/ 



esc -u = In I esc u — cot i 



1 

— = - arctan 

a z + vr a va< 



du 

\J a 2 — u 2 



arcsm 



(;) 



117. 



y udv = uv — J v 



Integrals involving au + b 



118. 



119. 



,v» , (au + 6)" +1 

(au + 6)"du=^ — , n^-1 

(n + l)a 

1, , . 

m \ au + b\ 



au + b a 



120. 



udu u b 

— = 2 In \au + b\ 

au + b a a z 



121. 



122. 



udu b 1 

(au + b) 2 a 2 (au + b) a 



H — - In \au + b\ 



du 1 

~i TV = T nl 

u(au + b) b 



u 

au + b 
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123. 



124. 



125. 



126. 



ly/au + b du 



u 2 \J au + bdu 



\J au-\-b 
2 



105a 3 
2 du 2 



y/au + b 15a 3 



2(3au - 26) 



15a 2 

u du 2(au — 2b) 



(au + bf /2 



3a 2 



Vau+~b 



- \2abu + 15aV) (au + bf /2 
Aabu + 3a 2 tt 2 ) V au + b 



Integrals involving u 2 ± a 2 



127. 



128. 



129. 



130. 



131. 



du 1 

= 7T- In 



u 2 — a 2 2a 
udu 



u 2 ± a 2 
u 2 du 



u + a 
i/ 2 lnU 2 ±a 2 



a i 

— a z 2 



u + a 



u au (u 
u — a arctan 



u z + a 2 



du 



u(u 2 ± a 2 ) 2a 2 



it 2 ± a 2 



Integrals involving y/u 2 ± a 2 



132. 



133. 



134. 



135. 



-\Af 2 ± a 2 



= \/ u 2 ± a 2 



/ U ^ 2 ± a 2 = ^ (« 2 ± a 2 ) 37 ' 



(/(/ 



Vu 2 ± a 2 



In 



+ \Ju 2 ±i 



Vu 2 ± a 2 2 



= — V u 2 ± a 2 =F — In u + \At 2 ± 
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136. 



137. 



138. 



139. 



/ 



du 1 

= - In 



u\J u 2 + a 2 a 

du 1 



a + \J u 2 + a 2 



— -j^=^= = — arcsec 
uy/u 2 - a 2 a 



du 



\Ju 2 ± a 2 



i 2 y/u 2 ± a 2 a 2 u 

,2 



/U CI 
\/u 2 ±a 2 du = — V u 2 ± a 2 ± — In u + \/u 2 ±a 2 



140. 



141. 



142. 



143. 



1 \Ju 2 ±a 2 du — -(u 2 ±a 2 ) =F -g- \/ u 2 ± a 2 - — In u + \/ u 2 ± , 



f Vu 2 +a 2 / 

I du = v u + a a In 

J u 

I 



+ V 'u 2 + a 2 



V u 2 — a 2 



du = \f' 



u 2 — a 2 — a arcsec 



Vu 2 ^* 



du 



\/u 2 ± a 2 



In 





u + y/ u 2 ± a 2 



Integrals involving \J a 2 — u 2 



144. 



J \J a 2 — u 2 du = — \/ a 2 — u 2 + — arcsin ^ 



145. 



146. 



147. 



148. 



149. 



u du 



V 'a 2 — u 2 



■\J a 2 — u 2 



i\J a 2 — u 2 du = — - (a 2 — u 2 ) 3 ^ 2 
o 



/\^Q 2 — U 2 / 
du = v a 2 — u 2 — a In 
it 



+ \/a 2 — u' 2 



du 



u\J a 2 — u 2 a 



In 



+ Va 2 



! V 'a 2 -u 2 du = -~ (a 2 - u 2 ) 3/2 + \j a 2 - u 2 + °— 
4 8 8 



arcsm 
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152. 



161. 



162. 



V a 2 — u 2 v a 



2 — -1/2 



U^ / U 

150. / s an = arcsin — 

J u 2 u \a 

f u 2 du u /—„ a 2 . (u\ 

151. / . = — — v a — n • + it arcsin — 
J Va 2 - u 2 2 V 2 \aJ 



du V a 2 — u 2 



u 2 \/a 2 - u 2 a 2 u 
Integrals involving trigonometric functions 

u sin(2au) 



153. / sin (ail) du = — 



2 4a 



U sin(2au) 



154. / cos (au) du = - 

J ' 2 4a 

/ 3 If cos 3 (au) 
sin (au) au = — cos(au) 
a V 3 

f 3/ \ j 1 /" • / \ sin 3 (au) 
156. / cos (au J au = - I sm(au) 

/u 1 
sin (au) cos 2 (au) du = — — — — sin(4au) 
8 32a 

158. / tan 2 (au) du = - tan(au) — u 



159. J cot 2 (au) du = — cot (au) — u 

/I i 
sec 3 (au) du = — sec(au) tan(au) H In I sec(au) + tan(au) 
2a 2a 



f 1 1 

/ csc 3 (au) du = —— csc(au) cot(au) + — In | csc(au) — cot(au) 

f usin(au)du = — (sin(au) — aucos(au)) 
J a 2 



163. / wcos(au) du = — (cos(au) + ausin(au)) 
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J u 2 sin(ait) du = (2au sin(au) — ( a 2 u 2 — 2 ) cos(au)) 

/u 2 cos(au) du — — r (2au cosiau) + ( a 2 u 2 — 2 ) sin(au)) 

f . . sin(a — o)u sin(a + 6)u 2 , 2 

y S1 n(au)sm(H^=^36) 2^T&P ° * 6 

/" . sin(a - 6)it , sin(a + b)u 2 2 

/ cos aw cos(6u du = — -— + — — — — , a ^ b 

J Z(a — b) 2{a + o) 

/. , . s , cos(a — 6)« cosfa + o)it 2 , 

smWcosHdM = __ _____ a ^ 6 

/ sin" itdu = — — sin" -1 ucosu + — — - f sin" -2 udu 
J n n J 

Integrals involving hyperbolic functions 

/sinh(au) du = - cosh(cra) 
a 

/I u 
sinh (au) du = — sinh(2au) — — 
4a _ 

f 1 

/ cosh(au) du = — sinh(au) 
J a 

f u 1 

/ cosh (au) du = — + — sinh(2em) 
J v ' 2 4a 

f . ,, , , /, n , cosh ((a + b) u) cosh ((a- 6) u) 

/ sinh(au) cosh(6u) = ^ H ^ — ^ — - 

J 2(a + 6) 2(a-&) 

f sinh(ait) cosh(au) du = — cosh(2au) 
J 4a 

J tanhudu = ln(coshu) 
J sechudu = arctan(sinhit) = 2axctan(e") 
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Integrals involving exponential functions 



ue au 



du = — K-(au — 1) 



/ au 
u 2 e au du = -3- ( a 2 u 2 - 2au + 2 ) 



e au du = -u n e au - - I u n - l e au du 
a a 



e au siaipu) du = (asinfbu) — bcos(bu)) 

a 2 + tr 



e au cos(bu) du = — —r (acos(6u) + bsintbu)) 
a z + tr 



Integrals involving inverse trigonometric functions 



csin ^ — ^ du = uarcsin ^ — ^ + \fc 



1 



)s ^ — ^ du = uarccos ^— ^ — \fo 



arctan ( — J du = u arctan ^— j — — In + u J 



Integrals involving inverse hyperbolic functions 



/ arcsinh ^ — ^ du — u arcsinh ^ — ^ — \f- 



u 1 + a 2 



arccosh^ — J du — u arccosh (^) — \Ju 2 — a 2 arccosh (^) > 
= u arccosh (r) + Vu 2 — a 2 arccosh (^) < 
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J arctanh ^— ^ du = it arctanh ^— ^ + — In (a 2 — u 2 ) 



Integrals involving logarithm functions 



J \nudu = u(\nu — 1) 



/ u n In u du = u 



Wallis' Formulas 



In u 



71 + 1 (n+1) 2 



n ^ -1 



/ sin" 1 xdx = J^ 2 cos" 1 x dx 
Jo 



(m-l)(m-3)...(20ri) , 
m(m-2)...(3 0r2) ^' 



where fc = 1 if m is odd and k = n/2 if m is even. 



tt/2 

sin" 1 x cos" x = 

o 

(m - l)(m - 3) . . . (2 or 1) (re- l)(n- 3) . . . (2or 1) 
(m + n)(m + 7i — 2) . . . (2 or 1) 

where k = tt/2 if both m and n are even and fc = 1 otherwise. 
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